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ABSTRACT 

An  echosounder  is  used  to  probe  various  atmospheric 
parameters.   An  acoustic  wave  is  transmitted  into  the 
atmosphere  and  information  deduced  from  the  backscattered 
energy. 

This  thesis  seeks  to  understand  the  range  limitations  of 
the  echosounder  and  to  explore  methods  to  quantify 
atmospheric  turbulence  parameters  at  a  given  range.   The 
propagation  of  the  acoustic  energy,  including  the  effects  of 
excess  attenuation,  are  modeled  to  predict  the  performance 
of  an  echosounder  when  various  parameters  are  changed.   The 
electronics  of  an  existing  echosounder  are  investigated  to 
understand  inherent  or  design  limitations. 
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I.  INTRODUCTION 

Many  atmospheric  turbulence-dependent  processes  take 
place  in  the  troposphere,  the  lowest  ten  to  twenty 
kilometers  above  the  earth's  surface.   Acoustic  energy 
interacts  with  irregularities  in  the  atmosphere  more 
strongly  than  electromagnetic  waves  and  is  potentially  a 
better  probe  for  these  irregularities.   Acoustic  echoes  have 
been  used  to  detect: 

1)  wind  speed  and  direction  profiles, 

2)  humidity  profiles, 

3)  strength  and  location  of  temperature  inversions, 

4)  temperature  inhomogenei ties , 

5)  mechanical  turbulence. 

The  acoustic  sounder,  also  known  by  the  names 
echosounder,  echosonde,  sodar,  and  acoustic  radar,  transmits 
a  pulse  of  acoustic  energy  into  the  atmosphere.   The  various 
atmospheric  parameters  can  be  determined  based  on  the 
intensity  and  frequency  of  the  scattered  energy.   The  range 
to  the  scattering  center  is  determined  from  the  elapsed  time 
between  transmission  and  return  of  the  scattered  energy. 

In  this  thesis  we  are  attempting  to  improve  our 
understanding  of  the  fundamental  range  limitations  of  the 
acoustic  sounder.   In  general,  the  acoustic  sounder  suffers 


from  multiple  scattering.   Turbulence  at  shorter  ranges 
alters  the  phase  front  of  the  propagating  pulse  and  reduces 
the  magnitude  of  the  return  from  longer  ranges. 

We  are  also  exploring  the  use  of  return  from  shorter 
ranges,  along  with  theory  for  the  degradation  of  the  energy 
along  the  path,  to  compensate  for  short  range  degradation  in 
a  boot-strap  fashion  and  quantify  the  return  from  longer 
ranges.   Acoustic  sounders  at  present  indicate  the  presence 
of  inhomogeneities  at  a  specific  altitude  but  the  magnitude 
of  the  inhomogeneities  may  have  errors  of  a  factor  of  four 
or  more . 

The  following  factors  affect  the  range  of  the 
echosounder : 

1)   Atmospheric: 

Pressure , 

Temperature , 


Temperature  Structure  Parameter  CT  , 

2 
Velocity  Structure  Parameter  Cv   , 

Water-vapor  Pressure,  and 

Ambient  Noise. 


2)  Echosounder: 


Antenna  Aperture  Factor, 
Antenna  Diameter, 
Efficiency, 
Frequency, 
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e)  Power  Transmitted,  and 

f)  Pulse  Length. 
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II.  BACKGROUND 

Acoustic  energy  propagates  in  the  atmosphere  as  a 
longitudinal  wave  of  pressure  variations.   The  energy  is 
scattered,  attenuated,  and  refracted.   Energy  which  is 
scattered  constitutes  the  returned  signal.   The  attenuation 
decreases  the  energy  ensonifying  a  given  volume  and 
decreases  the  returned  signal.   Energy  that  is  scattered 
through  small  angles  or  that  is  refracted  results  in  further 
degradation  or  excess  attenuation  of  the  signal. 

A.   SCATTERING 

The  echosounder  transmits  a  pulse  of  acoustic  energy 
that  is  scattered  by  temperature  and  velocity 
inhomogeneities .   Information  about  the  inhomogeneities  is 
then  based  on  elapsed  time  between  transmission  and  receipt 
of  the  return  signal,  the  strength  of  the  returned  signal, 
the  equations  for  scattering,  and  the  Doppler  shift. 
Variations  in  the  propagation  velocity  of  the  wave,  which 
are  a  function  of  the  temperature  and  velocity  variations  in 
the  air  mass,  produce  the  scattering.   In  addition  to 
scattering,  the  temperature  variations  cause  refraction 
while  the  velocity  turbulence  causes  both  a  shift  in  phase 
and  direction  of  propagation  [Ref.  l:p.  60]. 
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Random  temperature  or  wind  structure  caused  by 
turbulence  and  uniform  gradients  of  temperature  or  velocity 
contribute  to   the  scattering.   The  gradient  must  change  in 
a  scale  size  that  is  comparable  to  or  smaller  than  the 
acoustic  wavelength  in  order  to  contribute  to  the 
scattering.   It  appears  that  the  effects  of  uniform 
gradients  are  limited  to  beam  bending  and  possible  specular 
reflection  for  the  acoustic  frequencies  of  interest  here. 
This  refractive  structure  of  the  troposphere  and 
stratosphere  would  be  better  probed  with  what  has  come  to  be 
called  infrasound,  sound  of  frequencies  below  twenty  Hertz. 
[Ref.  l:p.  52] 

The  turbulent  atmospheric  temperature  and  velocity 

fluctuations,  at  high  wave  number,  follow  the  Kolmogorov 

-5/3 
K   '   power  spectral  density  law  (in  one  dimension).   The 

fluctuations  are  correlated  spatially  on  the  order  of  one 
centimeter  to  tens  of  meters.   The  nonzero  correlation 
lengths  and  the  declining  power  spectral  density  imply  a 
different  scattering  process  than  for  random  point 
scatterers,  even  though  the  turbulence  is  a  stochastic 
process  [Ref.  l:p  61]. 

An  expression  for  the  power  scattered  from  a  unit  volume 
per  unit  incident  flux  into  a  unit  solid  angle  is  [Ref.  2:p 
84] 

(J    =,055  cos2o  (£y2  cos2(i)  +  .13  C±2 ) ( sin(|)  )"11/3  , 
y/3         c2       2  T2       2 

IO 


where 

is  the  wavelength  of  the  transmitted  sound, 
T  is  the  average  atmospheric  temperature, 
0  is  the  angle  of  scattering  from  the  direction  of 

propagation, 

2 
CV"  is  the  velocity  structure  parameter, 


.2=  „VI..S      _  „/..i-»v2 


<VX3L  T/y(x+r)>^  ' 


r" 

and  Ct"1  is  the  temperature  structure  parameter, 
Ct2=  <T(x)  -  T(x+r)>2  . 

This  scattering  equation  indicates  that  backscatter 
(&=tt)    is  only  a  function  of  the  temperature  structure.   The 
velocity  at  any  point  is  the  vector  sum  of  the  phase 
velocity  of  the  sound  and  the  particle  velocity  of  the 
turbulence.   This  vector  sum  is  always  in  the  forward 
direction  because  the  particle  velocity  is  always  less  than 
the  speed  of  sound  c.   Therefore  the  backscatter  is  only  a 
function  of  the  temperature  structure  [Ref .  l:p  60].   The 
vertical  turbulent  velocity  does  cause  a  Doppler  shift  of 
the  frequency  of  the  backscattered  radiation. 

Scattering  over  a  region  with  correlated  scattering 
centers  produces  constructive  and  destructive  interference. 
The  backscattered  waves  are  partially  coherent.   This 
results  in  much  greater  intensities  than  would  be  received 
from  incoherent  scattering  such  as  Rayleigh  scattering  [Ref. 


l:p  61].   Also  the  scattering  with  the  interference,  over 
the  nonzero  coherence  length,  acts  like  an  array  of 
scattering  centers  in  a  regular  crystal  lattice  of  spacing 
L.   The  Bragg  condition 

L=  ^  /sin(f)/2), 
is  satisfied.   The  dominant  scattering  is  for  scale  sizes  of 

X/2.   [Ref.  l:p  61] 

B .   ATTENUATION 

The  atmosphere  absorbs  some  acoustic  energy  that 

propagates  through  the  atmosphere  and  reradiates  this  energy 

at  different  frequencies.   Quantitatively  the  power  lost 

(P-,)  over  a  path  length  1  is  given  by 

P-j_=  e    , 

where    is  the  attenuation  coefficient  in  Nepers  per  meter. 

Historically  the  attenuation  coefficient  has  been  considered 

to  be  the  sum  of  several  terms; 

c<   =  classical  viscous  losses, 
t?v  cl 

x  r=  molecular  rotation  losses, 

^vib=  molecular  vibrational  losses,  N2   and  C^* 

Classical  and  rotational  loses  are  negligible  below 
about  three  kilohertz,  the  region  of  interest  for 
echosounder  operation  [Ref.  l:p  54  and  Ref.  3:p.  18-2], 
consequently  vibration  of  N2  and  02  produce  most  of  the 
attenuation. 
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Rotational  and  vibrational  losses  are  referred  to  as 
relaxation  processes.   The  acoustic  energy  excites  internal 
energy  modes  of  the  N2  and  02  molecules.   The  rate  of 
collisions  with  water  vapor  determines  the  rate  and 
efficiency  of  conversion  of  the  energy  into  translational 
energy  ( heat) . 

The  phase  is  shifted  due  to  the  relaxation  processus. 
This  is  one  of  the  reasons  some  of  the  microwave  radar  pulse 
compression  techniques  cannot  be  used  in  echosounders .   A 
number  of  the  pulse  compression  techniques  rely  on  the  phase 
not  changing  during  propagation. 

The  dependence  of  the  attenuation  on  the  water  vapor 
pressure  is  believed  to  be  due  to  a  resonance  process 

between  the  lowest  vibrational  states  of  the  O2  and  N2 
molecules  with  the  water  molecules.   Henderson  and  Hertfeld 
in  Reference  4  [p. 986]  state  that  the  lowest  vibrational 
states  of  C>2  and   water  vapor  are  only   thirty-nine 
wavenumbers  (56  K)  apart  at  1600  cm"  .   For  this  reason 
oxygen  was  thought  to  be  primarily  responsible  for  the 
humidity  dependent  absorption  of  sound.   Henderson  and 
Herzfeld  [Ref.  4],  and  many  others,  assumed  nitrogen  was  an 
inactive  dilutent  having  no  effect  on  the  relaxation 
processes  [Ref.  4:p.  987].   Unfortunately  the  theory  did  not 
agree  with  the  data  for  low  frequencies  and  high 
humidities . 
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Theory  incorporating  the  relaxation  processes  involving 
nitrogen  and  water  at  low  frequencies  and  relative 
humidities  greater  than  twenty-five  percent  bring  the  theory 
into  agreement  with  the  data.   At  high  frequencies  and  low 
relative  humidities,  theory  and  data  match  well  with  oxygen 
making  the  main  contribution  and  nitrogen  acting  as  an  inert 
dilutent.   At  low  frequencies  and  high  humidities,  nitrogen 
seems  to  make  the  main  contribution  and  oxygen  acts  as  an 
inert  dilutent.   For  the  range  of  relative  humidities  in 
which  many  people  live  and  over  a  good  part  of  the  audible 
frequency  range  it  appears  nitrogen  is  the  major  contributer 
to  the  relaxation  processes  [Ref.  5:p.  165].   This  is  in 
contrast  to  the  previous  assumption  that  oxygen  was 
responsible  [Ref.  6:p.  604], 

An  expression   for  attenuation  due  to  the  molecular: 

absorption  is  given  by  [Ref.  7:p.  34] , 

<*=  °<max  ^018*f  \2  +fl*  {£/£m)2\  2       L/2  , 
304. 8^    fmj   Ll  +(f/fm)  /   „ 

in  dB/m, 

where 

f  is  the  frequency  of  the  transmitted  pulse, 

f   is  the  Napier  frequency, 

fm    =(10+6600    h+4400   h2)    (      P    \   /         519       ^  °'8    , 

1014      V1.3T+492 


o^  =  .0078    f J  1.8    T+4 


m 


92V2*5    expf7.77(l-    519  \\ 

;  V  V      1.8T+4927; 
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where 

T  is  temperature  in  degrees  Celcius, 
h  is  the  percent  mole  ratio  of  water  vapor, 
h=  100  (e/p) , 

e=  water  vapor  pressure  in  mb , 
P=  atmospheric  pressure  in  mb, 
To  convert  from  dB/m  to  Nepers/m  note  that 
10  log(I(x)/I(0)  )=  10  log(exp(-c<x)  , 
==>  4.34  X(Nepers/m)=  X  (dB/m) . 

The  Napier  frequency  is  the  frequency  of  maximum 

absorption  per  wavelength  and  ^max  is  the  attenuation  at 
the  Napier  frequency.   The  Napier  frequency  is  shifted  to 
higher  frequencies  by  even  small  amounts  of  water  vapor. 
Plots  of  attenuation  versus  water-vapor  pressure  in 
millibars  (Figure  1)  and  versus  relative  humidity  (Figure  2) 
are  shown  for  a  range  of  frequencies. 

For  a  given  temperature,  the  attenuation  has  a  maximum 
and  decreases  for  higher  or  lower  humidities.   With  higher 
temperatures,  the  maximum  attenuation  increases  and  the 
relative  humidity  at  which  that  maximum  occurs  decreases. 
Plots  of  attenuation  versus  water-vapor  pressure  in 
millibars  (Figure  3)  and  versus  relative  humidity  (Figure  4) 
are  shown  for  a  range  of  temperatures. 
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Fig.  1.  Attenuation  versus  Water-vapor  Pressure  (mb.) 
for  Frequencies  around  1.6  KHz. 
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Fig.  2.  Attenuation  versus 
Relative  Humidity  (%) 
for  Frequencies  around  1.6  KHz. 
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3.    Attenuation   versus  Water-vapor  Pressure    (mb.) 
for  Temperatures   Around    20   Celsius. 
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C.   EXCESS  ATTENUATION 

The  atmosphere  scatters  and  attenuates  the  acoustic 
energy.   In  the  previous  section  on  scattering,  we  were  only 
concerned  with  the  energy  scattered  from  a  given  volume. 
Excess  attenuation  accounts  for  scattering  and  energy  loss 
for  the  round  trip  in  the  atmosphere  to  and  from  a  given 
volume  of  interest.   Excess  attenuation  arises  because  the 
atmosphere  degrades  the  mutual  coherence  of  the  acoustic 
wave.   The  divergent  solid  angle  of  the  acoustic  wave  is 
larger  than  would  occur  for  a  coherent,  diffraction  limited 
wave. 

The  excess  attenuation,  Ze  ,  accounts  for  this  energy 
lost  due  to  small-angle  scattering.   This  turbulent  beam 
broadening  reduces  the  on  axis  intensity.   Clifford  and 
Brown  in  Reference  8  [p.  1972]  develop  the  equation 

Ze  =  1/Q+N)    for  N<1, 

Ze  =1.5/(1+N)   for  N>1, 
where 

N  =  (Do/^e>2  ' 

Dn=  antenna  diameter, 

y^Qg  is  the  atmospheric  acoustic  coherence  length 

=  1.46*k2  *0j^R0  ds*Cne2(s)*(^lzS/)5/3  +  /:  sj5/3)~3/5  . 

2R0       2R0 
The  term  varies  from  a  value  of  one    for  no  excess 
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attenuation  to  an  asymptote  of  zero,  implying  the  energy 
would  be  spread  over  a  2  TTsolid  angle. 

There  is  a  step  in  the  functional  dependence  of  the 
excess  attenuation  when  N=l  that  can  be  seen  in  the 
equations  above  and  in  the  plots.   There  clearly  is  no 
physical  discontinuity  but  rather  a  transition  between  the 
theoretical  dependence  between  two  asymptotic  regions. 

Figures  5,6,  and  7  are  plots  of  excess  attenuation 
versus  range  for  antenna  sizes  of  .5,  1,  and  1.5  meters, 
respectively,  for  several  frequencies.   As  can  be  seen  from 
these  plots  or  the  equations  above,  the  relationship  of  the 
antenna  size  to  the  coherence  length  has  a  significant 
effect  on  the  excess  attenuation.   A  larger  antenna  will 
not,  by  itself,  produce  greater  range.   The  antenna  size 
must  be  matched  to  the  transmission  frequency  in  terms  of 
antenna  design  guidelines  and  effects  upon  the  excess 
attenuation. 
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Ill  EXPERIMENT 

The  acoustic  sounder  is  an  excellent  tool  for  probing 
the  lower  troposphere.   It  suffers  from  several 
shortcomings: 

1)  The  range  has  typically  been  limited  to  a  few  hundred 
meters . 

2)  It  has  proven  difficult  to  quantify  the  measurements 
accurately  for  a  given  range. 

In  this  thesis,  we  approached  the  problem  from  several 
directions.   We  looked  at  the  hardware  to  understand  the 
signal  to  noise  limitations  and  to  understand  how  the 
magnitude  of  the  returned  signal  could  be  calibrated.   We 
developed  a  software  model  which  would  allow  one  to  estimate 
the  power  backscattered  from  the  atmosphere  at  a  given 
range,  based  on  profiles  of  atmospheric  characteristics  and 
input  parameters  for  the  acoustic  radar. 

A .   HARDWARE 

The  acoustic  sounder  we  were  working  with  was  the 
Aerovironment  Model  300.   It  consists  of  an  electronic 
module  which  generates  a  1600  Hertz  electrical  pulse.   The 
pulse  is  converted  to  acoustic  energy  by  a  transducer  which 
feeds  a  1.25  meter  parabolic  reflector.   Energy 
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backscattered  by  the  atmosphere  is  then  received  by  the 
reflector  and  transducer.   The  electrical  signal  is  then 
filtered  and  amplified.   In  addition,  a  ramp  amplifier 
compensates  for  the  1/r  decreasing  signal  amplitude  with 
range  to  decrease  the  dynamic  range  requirements  of  the 
present  data  display,  a  strip  chart  recorder. 

We  replaced  the  various  integrated  circuit  amplifiers 
and  filters  with  more  current  designs  with  lower  noise.   We 
replaced  the  pre-amplif ier  with  an  OPA  111  and  the  rest  with 
LF  356  BN  devices. 

We  traced  the  amplifier  and  filters  of  the  receiver 
board  to  understand  the  undocumented  choices  the 
manufacturer  had  made.   Figures  8  and  9  show  a  pre- 
amplifier, high  and  low  pass  filters  and  two  stage 
amplification.   The  filters  are  bi-quad  configured  with 
notched  outputs.   Figure  10  shows  the  notching.   Figure  10 
represents  the  frequency  spectrum  output  of  the  receiver 
board  as  measured  with  the  HP  3561A  Frequency  Spectrum 
Analyser  with  random  noise  from  the  HP  3561A  providing  the 
input  signal  before  the  filters.   Figure  11  represents 
frequency  spectrum  of  the  receiver  board  with  the  random 

noise  across  the  transducer.   A  105  ohm  resistor  was  used  to 
match  impedences  as  shown  in  Figure  12. 


23 


<Hh 


Fig. .8, , Schematic  of  the  Receiver  Board 
of  the  Aerovironment  Model  300. 
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Fig.  9.  Schematic  of  the  High  and  Low  Pass  Filters  of 
the  Receiver  Board  of  Aerovironment  Model  300. 
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Fig.  10.  Frequency  Spectrum  of  the  Receiver  Board  of 
the  Aerovironment  Model  300.  (Input  before  filters) 
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Fig.  11.  Frequency  Spectrum  of  the  Receiver  Board  of 
the  Aerovironment  Model  300.  (Input  across  transducer) 
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Fig.  12.  Connection  of  Noise  Source  to  Transducer 
The  manufacturer  had  a  filter-oscillator  board  after  the 
receiver  board  with  a  switch  for  three  pass  bands;  narrow, 
medium,  and  wide.   With  the  same  input  as  in  Figure  10  we 
measured  the  frequency  output  for  the  three  settings  of  the 
filter-oscillator  board  (Figure  13).   As  can  be  seen,  there 
is  little  difference  between  the  three  settings. 

From  previous  measurements,  the  resistor  labeled  R--  in 
the  schematic  (Figure  8)  of  the  receiver  board  was  found  to 
load  down  the  input.   The  resistor  was  initially  25  kilo- 
ohms.   Figures  14,  15,  and  16  represent  the  signal  from  the 
input  transformer  with  R?c.    =25K,  100K,  and  ^S  ohms 
respectively.   Random  noise  from  the  HP  3561A  was  input 
across  the  transducer  as  shown  in  Figure  12.   The  Q  of  the 
input  transformer  was  improved  by  increasing  the  resistance. 
The  lower  curve  in  each  Figure  is  with  no  noise  signal  input 
from  the  HP  3561A. 

It  appeared  the  Q  of  the  filters  could  also  be  increased 
by  adding  the  resistors  labeled  RQ  in  Figure  17.   The  upper 
curves  in  Figures  18,  19,  and  20  are  the  output  of  the 
receiver  board  with  each  Rq  =  1**,    750,  and  560  ohms 
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Pig.  13.  Frequency  Spectrum  of  the  Filter  Oscillator  Board 

of  the  Aerovironment  Model  300. 

(Input  before  the  filters  of  the  receiver  board) 
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Fig.    14.    Frequency   Spectrum  of    the   Input  Transformer 
(Input   across    transducer  with   R2c_25t_      \ 
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Fig.  15.  Frequency  Spectrum  of  the  Input  Transformer 
(input  across  transducer  with  1*25=100*   ) 
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Fig.  16.  Frequency  Spectrum  of  the  Input  Transformer 
(input  across  transducer  with  1*25=  ^) 
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Plg.._17;  Schematic- of  the  Receiver  Board  with  RQ 
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Fig.    18.    Frequency   Spectrum  of    the   Receiver  Board   with   RQ=  •>£> 

(input  before    filters) 
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Fig.    19.    Frequency   Spectrum  of    the   Receiver   Board   with 
r0=750  -A  (input  before    filters) 
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Fig.    20.    Frequency   Spectrum  of    the   Receiver  Board   with 
rq=560  JX  (input  before   filters) 
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respectively.   The  input  was  random  noise  from  the  HP  3561A 
after  the  pre-amplif ier  but  before  the  filters  of  the 
receiver  board.   The  effect  of  decreasing  the  resistance 
was  an  increase  in  gain  but  the  pass-band  was  no  longer  flat 
on  top.   However,  for  the  range  of  vertical  wind  velocities, 
about  two  or  three  meters  per  second,  a  pass-band  of  fifty 
or  sixty  Hertz  is  adequate  to  allow  for  the  Doppler  shift. 
As  can  be  seen,  the  filters  could  be  tuned  for  a  narrow 
pass-band  which  would  have  better  shape.   We  would  gain 
amplification  and  still  have  a  pass-band  that  is  wide  enough 
to  meet  our  needs. 

The  lower  curves  in  Figures  18,  19,  and  20  are  with 
notching,  the  upper  without.  The  effect  of  the  notching  was 
a  sharper  cut  off  near  the  bandpass  frequency  (1600  Hz)  but 
a  loss  at  higher  and  lower  frequencies. 

B .   SOFTWARE 

The  signal  received  by  an  acoustic  sounder  from  a  given 
range  tells  us  that  scattering  centers  exist  but  we  can  say 
little  about  the  size  and  magnitude  of  the  scattering 
centers.   We  don't  know  how  much  energy  was  incident  on  the 
scattering  volume  nor  how  much  the  signal  was  degraded  on 
the  return  path.   If,  on  the  other  hand,  previous  signals 
are  used  to  determine  the  atmospheric  characteristics  then 
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an  estimate  of  the  degradation  of  the  incident  power  and  the 
degradation  of  the  return  signal  can  be  made. 

The  software  programs  I  have  written  here  go  in  the 
opposite  direction.   Given  certain  atmospheric 
characteristics,  the  power  returned  is  estimated.   With  the 
modules  tested,  it  would  then  be  a  matter  of  turning  it 
around  to  take  actual  data  and  estimate  atmospheric 
conditions.   In  the  present  form  they  serve  our  purposes  in 
allowing  us  to  explore  the  effects  of  parameter  changes. 

All  the  programs  are  fundamentally  built  around  the 

echosonde  equation,  also  sometimes  referred  to  as  the  radar 

equation  in  meteorology  [Ref.  9  and  Ref.  7:p.  3]. 

Pr  =  Er*(Pt*Et)*(exp(-2^R)  )*  (J"0(R,f  )*(q^T  )*(AjL^)*Ze2, 

2      R 

where 

Pr  is  the  power  returned  from  a  range  R, 

Pt  is  the  power  transmitted  at  frequency  f, 

Er  is  the  efficiency  of  conversion  of  acoustic  power 

to  electrical  power  by  the  transducer  , 

Et  is  the  efficiency  of  conversion  of  electrical  power 

to  acoustic  by  the  transducer, 

c   is  the  speed  of  sound  in  m/sec, 

T  is  the  pulse  length, 

A   is  the  area  of  the  antenna, 

R   is  the  range,  and 

G   is  the  effective-aperture  factor  of  the  antenna. 
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Transducer  efficiencies  must  be  measured  for  each  driver 
and  typically  range  from  a  few  percent  [Ref.  7:p.  3]  to 
twenty  five  percent  [Ref.  10:p.  11-10] . 

The  scattering  cross  section  per  unit  volume  (0  q'  i-s 

the  fraction  of  incident  power  backscattered  per  unit 

distance  into  a  unit  solid  angle  at  a  given  frequency.   From 

Reference  7  [p.  4]  and  Reference   2, 

if  Q    =  .0039  *  k1/3  *  Ct2  , 

i0 
where 

k  is  the  wavenumber  =  2  ^/wavelength, 

Tn  is  the  local  mean  temperature  in  Kelvin, 

Ct2  is  the  temperature  structure  parameter, 

The  power  scattered  from  a  scattering  volume  is 

Pb(I)  =  (Pt*Et-Pb(I-l)  )*exp(-2'XR)*(cc0/2)*Ze2*  X  , 

The  power  returned  to  the  antenna  is 

Pr=Pb*A*G*Er/R2, 
where  the  return  path  attenuation  was  already  included  in 
Ps .   The  excess  attenuation  Ze2  was  discussed  in  the 
background  section. 

The  dependence  with  height  of  the  velocity  structure 
parameter  Cv2  was  needed  for  the  calculation  of  the  excess 
attenuation.   Reference  2  and  Reference  11  [p. 149]  give 

Cv2  =  2  *^  2/3, 
where  t£  is  the  dissipation  rate  of  turbulent  kinetic  energy. 


3 


"> 


Reference  11  [p.  154]  and  Reference    (p  194)  give 

<£  =  (.33m/s)3   *  (1  +  .07  *  R3/5)3/2, 
k  *  R 

for  a  stable  surface  layer.   Figure  21  is  a  plot  of  range 

versus  this  Cv  . 


The  programs  have  four  temperature  structure  parameter 
(Cm2)  profiles.   The  operator  must  choose  one.   The  first  is 
based  on  data  as  presented  in  Reference  13  [p.  398]  for 
midday  clear  weather  above  the  Tularosa  Basin  desert  in  New 
Mexico.   The  second  choice  is  for  the  same  data  multiplied 
by  a  factor  two  to  approximate  looking  up  a  convective 
plume.   The  third  case  for  a  nocturnal  atmosphere  assumes  a 
dependence  proportional  to  the  negative  exponential  of  the 
range  as  presented  in  Reference  13  [p.  399]  with  tower  data 
from  the  same  reference  used  for  the  first  sixty-five 
meters.   The  fourth  case  assumes  a  dependence  proportional 
to  height  to  the  -4/3  and  a  surface  vertical  heat  flux  of 
.095  [Ref.  7:p.  7],   Also  case  four  allows  for  the  operator 
to  input  the  height  of  an  inversion  layer  with  Cm   being 
proportional  to  height  to  the  -4/3  above  the  inversion 
layer.   Figures  22,  23,  24,  and  25  are  plots  of  these  four 
profiles  . 

„  2        o 

Cv   and  Ctz  were  then  used  to  calculate  the  acoustic 

refractive  index  factor  [Ref.  14:p.  119] 
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Fig.    21.    Velocity   Structure   Parameter    (C    2)    Profile. 
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Fig.  22.  Range  versus  Temperature  Structure  Parameter  (CT2) 

for  Profile    1 
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Fig.  23.  Range  versus  Temperature  Structure  Parameter  (C  2) 

for  Profile  2 
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Fig.  24.  Range  versus  Temperature  Structure  Parameter  (C_,2) 

for  Profile  3 
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Fig.  25.  Range  versus  Temperature  Structure  Parameter  (CT2) 

for  Profile  4 
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Cnez  =  Ctz   +  Cv| 
4T02    c0 

Figure  26  is  a  plot  of  Cne^  using  the  third  profile  for  Ct  . 

Figure  27  represents  an  input  flow  chart.   It  summarizes 
the  effect  of  each  variable  or  atmospheric  parameter. 
Figures  28  and  29  are  plots  of  the  range  versus  the  power 
returned  to  the  echosounder  for  profiles  3  and  4. 

The  segments  of  the  programs  were  tested  against 
existing  data  to  verify  proper  operation.   Based  on  input 
temperature,  atmospheric  pressure,  water-vapor  pressure,  and 
the  frequency  of  operation  of  the  acoustic  sounder  the 
programs  calculate  the  attenuation.   If  the  operator 
desires,  the  first  program  will  plot  the  attenuation  as  a 
function  of  absolute  water-vapor  pressure  and/or  relative 
humidity  for  frequencies  at  one-third  octaves  around  the 
input  frequency  and  then  again  for  temperatures  at  ten 
degree  Celcius  intervals  around  the  input  temperature. 
These  plots  were  used  to  check  the  attenuation  against  data 
[Ref .  15]  and  [Ref .  5] . 

Ambient  noise  levels  of  acoustic  sounders  are  found  to 
be  about  ten  to  forty  dB  above  the  theoretical  Johnson  noise 
limit  [Ref.  16:p.  19-4].   This  noise  level  determines  the 
maximum  range.   This  maximum  range  is  comparable  to  range 
capabilities  of  the  Aerovironment  System  300  when  the 
operating  parameters  of  the  Aerovironment  are  used  in  the 
program. 
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Fig.    26.    Acoustic   Index   Structure   Parameter    (C      2)    Profile 

using   Ct     Proflle    3 
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Fig.  27.  Input  Flow  Chart  for  Computer  Model 
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Fig.  28.  Range  versus  Power  Returned  to  the  Antenna 
Profile  3,   Parameters  of  the  Aerovironment 
Model  300,  and  S.T.P. 
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Fig.  29.  Range  versus  Power  Returned  to  the  Antenna 

for  CT2  Profile  4,  Parameters  of  the  Aerovironment 

Model  300,  and  S.T.P. 
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Trends,  as  input  parameters  where  changed,  were  also 
used  to  check  the  program  output.   As  the  frequency 
increased,  the  backscatter  of  the  turbulence  in  the 
atmosphere  decreased  slowly  and  the  absorption  of  the 
atmosphere  increased  rapidly  in  agreement  with  Reference  17 
[p.  1-10]. 

The  calculation  for  excess  attenuation  was  difficult  to 
check.   At  ranges  of  about  400  to  500  meters  the  excess 
attenuation  of  a  typical  acoustic  radar  has  been  found  to  be 
about  .25,  in  agreement  with  the  program  calculation. 
PROGRAM  1 

The  first  program  takes  atmospheric  parameters  and  the 
echosounder  parameters  as  input  and  outputs  nine  plots.   The 
inputs,  by  the  operator  at  the  keyboard,  are: 
Atmospheric  parameters; 

Atmospheric  pressure  in  millibars, 

The  profile  of  Ct   from  four  options, 

Temperature  in  degrees  Celsius, 

Water-vapor  pressure  in  millibars, 
Echosounder  data; 

Antenna  diameter  In  meters, 

Frequency  of  the  echosounder  in  Hertz, 

Power  transmitted  by  the  echosounder  in  Watts,  and 

Pulse  length  of  the  transmitted  acoustic  energy. 
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The  program  outputs  the  following  plots: 

1.  Attenuation  (1/m)  versus  water-vapor  pressure  (mb)  for 
five  frequencies  at  one-third  octaves  around  the  input 
frequency. 

2.  Attenuation  (1/m)  versus  relative  humidity  (%)  for  five 
frequencies  at  one-third  octaves  around  the  input 
frequency. 

3.  Attenuation  (1/m)  versus  water-vapor  pressure  (mb)  for 
five  temperatures  at  ten  degree  intervals  around  the 
input  temperature. 

4.  Attenuation  (1/m)  versus  relative  humidity  (%)  for  live 
temperatures  at  ten  degree  intervals  around  the  input 
temperature . 

5.  Range  (m)  versus  excess  attenuation. 

9 

6.  Range  (m)  versus  the  temperature  structure  factor  fit    . 

7.  Range  (m)  versus  the  velocity  structure  factor  CV- . 

8.  Range  (m)  versus  the  acoustic  refractive  index 

2 
structure  factor  Cne  . 

9.  Range  (m)  versus  the  power  backscattered  to  the 
echosounder . 

All  of  the  programs  prompt  the  operator  for  inputs  and 
with  a  series  of  yes/no  questions  allows  the  operator  to  re- 
run with  the  same  inputs  or  change  the  inputs.   Figure  30  is 
a  flow  chart  of  program  one  and  appendix  1  is  a  copy. 
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Fig.  30.  Flowchart  of  Computer  Model  One 


46 


Figures  1,  2,  3,  4,  21,  22,  23,  24,  25,  26,  28,  and  29  are 
output  plots  of  program  one. 

The  second  program  takes  the  same  atmospheric  parameters 
used  in  the  first  program.   With  the  exception  of  the 
frequency,  it  takes  the  same  echosounder  parameters  also. 
This  program  outputs  a  plot  of  range  versus  frequency  and  a 
plot  of  range  versus  excess  attenuation  for  various 
frequencies.   Figure  31  is  a  flow  chart  of  the  program  and 
appendix  2  is  a  copy.   Some  output  plots  of  the  program  are 
included  in  the  conclusions  section. 

The  third  program  takes  the  same  atmospheric  and 
echosounder  parameters  as  the  first.   The  program  outputs  a 
plot  of  range  as  a  function  of  efficiency  of  the  transducer, 
assuming  the  same  efficiency  for  transmit  and  receive.  All 
the  other  programs  and  plots  in  this  thesis  assume 
efficiencies  of  25%,  which  is  on  the  high  side  of  typical 
performance.   Appendix  3  is  a  copy  and  an  output  plot  is 
presented  in  the  conclusions  section. 

The  fourth  program  has  the  same  inputs  as  the  first  with 
the  exception  of  antenna  area,  which  is  the  dependent 
variable  for  the  output  plot.   It  plots  the  range  as  a 
function  of  antenna  area  for  several  frequencies  and  is 
presented  and  discussed  in  the  conclusions.   The  model  does 
not  include  different  efficiencies  based  on  optimum  antenna 
design  for  different  frequencies.   The  output  reflects  the 
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Fig.  31.  Flowchart  of  Computer  Model  Two 
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effect  of  different  antenna  diameters  on  the  excels 
attenuation  and  the  echosonde  equation.   I  am  uncertain  ho// 
much  these  drive  the  considerations  for  optimum  antenna 
design.   Appendix  4  is  a  copy  of  the  program. 

The  fifth  program  has  the  power  to  the  transducer  as  the 
dependent  variable  plotted  against  range  for  several 
frequencies.   Appendix  5  is  a  copy  and  the  output  plot  is 
presented  and  discussed  in  the  conclusions  section . 

The  sixth  program  has  the  ambient,  atmospheric 
background  noise  as  the  dependent  variable.   Appendix  6 
contains  a  copy  of  the  program  and  the  output  plot  is 
presented  and  discussed  in  the  conclusions  section. 

These  programs,  taken  together,  allow  a  parametric 
analysis  of  the  effect  of  different  parameter  changes. 
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IV  CONCLUSIONS 

Reference  10  [p.  II-l]  indicates  that  as  the  frequency 
increases: 

1)  Background  and  wind  noises  decrease  except  for  marked 
peaks  due  to  fans,  etc.   This  relationship  between 
noise  and  frequency  is  not  included  in  the  this  model. 

2)  The  reflectivity  of  the  turbulence  in  the  atmosphere 
decreases  slowly  and  the  absorption  of  the  atmosphere 
increases  rapidly.   The  model  agrees  with  this  for  low 
values  of  the  water-vapor  pressure.   The  effects  of 
water-vapor  pressure  on  the  attenuation  are  shown  in 
Figures  1  to  4 . 

3)  Wildlife  sounds  tend  to  increase  with  increasing 
frequency  and  dominate  the  background  noise  at  about 
3000  Hertz.   This  effect  is  not  included  in  the  model. 

4)  The  Doppler  shift  requires  the  receiver  bandwidth  to  be 
increased.   The  model  doesn't  take  this  into 
account.   The  bandwidth  is  used  to  calculate  the 
Johnson  noise  but  Johnson  noise  is  not  significant. 

The  noise  really  should  be  scaled  with  frequency  to 
depict  the  gain  in  range  accurately  for  lower 
frequencies . 
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Transducer  efficiencies  vary  with  frequency.   This  is 
not  included  in  the  model.   An  increase  in  efficiency  will 
have  a  much  greater  impact  on  the  potential  range  than 
increases  in,  say,  power.   An  increase  in  efficiency  will 
help  both  with  transmission  and  return,  increasing  the 
transmitted  power  and  the  returned  electrical  signal 
strength.   Figures  32  is  a  plot  of  the  maximum  range  for 
efficiencies  of  .05  to  .5.   This  figure  is  the  output  of 
program  3.   The  range  increases  quickly  with  improvements  in 
efficiency.   The  efficiency  could  be  improved  by  using 
better  designed  horns,  such  as  catenoidal  or  exponential. 

Optimum  antenna  diameters  vary  with  frequency.   The 
output  of  program  4,  Figure  33,  demonstrates  this  effect. 
Not  included  are  the  effects  on  the  antenna  effective 
aperture  factor  G,  except  to  the  extent  it  may  be  effected 
by  the  excess  attenuation.   The  discontinuities  in  the 
curves  are  due  to  the  discontinuities  in  the  equations  for 
the  excess  attenuation.   For  a  given  frequency,  the  range 
increases  arid  then,  as  antenna  diameters  increase  beyond  an 
optimum,  the  range  decreases.   The  decrease  for  larger 
antenna  diameters  is  the  result  of  the  excess  attenuation. 

One  other  consideration  as  to  choice  of  frequency  is 
that  the  resolution  is  increased  with  increasing  frequency. 
A  specific  need  for  the  echosounder  may  drive  this 
constraint  and  is  not  considered  in  this  model. 
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Figure  34  is  a  plot  of  range  as  a  function  of  frequency 
for  an  antenna  size  of  1.5  meters.   The  water-vapor  pressure 
input  was  ten  millibars.   Figure  35  is  for  the  same  antenna 
size  but  with  a  water-vapor  pressure  of  2  millibars.   The 
step  in  each  curve  is  not  a  physical  effect  but  rather  due 
to  the  step  of  1.5  in  the  theoretical  equation  for  the 
excess  attenuation.   Figure  1  showed  the  greater  value  of 
the  attenuation  at  around  2  mb.  water-vapor  pressure  for  the 
range  of  frequencies  we  are  dealing  with.   These  previous 
two  plots  demonstrate  the  effect.   Figures  34  and  35  are 
some  of  the  output  plots  of  program  2.   Both  figures 
demonstrate  the  greater  attenuation  of  the  acoustic  wave  as 
the  frequency  increases. 

Significant  increases  in  range  can  be  achieved  by 
judicious  choice  of  frequency.  Cons  iterations  are  the 
water-vapor  pressure  (Figures  1  to  4),  the  excess 
attenuation  (Figures  5  to  7),  which  depends  on  the  antenna 
diameter  (Figure  33),  the  frequency  (Figure  33  to  35),  and 
the  anticipated  frequency  spectrum  of  the  background  noise. 

Figure  36  models  the  effect  of  increasing  the 
transmitted  power  to  increase  the  range.   This  is  the  output 
plot  of  program  5.   The  slope  of  the  curve  is  not  very  steep 
and  becomes  less  so  for  higher  frequencies.   A  considerable 
increase  in  power  is  required  to  double  the  range.   Also  the 
model  does  not  include  nonlinear  transducer  effects  which 
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arise  when  the  power  is  increased.   Reference  1  [p.  57] 
points  oat  that  at  some  point  the  acoustic  wave  becomes 
distorted,  which  implies,  from  Fourier  analysis,  there  is  a 
flow  of  energy  out  of  the  fundamental  frequency  into  the 
higher  harmonics.   Since  the  higher  frequencies  are 
attenuated  more  quickly  we  will  soon  reach  a  point  of 
saturation,  which  Brown  [Ref.  1  :p  57]  called  nonlinear 
saturation  to  distinguish  it  from  another  saturation  effect 
he  discusses. 

Figure  37  is  the  output  of  program  6.   This  shows  the 
effect  of  decreasing  the  ambient  background  noise  level. 
This  decrease  in  noise  might  be  achieved  by  using  digital 
processing  and  fast  Fourier  transforms  to  achieve  a  narrower 
bandwidth.   The  bandwidth,  and  therefore  the  noise,  can  be 
much  smaller. 

These  plots  -demonstrate  which  parameter  changes  might 
best  improve  the  range  of  an  echosounder,  which  was  one  of 
our  goals  in  this  thesis.   The  modeled  performance  changes 
should  allow  for  intelligent  decisions  of  the  necessary 
parameters  for  the  expected  uses  of  the  echosoundei: . 

We  also  sought  to  explore  the  ability  to  use  the 
returned  signal  to  quantify  atmospheric  parameters 
accurately  at  a  given  range.   As  can  be  seen  from  Figure  27 
(the  input  flow  chart),  for  a  single  acoustic  radar  return 
you  could  calculate  CT   based  on  an  assumed  profile  of  Cv-. 
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However  this  need  to  assume  a  Cv2  profile  could  be 
eliminated  with  the  use  of  two  or  more  echosounders . 
Reference  11  and  many  others  describe  techniques. 

Another  approach  would  be  to  measure  the  Doppler  width. 
The  Doppler  width  is  the  spread  of  frequency  around  the 
Doppler  shifted  frequency.   Epsilon  and  therefore  Cv2  can   be 
related  to  the  Doppler  width.   With  the  use  of  fast  Fourier 
transforms  and  digital  processing  Cv   and  CT   could  be 
measured  simultaneously  with  one  echosounder. 

Using  the  returned  signal  from  lesser  ranges  the  energy 
Incident  on  a  given  volume  could  be  estimated  and  the 
degradation  of  the  return  signal  could  be  estimated.   In 
this  way  the  computer  program  would  allow  one  to  essentially 
boot-strap  up  to  a  given  range  and  more  accurately  depict 
the  atmospheric  parameters  based  on  the  returned  signal. 
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APPENDIX    A 

COMPUTER   PROGRAM    1 
FULLER,    ROBERT  PROG    1 

10  SEP    85 


**##*##PURPOSE******* 

THIS  PROGRAM  UILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
ESTIMATE  THE  RANGE.   THE  FOLLOWING  INPUTS  ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1 )ATOMOSPHERIC  PRESSURE  in  millibars 

Z)SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 
a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 
3) TEMPERATURE  IN  DEGREES  CELSIUS 
4)UATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNDER  DATA 

5)ANTENNA  DIAMETER  IN  METERS 

G)FREQUENCY  OF  ECHOSOUNDER  IN  Hz 

7)P0WER  TRANSMITTED  BY  ECHOSOUNDER  IN  UATTS 

8JPULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOWING  PLOTS: 

1  )ATTENUATION(  1/m.  )  VERSUS  WATER-VAPOR  PRESSURE ( mb)  FOR  FIVE 

FREQUENCIES  AT  ONE-THIRD  OCTIVES  AROUND  THE  INPUT  FREQUENCY. 
2)ATTENUATION( 1/m. )  VERSUS  RELATIVE  HUMIDITY(Z)  FOR  FIVE  FREQUENCIES 

AT  ONE-THIRD  OCTIVES  AROUND  THE  INPUT  FREQUENCY. 
3)ATTENUATION( 1/n.  )  VERSUS  WATER-VAPOR  PRESSURE( mb)  FOR  FIVE 

TEMPERATURES  AT  TEN  DEGREE  INTERVALS  AROUND  THE  INPUT  TEMPERATURE. 
4>ATTENUATION( 1/m. )  VERSUS  RELATIVE  HUMIDITY(Z)  FOR  FIVE  TEMPERATURES 

AT  TEN  DEGREE  INTERVALS  AROUND  THE  INPUT  TEMPERATURE. 
5)RANGE(m.)  VERSUS  EXCESS  ATTENUATION 
G)RANGE<m.)  VERSUS  TEMPERATURE  STRUCTURE  FACTOR. 
7)RANGE(n.)  VERSUS  VELOCITY  STRUCTURE  FACTOR. 
8)RANGE(m.)  VERSUS  ACOUSTIC  STRUCTURE  FACTOR. 
9>RANGE(m. )  VERSUS  BACKSCATTERED  TO  ECHOSOUNDER. 


*****  VARIABLES***** 


Again      TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 
Ano_change  TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 

ANOUTHER  RUN. 
Ant  area   ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 
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Ant_diam 

Aton_pres 

Atten 

Att_freq 


! 

!   Att_max 

I 


A"t"t_terip 


Bn 


C1 

& 

CZ 

C3 

CneZ( 

*) 

CteZ<*> 
CveZ(») 


Dx 
Dy 

Epsilon 

Et 

Er 

Es 
Exc_att<*) 

F 

Fmax 

Freq 
Freq_con 


ANTENNA  DIAMETER. 
INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 
INPUT  OF  ATMOSPHERIC  PRESSURE  IN  nib. 
ATTENUATION  OF  ACOUSTIC  UAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
VARIABLE  USED  TO  DETERMINE  IF  OPERATOR 

WANTS  TO  PLOT  ATTENUATION  VERSUS  UATER_VAPOR 

PRESSURE  FOR  VARIOUS  FREQUENCIES  AROUND  THE  INPUT 

FREQUENCY.   IF  SO  THEN  PLOT  IS  DONE  IN  SUBPROGRAM 

Att_freq. 
VARIABLE  IN  SUBPROGRAM  ATTENUATION.   IT  IS  THE 

ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

ATTENUATION  FOR  THE  INPUT  CONDITIONS. 
VARIABLE  USED  TO  DETERMINE  IF  OPERATOR  UANTS 

TO  PLOT  ATTENUATION  VERSUS  UATER_VAPOR  PRESSURE 

FOR  VARIOUS  TEMPERATURE  AROUND  THE  INPUT  TEMPERATURE. 

THE  PLOT  IS  DONE  IN  THE  SUBPROGRAM  Att.temp. 
BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 

VERTICAL  VELOCITIES. 
VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 
VARIABLES  USED  TO  SCALE  INVERSION  HEIGHT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 
CONSTANT  USED  IN  CALCULATION  FOR  CteZ  PROFILE  4. 
ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 

PARAMETER.   CALCULATED  BASED  ON  SELECTION  OF 

PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 

AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 

PARAMETER. 
ARRAY  OF  VALUES  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 
ARRAY  OF  VALUES  OF  THE  VELOCITY  STRUCTURE   PfiRAMETER. 

VALUES  BASED  ON  CALCULATION  USING  ASSUMED  DISSAPATION 

RATE. 
STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 
STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 
DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 
TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 
EFFICIENCY  OF  ECHOSOUNDER  UHEN  RECEIVING  ACOUSTIC 

BACKSCATTER. 
SATURATION  VAPOR  PRESSURE  AT  GIVEN  TEMPERATURE. 
EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 
USED  TO  INSURE  VALUE  OF  FREQUENCY  PASSED  TO  SUBPROGRAM 

Att_freq  WAS  NOT  CHANGED. 
ANTENNA  EFFECTIVE  APERATURE  FACTOR. 
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H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POUER  UNTIL 

POUER  UAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 
Interval    MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 
Inver      HEIGHT  OF  INVERSION  LAYER. 
J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  UAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Mess_up    TEST  VALUE  TO  CHECK  RESPONSES  UHEN  OPERATOR  IS  INPUTTING 

RESPONSES. 
N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 
Neu_va     TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 

A  VARIABLE  BEFORE  A  NEU  RUN. 
Noise      ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 
Pou_back<*)  POUER  BACKSCATTERED  FROM  GIVEN  RANGE. 
Pow_ret(*)  POUER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 
Pou_trans   INPUT  OF  POUER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 
Profile    OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 
Pstar      VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 
Pulse      TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.   INPUT  BY 

OPERATOR. 
R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1         VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE 
Range(*)    ARRAY  OF  RANGE  VALUES. 
Remainder   REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 

NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 

ATTENUATION. 
Rge        VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 
Rho        CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 
Sigma<*)    FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 
Speed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 
Sumpou_back  SUM  OF  THE  BACKSCATTERED  ENERGY 
T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
Temp_con   USED  TO  INSURE  VALUE  OF  TEMPERATURE  PASSED  TO  SUBPROGRAM 

Att_temp  UAS  NOT  CHANGED. 
Temp$      VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Title*     STRING  PASSED  TO  SUBPROGRAM  PT  FOR  TITLE  OF  PLOT. 
Tstar      INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
Var        USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  WISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 
Uat_pres    ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 
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RESPONSE  TO  YES  OR  NO  QUESTION. 

Xlabel$  LABEL  ON  X  AXIS  PASSED  TO  SUBPROGRAM  Pt  FOR  PLOTTING. 

Xmax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xmin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xrange  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

Xvar<*)  VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Pt 

YlabelS  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Pi  FOR  PLOTTING. 

Ymax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Ynin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

Ze  INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


"DECLARE  VARIABLES 

INTEGER  I 

DIM  Poui_back(  1  500)  ,Poui_ret<  1500)  ,Range<  1500)  ,CteZ(  1500)  ,Sigma(  1500) 

DIM   CveZ( 1500),CneZ<  1 500) ,Exc_att( 1500) ,Xvar<  1500) 

DIM   Ti-tle$[50]  ,Xlabel$[50]  ,Ylabel$t  1G1 

I 

PLOTTER  IS  705,"HPGL" 
LINE  TYPE  1 


! INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" ,Temp 

INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" ,Atom_pres 

INPUT  "ENTER  WATER  VAPOR  PRESSURE  IN  MILLIBARS" , Wat.pres 


INPUT  ECHOSOUNDER  DATA 
INPUT  "ENTER  FREQUENCY  OF  ECHOSOUNDER  IN  HERTZ", Freq 
INPUT  "ENTER  ANTENNA  DIAM  IN  METERS" , Ant_diam 

INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS=" ,Pul se 
INPUT  "ENTER  POWER  TO  TRANSMITTER  IN  WATTS" ,Pou_trans 
Et=.Z5        "TRANSMIT  EFFICIENCY 
Er=.Z5        (RECEIVER  EFFICIENCY 
G=.40        ! ANTENNA  EFFECTIVE  APERATURE  FACTOR 


SELECT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 


THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 
YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 
PROFILE  TO  BE  USED." 
YOUR  SELECTIONS  ARE" 

1       A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS 
PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398  WHICH  ' 
GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

Z  THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 
WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE  LOOKING" 
UP  A  CONVECTIVE  PLUME" 
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PRINT  "3ft  TEMPERftTURE  STRUCTURE  PROFILE  BASED  ON  DftTft  ftS" 

PRINT  "  PRESENTED  IN  UflLTERS/KUNDEL  1981  PAGE  398  FOR" 

PRINT  "  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 

PRINT  "  IS  EXP<-.001*HEIGHT  ABOVE  B5  METERS  AND  HEIGHT" 

PRINT  "  TO  THE  -1.4G  BELOU  B5  METERS" 

PRINT  "  " 

PRINT  "     4      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 

PRINT  "  HEIGHT  TO  THE  -4/3  " 

I 

PRINT  "  " 

PRINT  "  " 

Mess_up=1 

WHILE  Mess_up=1 

INPUT  "ENTER  THE  DESIRED  PROFILE  (1  OR  Z  OR  3  OR  4)", Profile 
IF  Profile=4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" , Inver 
END  IF 
IF  Profiled  THEN 

Mes6_up=0 
ELSE 

IF  Profile=2  THEN 

Mess_up=0 
ELSE 

IF   Profile=3    THEN 

Mess_up=0 
ELSE 

IF  Profile=4  THEN 

Mess_up=0 
ELSE 

PRINT  Profile,"  UAS  NOT  ONE  OF  THE  OPTIONS!!!! 
Mess_up= 1 
END  IF 
END  IF 
END  IF 
END  IF 
END  UHILE 
OUTPUT  KBD;"K"; 
! 
I 

Again=1 

(INITIALIZE  ARRAYS  FOR  SUCCESSIVE  RUNS 
UHILE  Again=1 
FOR  J=1  TO  I 

Pou_back( J>=0 
Pou_ret(  J)=0 
Range( J)=0 
CteZ(J)=0 
Sigma< J )=0 
CveZ( J)=0 
CneZ( J)=0 
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Exc_att< J)-0 
NEXT  J 
FOR  J-1  TO  15 

PRINT  "  " 
NEXT  J 

PRINT  "DO  YOU  WISH  TO  SEE  A  PLOT  OF  ATTENUATION  VERSUS  HUMIDITY" 
PRINT  "FOR  VARIOUS  FREQUENCIES  AROUND  THE  FREQUENCY  YOU  INPUT?" 
LINPUT  "IF  YES  ENTER  Y  ,  IF  NO  ENTER  N",X$ 
Att_freq=FNYes(X$> 
IF  Att_freq=1  THEN 

OUTPUT  KBD;"K"; 

CALL  Att_f  req(  Atom_pres  ,Freq,  Temp  ) 
END  IF 

OUTPUT  KBD;"K"; 
! 
! 
FOR  J-1  TO  15 

PRINT  "  " 
NEXT  J 

PRINT  "DO  YOU  WISH  TO  SEE  A  PLOT  OF  ATTENUATION  VERSUS  HUMIDITY" 
PRINT  "FOR  VARIOUS  TEMPERATURES  AROUND  THE  TEMPERATURE  YOU  INPUT?" 
LINPUT  "IF  YES  ENTER  Y  ,  IF  NO  ENTER  N",X$ 
Att_tenp=FNYes(X$) 
IF  Att_temp=1  THEN 

OUTPUT  KBD;"K"; 

CALL  Att_tenp(  Atom_pres , Freq , Temp ) 
END  IF 

OUTPUT  KBD;"K"; 
! 

(CONVERT  TEMPERATURE  TO  KELVIN 

T=Temp+Z73 

! 

(CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE( CELCIUS ) 

Speed_sound=Z0.05»<  T)".5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 

C=Z0.05*Z73A.5 

! 

(CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
(EQUATION  IS  FROM  NEFF  1975 

CALL  Attenuat  ion(  Aton_pres , At ten, Freq, Temp, Uat_pres ) 
1=1     (INDEX  FOR  LOOP 

Pou_back(0)=0     (INITIALIZE  VARIABLE  FOR  POUER  BACKSCATERRED 
K=Z*PI*Freq/Speed_sound     (UAVENUMBER 
Ant_area=PI«( Ant_dian/Z ) AZ    (ANTENNA  AREA 
Interval=( Speed_sound*Pul se* 1 .E~3)/Z 

Bn=Z*Freq*( 1  -  1 /(  G/Speed_sound+ 1 ) ) (BANDUIDTH  FOR  3M/S  VERTICAL  VELOCITIES 
Noise=1 .38E-Z3»Bn»(T)+Z.E-14!MINIMUM  DETECTABLE  SIGNAL  AS 
JOHNSON  NOISE<NEGLIGABLE)4  ESTIMATED  BACKGROUND 

CALCULATE  RETURNED  SIGNAL  POUER  UNTIL  IT  IS  LESS  THAN  BACKGROUND. 
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Range(0)=0 
Suropou_back  =0 
REPEAT 

Range< I )=Range(  I-  1  HZ 
SELECT  Profile 
CASE  1 

!THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  UEATHER  ABOVE  THE 
ITULAROSA  BASIN,  NEU  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 
! AVERAGING  TIME. 

CteZ(I)=Z.1Z*Range(I)A<-1 .IB) 
CASE  Z 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  UEATHER  ABOVE  THE 
ITULAROSA  BASIN,  NEU  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 
•AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  UAS 
'INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  CONVECTIVE  PLUME 

CteZCI )-Z*Z. 1Z*Range(I )A(-1 . IB) 
CASE  3 

!THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
(PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
! NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
!IS  EXP(-. 001 "HEIGHT  ABOVE  65  METERS  AND  HEIGHT" 
!TO  THE  -1.4B  BELOU  G5  METERS" 
IF  RangeCIXBS  THEN 

CteZ(I)=7S.S*Range(ir<-1  .46) 
ELSE 

CteZ( I >=3.66E-Z*EXP( - .001«Range(  I ) ) 
END  IF 
CASE  4 

"CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!     THAT  IT'S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POUER 
!     THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  Cm. /sec 
!THE  DEPENDENCE  UITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
!MY  LEFT  EAR 

I  EQUATION  FROM  NEFF ,  1975 
SELECT  Range( I )/Inver 
CASE  <.9 

C3=((  .0Z4>*(Tr<  .667)) 
CteZd  )=C3*<Range<  I  )  )A<-1  .33) 
CASE  .9  TO  1 
'  CteZ(I)»CteZ<I-1 > 

'  R1=Range<I) 

)  C1=CteZ(I) 

>  CASE  1  TO  1.3 

)  CteZ< I )=10A( <LGT<C3*Range< I ) >-LGT(C3*R1 ) )+LGT(C1 ) ) 
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RZ=Range< I ) 
CZ-CteZ(I ) 
CASE  ELSE 

CteZ(I)-C3*(Range< I >*<  -  1 . 33)-RZ*< -  1  .  33 ) )  +  CZ 
END  SELECT 
END  SELECT 

! CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
(AND  TETARSKI 

Epsilon«( .Z866/Range< I  ))•<  1 + .07*< Range< I) )*< .6) )"(  1 .5) 
CveZ(  I  )-Z*Epsilon',<  .G67) 

•CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 
(FORMULA  FROM  TETARSKI  19G1 

CneZ(I)-<CteZ(I)/(Z.98E+5))+(CveZ<I)/(C»C>) 

! CALCULATE  THE  FRACTION  OF  POWER  BACKSCATTERED  FROM  INTERVAL. 
(EQUATION  FROM  NEFF .  1975 
SigmaU  )-(  .0039»<K*(  1/3)  )*CteZ(  I))/(T)*Z 
! 

(CALCULATE  THE  EXCESS  ATTENUATION.   EQUATION  FROM  CLIFFORD,  1980 
(THE  EXCESS  ATTENUATION  IS  Exc_att( I > 

!THE  MODULO  STRUCTURE  IS  TO  SKIP  SOME  OF  THE  INTEGRALS  ONCE  THE 
(RESOLUTION  IS  LESS  IMPORTANT. 
IF  K50  THEN 

Remainder=0 
ELSE 

IF  K150  THEN 

Remainder=Range( I  )  MODULO  10 
ELSE 

Remainder=Range< I )  MODULO  Z0 
END  IF 
END  IF 

IF  Renainder=0  THEN 
Rho=0 
L=0 
H=0 

Rge=Range(I)     (CONSTANT  IN  INTEGRAL 
R=0 
FOR  J=0  TO  Z»I 

F=CneZ(INT( J/Z+1  >) 

F=F«<  1-R/<Rger(  1  .  67)  +  (  R/(  Rge)  )  "  (  1  .  G7  )  ) 

IF  J>0  THEN 

IF  J<Z*I  THEN 

IF  INT(J/Z)=J/Z  THEN 
L  =  L+F 
F=0 
ELSE 

H=H+F 
F=0 
END  IF 
END  IF 
END  IF 
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Rho=Rho+F 
R=R+1 
NEXT    J 

Rho=Rho+4*L+Z*H 

Rho=< ( <  Rho* . 33 )*K*K* 1 . 46 )*  < - . 6 ) > 
N=<Ant_diam/RhorZ 
IF   N<=1    THEN 

Ze=1/< 1+N) 
ELSE 

Ze=1 .5/< 1+N) 

!STEP  OF  1.5==>SEE  CLIFFORD  1980 
END  IF 

Exc_att(I)=Ze*Ze 
ELSE 

Exc_att(  I )=Exc_att< I- 1 ) 
END  IF 
! 

! CALCULATE  THE  POWER  BACKSCATTERED 

Pou_back( I )  =  ( Pou_trans*Et-Sumpou_back )*EXP(  -At ten* Range (  I ) ) 

Pou_back< I )=Pou_back< I )*Interval*Exc_att(  I )*Signa<  I ) 

Sui»ipou_back  =  Sunpou_back+Pou_back  (  I  ) 

! CALCULATE  THE  POWER  BACKSCATTERED  TO  THE  ANTENNA 

!TWO  WAY  PATH  ATTENUATION  ACCOUNTED  FOR  ABOVE. 

Pou_ret( I )=Pou_back( I )*EXP( -Atten*Range(  I ) )*Ant_area*G*Er/Range<  I )*Z 

PRINT  "RANGE  =  ",Range(  I) 

PRINT  "POWER  RETURNED=" ,Pou_ret(I) 

1  =  1  +  1 
UNTIL  Pou_ret(  1-1 ) >0  AND  Pou_ret(  I- 1 )<=Noise 
I 
! 

OUTPUT  KBD^'K"; 
l 

PRINT  "INPUT  CONDITIONS" 


PRINT 

•■       it 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

"TEMPERATURE= 
"ATMOSPHERIC  PRESSURE= 
"WATER  VAPOR  PRESSURE= 
"PULSE  LENGTH= 
"TRANSMITTED  FREQUENCY= 
"ANTENNA  DIAMETER= 
"POWER  TRANSMITTED= 
PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED 
IF  Profile=4  THEN 

PRINT  "INVERSION  HEIGHT  ",Inver 
END  IF 

SELECT  Profile 
CASE  1 

PRINT  "      1       A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 
PRINT  "  AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 


,Tepip,"CELCIUS" 
, Atom_pres , "nb" 
, Wat_pres , "mb" 
, Pulse , "ns" 
,Freq,"  Hz" 
, Ant_diam , "  n." 
,Pou_trans,"  WATTS" 
" .Profile 
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PRINT  "  UHICH  GIVES  A  HEIGHT  TO  THE  -LIE  PROFILE" 

CASE  2 

PRINT  Z      R  TEMPERARURE  STRUCTURE  PROFILE  BfiSED  ON  DATA- 
PRINT  "  AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 
PRINT  "            UHICH  GIVES  A  HEIGHT  TO  THE  -LIB  PROFILE- 
PRINT  "  BUT  WITH  A  FACTOR  OF  TUO  TO  APPROXIMATE- 
PRINT  "            LOOKING  UP  A  CONVECTIVE  PLUME- 
CASE  3 

PRINT  "  3      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
PRINT  "  AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398" 

PRINT  "  FOR  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH" 

PRINT  "  HEIGHT  IS  EXP( - . 001 ^HEIGHT  ABOVE  65  METERS" 

PRINT  M  AND  HEIGHT  TO  THE  -1.4G  BELOU  G5  METERS- 

CASE  4 

PRINT  4      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
PRINT  "  HEIGHT  TO  THE  -4/3  " 

END  SELECT 
I 

PRINT  USING  "K'V'THE  MINIMUM  DETECTABLE  SIGNAL  UAS  SET  AT"  ,  Noise ,  "UATTS' 
! 

PRINT  "  " 

PRINT  "OUTPUT  CONDITIONS" 

! 

PRINT  USING  "K";"RANGE=" .Range(I-Z) ,"  n." 
| 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 
PAUSE 

OUTPUT  KBD; "K" ; 

PRINT  "THE  FOLLOUING  GRAPHS  UILL  NOU  BE  PLOTTED  UITH  RANGE  VERSUS" 
PRINT  "EXCESS  ATTENUATION- 
PRINT  "CTEZ  (TEMPERATURE  STRUCTURE  PARAMETER)" 
PRINT  "CVEZ  (VELOCITY  STRUCTURE  PARAMETER)" 

PRINT  "CNEZ  (ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  PARAMETER)" 
PRINT  "POUER  RETURNED" 

PRINT  "UHEN  YOU  ARE  READY  FOR  THE  FIRST  GRAPH  HIT  CONTINUE" 
PAUSE 

OUTPUT  KBD;"K"; 
I 

GRAPHICS  ON 
VIEUPORT  15, 120, 10,70 
! 
! 

(PLOT  RANGE  VERSUS  EXCESS  ATTENUATION 
Xmin=0 
Xmax=1 
Ynin=0 
IF  Range(I-Z)<=1000  THEN 

Ymax=3 
ELSE 

Ynax=4 
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END  IF 

Dx  =  .1 

Dy=1 

Xrange=Xnax-Xnin 

Yr  ange= Yniax  -  Yi»i  i  n 

UINDOU  Xmin.Xroax ,Ynin, Ynax 

AXES  Dx.Dy.Xmin.Ymin, 1 ,1 

CLIP  OFF 

! 

•LABEL  PLOT 

CSIZE  4, .6 

LDIR  0 

LORG  5 

MOVE  .5,Ynax*1 . 1 

LABEL  "RANGE  VERSUS  EXCESS  ATTENUATION" 

! LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  5 

FOR  J=0  TO  Xmax  STEP  Dx 

MOVE  J,-. 1*Dy 

LABEL  J 
NEXT  J 

MOVE  .5*Xrange,-.3*Dy 
CSIZE  4..G 

LABEL  "EXCESS  ATTENUATION" 
(LABEL  VERTICAL  AXES 
LORG  8 
FOR  J=0  TO  Ymax  STEP  Dy 

CSIZE  4, .6 

MOVE  -.3*Dx,J 

LABEL  "10" 

CSIZE  Z 

MOVE  -. 1*Dx,J+.05*Dy 

LABEL  USING  "K";J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  -Dx, .5*Yrange 
CSIZE  4, .6 

LABEL  "RANGE  (METERS)" 
CLIP  ON 
FOR  J=1  TO  1-1 

PLOT  Exc_att(  J) ,LGT< Range< J)) 
NEXT  J 
PRINT  "HIT  CONTINUE  FOR  NEXT  PLOT" 


PLOT  RANGE  VERSUS  TEMPERATURE  STRUCTURE  PARAMETER 
PAUSE 
GCLEAR 


ID 


! 

IF  Range(I-Z)<-1000  THEN 

Xnin=-4 

Ymax"3 
ELSE 

Xnin=-5 

Ynax=4 
END  IF 
Xnax*© 
Ynin=0 

Title*="RANGE  vs  TEMPERATURE  STRUCTURE  PARAMETER" 
Xlabel$="TEMPERATURE  STRUCTURE  PARAMETER" 
Ylabel$=" RANGE  (METERS)" 

CALL  Log_log(Xnin,Xmax,Ypiin,Yroax,CteZ(  *  )  ,1-1  ,Range<  *  )  ,1-1  .Title*  ,  Xlabel  $ 
1$) 
l 
! 

IPLOT  RANGE  VERSUS  VELOCITY  STRUCTURE  PARAMETER 
PAUSE 
GCLEAR 
! 

Xmin=-Z 
Xmax=0 
Ymin=0 
IF  Range( 1-2 )<= 1000  THEN 

Ynax=3 
ELSE 

Ymax=4 
END  IF 

Ti"tle$=" RANGE  VERSUS  VELOCITY  STRUCTURE  PARAMETER" 
Xlabel$="VELOCITY  STRUCTURE  PARAMETER" 

CALL  Log_log(Xnin,Xmax,Ymin,Ynax,CveZ( *) ,1-1 ,Range(*) ,  I- 1  .Title* , Xlabel* 
il$) 
! 
! 

IPLOT  RANGE  VERSUS  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  PARAMETER 
PAUSE 
GCLEAR 
Xnin=-7 
Xmax=-5 
Ymin=0 
IF  Range( I-Z)<=1000  THEN 

Ymax=3 
ELSE 

Ymax=4 
END  IF 

T it le$=" RANGE  vs  ACOUSTIC  INDEX  STRUCTURE  PARAMETER" 
Xlabel*="ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  PARAMETER- 
CALL  Log_log(Xnin,Xnax,Ymin,Ynax ,CneZ<  » ) ,1-1  ,Range( * ) ,1-1  , Title* , Xlabel* 
el$) 
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PLOT  RANGE  VERSUS  POWER  RETURNED  TO  ANTENNA 
PAUSE 
GCLEAR 
IF  Range* I -Z><= 1000  THEN 

Xmin=-15 

Xroax=-5 

Yi»iax  =  3 
ELSE 

Xmin=- 17 

Xroax=-4 

Yroax=4 
END  IF 
Ymin=0 

T it le$=" RANGE  VERSUS  POWER  RETURNED  TO  ANTENNA" 
Xlabel$="POUER  RETURNED  TO  ANTENNA" 
Ylabel$="RANGE  (METERS)" 

CALL  Log_log(  Xnin.Xmax  , Ymin, Ymax  ,Pou_re-t(  »  )  ,1-1  ,Range<  *  )  ,1-  1  ,Title$,Xlab 
VlabelS) 
PAUSE 
GCLEAR 


LINPUT  "UOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again=FNYes(X$) 
Mess_up=1 
WHILE  Mess_up=1 
SELECT  Again 
CASE  1 

LINPUT  "UOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X$ 
Neu_va=FNYes(X$) 
Ano_change= 1 
WHILE  Ano_change=1 
SELECT  Neu_va 
CASE  1 
PRINT  "       VARIABLE  CURRENT  VALUE" 

PRINT  USING  "K";"1   TEMPERATURE        ", Temp , "CELSIUS" 
PRINT  USING  "K'V'Z   ATOMOSPHERIC  PRESSURE  " , Aton_pres , "mb" 
PRINT  USING  "K";"3   WATER  VAPOR  PRESSURE  " , Wat_pres , "nb" 
PRINT  "  " 

PRINT  USING  "K";"4   FREQUENCY  OF  ECHOSOUNDER  ",Freq,"  Hz" 
PRINT  USING  "K";"5   ANTENNA  DIAMETER     " , Ant_dian, "  n. " 
PRINT  USING  "K";"G  PULSE  LENGTH         ".Pulse,"  ns" 
PRINT  USING  "K";"7   POUER  TRANSMITTED  " ,Pou_trans , "  WATTS" 
PRINT  USING  "K";"8   ATMOSPHERIC  PROFILE   ".Profile 
PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 
PRINT  "CHANGE" 
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INPUT  Var 
SELECT  Var 

CASE  1 

INPUT  "TEMPERATURE*" .Temp 
CASE  Z 

INPUT  "ATMOSPHERIC  PRESSURE  IN  mb  =  " , Aton_pres 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  rib-"  ,  Uat_pres 
CASE  4 

INPUT  "FREQUENCY  IN  Hz=",Freq 
CASE  5 

INPUT  "ANTENNA  DIAMETER  IN  m.  ="  ,  An"t_diam 
CASE  G 

INPUT  "PULSE  LENGTH  IN  r*is  =  ", Pulse 
CASE  7 

INPUT  "POUER  TRANSMITTED  IN  UATTS=" ,Pou_trans 
CASE  8 

PRINT  "  1   ==>  CtAZ  PROFILE  OF  ZA(-1.I6>" 
PRINT  "        FROM  UALTERS/KUNDEL  1981" 
PRINT  "  Z   ==>  SAME  AS  ONE  BUT  UITH  FACTOR" 
PRINT  "        OF  TUO  TO  APPROXIMATE  LOOKING" 
PRINT  "        UP  A  CONVECTIVE  PLUME" 
PRINT  "  3   =  =  >  CfZ  PROFILE  OF  EXP< - .001 »Z )" 
PRINT  "        FROM  UALTERS/KUNDEL  1981." 
PRINT  "  4   ==>  CT*Z  PROFILE  OF  ZA(-4/3)" 
INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" , Prof i 1 e 
IF  Profile=4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS*"  ,  I  river 
END  IF 
CASE  ELSE 

PRINT  Var,"IS  NOT  ONE  OF  THE  OPTIONS- 
END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE( ENTER  Y  OR  N)?",X$ 
Neu_va-FNYes(X$) 
Mess_up=Z 
CASE  Z 

Mess_up=Z 
Ano_change=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  I  Ml!" 
Ano_change=  1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)7"  ,)<$ 
Neui_va=FNYes(X$) 
END  SELECT 
END  WHILE 
CASE  Z 

Mess_up=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N   !!!!!!" 
Mess_up= 1 
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LINPUT  "MAKE  ANOTHER 
Again=FNYes<><$> 
END  SELECT 
END  WHILE 
END  UHILE 

PRINT  "THAT'S  ALL,  FOLKS" 
end 
l 


RUN  (ENTER  Y  OR  N)?"  ,X3: 


ICALCULTE  THE  ATTENUATION 

sub  Attenuat ion( Aton_pres , Atten.Freq.Tenp, Uat_pres ) 
THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 

INPUT   ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 
TEMPERATURE  IN  DEGREES  CELCIUS 
UATER-VAPOR  PRESSURE  IN  MILLIBARS 

OUTPUT   ATTENUATION  IN  1 /METERS 


VARIABLES 
Atoi»i_pres 
Atten 

Att  max 


Froax 

Freq 

H 

Pstar 

Temp 
Tstar 


INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 
ATTENUATION  OF  ACOUSTIC  UAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
VARIABLE  IN  SUBPROGRAM  ATTENUATION.   IT  IS  THE 

ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

ATTENUATION  FOR  THE  INPUT  CONDITIONS. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 
INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 


IN  SUBPROGRAM  ATTENUATION. 
Ua-t_pres    ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
H=100*Wat_pres/Atom_pres 
Tstar=( 1 .8*Tenp+49Z  >/519 
Pstar=A-tom_pres/1014 

Fnax  =  (  10+6600*H+44400*H*H>*Pstar/Tstar'>  .8 
Att_i»iax  =  .0078*Fmax*TstarA(-2.5)*EXP(7.77*(  1-1 /Tstar)  ) 
F=Freq/Fmax 

Atten=(Att_i«iax/304.8)*(  (  .  1  8*F  )  "Z  +  <  Z*F*F/(  1  +F*F  )  )  AZ  )  "  .  5 
Atten=(A-tten+1.74E-10*Freq*Freq)/4.35 
SUBEND 
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I 

def  FNYes(XS) 

ITHIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 


INPUT 


t 


X$ 


FNYes 


!     OUTPUT 

! VARIABLES 

!   Temp$      VARIABLE  STRING  USED  IN  FUNCTION  YES. 

!   XS         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

!  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  Tenp$[ 13 

Tenp$t 1 , 11=TRIM$(X$) 

SELECT  Tenp$ 


CASE 


CASE 


CASE 


CASE 


END 


"Y'W 
RETURN 
"N","n' 
RETURN 

RETURN 
ELSE 
RETURN 
SELECT 


-Z 


FNEND 


sub  Att_f req( Aton_pres ,Freq,Tenp ) 

THIS  SUBPROGRAM  PLOTS  THE  ATTENUATION  OF  SOUND  VERSUS  UATER- 
VAPOR  PRESSURE  FOR  FIVE  DIFFERENT  FREQUENCIES  AT  1/3  OCTIVE 
INTERVALS  AROUND  THE  INPUT  FREQUENCY. 

INPUT    ATOMOSPHERIC  PRESSURE  IN  MILLIBARS 
FREQUENCY  IN  HERTZ 
TEMPERATURE  IN  CELCIUS 

OUTPUT   PLOT  OF  ATTENUATION  VERSUS  UATER-VAPOR  PRESSURE 


VARIABLES 
Aton_pres 
At  ten 

Att_freq 


Es 

Freq 

Freq_con 


INPUT  OF  ATMOSPHERIC  PRESSURE  IN  nb. 

ATTENUATION  OF  ACOUSTIC  UAVE.   CALCULATED  IN 
SUBPROGRAM  ATTENUATION. 

VARIABLE  USED  TO  DETERMINE  IF  OPERATOR 
UANTS  TO  PLOT  ATTENUATION  VERSUS  UIATER_VAPOR 
PRESSURE  FOR  VARIOUS  FREQUENCIES  AROUND  THE  INPUT 
FREQUENCY.   IF  SO  THEN  PLOT  IS  DONE  IN  SUBPROGRAM 
Att_-freq. 

SATURATION  VAPOR  PRESSURE  AT  GIVEN  TEMPERATURE. 

INPUT  FREQUENCY  OF  ECHOSOUNDER. 

USED  TO  INSURE  VALUE  OF  FREQUENCY  PASSED  TO  SUBPROGRAM 
Att_freq  UAS  NOT  CHANGED. 


75 


J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Temp_con   USED  TO  INSURE  VALUE  OF  TEMPERATURE  PASSED  TO  SUBPROGRAM 

Uat_pres    ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 
Ywax       LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Freq_con=Freq 
Temp_con=Temp 
Yrcax=0 

Freq=DROUND(Freq*Z"(Z/3) ,3) 
FOR  J=0  TO  4  STEP  . 1 

Uat_pres=J 

CALL  Attenuai  ion<  Atom_pres , Atten , Freq , Temp , Uat_pres ) 

IF  A"tten>Ymax  THEN 
Ymax=Atten 

END  IF 
NEXT  J 

Ynax  =PROUND( Ymax  + . 005 ,  -  2  ) 
GRAPHICS  ON 
VIEWPORT  15,120,10,70 
WINDOW  0, 15,0,Ymax 
AXES  .5,. 005 ,0,0,2,2 
CLIP  OFF 
(LABEL  PLOT 
CSIZE  4, .5 
LDIR  0 

MOVE  7.5,Ymax*1 . 1 
LORG  5 

LABEL  "ATTENUATION  VERSUS  WATER-VAPOR  PRESSURE" 
CSIZE  4, .5 

! LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  G 
FOR  J=0  TO  IS 

MOVE  J,0 

LABEL  J 
NEXT  J 

MOVE  7.5,-Ymax*. 1 
CSIZE  4, .6 

LABEL  "WATER-VAPOR  PRESSURE    mb" 
(LABEL  VERTICAL  AXES 
LORG  8 
CSIZE  4, .6 
FOR  J=0  TO  Ynax  STEP  .01 

MOVE  .25. J 
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LABEL  J 
NEXT  J 
LDIR  PI/Z 
LORG  6 

MOVE  -Z,Ynax«.5 
CSIZE  4, .6 

LABEL  "ATTENUATION    1/m" 
! PRINT  INPUT  CONDITIONS 
LDIR  0 
LORG  1 

MOVE  5.5,Ymax 
CSIZE  4, .6 

LABEL  USING  "K"  ;  "ATMOSPHERIC  PRESSURE  =  "  ,  Atom_pres  ,  "  mb" 
MOVE  5.5,Ymax*.95 
CSIZE  4, .6 

LABEL  USING  "K" ; "TEMPERATURE=" ,Temp, "  CELSIUS" 
J=0 
FOR  L=1  TO  S 

FOR  X=0.  TO  B.5+J  STEP  .05 
Uat_pres=X 

CALL  Attenuation*  Atom_pres , Atten.Freq , Temp , Uat_pres ) 
IF  INT(L/Z)=L/Z  THEN 

LINE  TYPE  Z 
ELSE 

LINE  TYPE  1 
END  IF 

PLOT  X.Atten 
NEXT  X 
LDIR  0 
LORG  Z 
CSIZE  3, .6 
MOVE  X.Atten 
LINE  TYPE  1 

LABEL  USING  "K" ; Freq, "Hz" 
Freq=DROUND(Freq*Z*(-1/3) ,3) 
J=J+1.5 
NEXT  L 
I 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

FOR  J=0  TO  IS 

PRINT  "  " 
NEXT  J 

PRINT  "UOULD  YOU  LIKE  TO  SEE  THE  SAME  PLOT  ONLY- 
PRINT  "WITH  ATTENUATION  VERSUS  RELATIVE  HUMIDITY?" 
LINPUT  "IF  YES  ENTER  Y,  ELSE  ENTER  N",X$ 
Att_freq=FNYes(X$) 
IF  Att_freq=1  THEN 
OUTPUT  KBD;"K"; 
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UINDOU  0,1 00,0, Ynax 

CLIP  ON 

AXES  B.. 005, 0,0,2,2 

CLIP  OFF 

! LABEL  PLOT 

CSIZE  4,.G 

LDIR  0 

MOVE  50,Ymax*1 . 1 

LORG  B 

LABEL  "ATTENUATION  VERSUS  RELATIVE  HUMIDITY" 

CSIZE  4, .6 

! LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  B 

FOR  J=0  TO  100  STEP  10 

MOVE  J,0 

LABEL  J 
NEXT  J 

MOVE  50,-Ynax*. 1 
CSIZE  4,.G 

LABEL  "RELATIVE  HUMIDITY    %" 
(LABEL  VERTICAL  AXES 
LORG  8 
CSIZE  4..B 
FOR  J=0  TO  Ymax  STEP  .01 

MOVE  1  ,J 

LABEL  J 
NEXT  J 
LDIR  PI/2 
LORG  G 

MOVE  -15,Ymax*.5 
CSIZE  4, .6 

LABEL  "ATTENUATION    \/m" 
! PRINT  INPUT  CONDITIONS 
LDIR  0 
LORG  1 

MOVE  40,Ymax 
CSIZE  4, .6 

LABEL  USING  "K" ; "ATMOSPHERIC  PRESSURE=" , Atom_pres , "  mb" 
MOVE  40,Yriax*.95 
CSIZE  4,.G 

LABEL  USING  "K" ; "TEMPERATURE=" ,Temp, "  CELSIUS" 
CLIP  ON 
T=Temp+273 
Es=10A(9.4-23S3/T) 
J=0 

Freq=DROUND(Freq_con*2A<-2/3) ,3) 
FOR  L=1  TO  5 

FOR  X=0.  TO  40+ J  STEP  .5 

Ua-t_pres  =  Es*Atom_pres*X/(  (  Aton_pres-Es*(  1+X/100)  )*100) 
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CALL  Attenuat  ion(  Atom_pres  ,  At  ten,  Freq,  Temp,  Uat_pres  ) 
IF  INT<L/Z)=L/Z  THEN 

LINE  TYPE  Z 
ELSE 

LINE  TYPE  1 
END  IF 

PLOT  X.Atten 
NEXT  X 
LDIR  0 
LORG  Z 
CSIZE  3, .6 
MOVE  X.Atten 
LINE  TYPE  1 

LABEL  USING  "K" ; Freq, "Hz" 
Freq=DROUND<Freq*Z*(  1/3) ,3) 
J=J+1 1 
NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 
PAUSE 
GCLEAR 

OUTPUT  KBD;"K"; 
ELSE 

GCLEAR 

OUTPUT  KBD;"K"; 
END  IF 

Temp=Temp_con 
Freq=Freq_con 
SUBEND 
! 
I 

sub  Att_temp(  Atom_pres ,Freq, Temp) 
I 

I 

THIS  SUBPROGRAM  PLOTS  THE  ATTENUATION  OF  SOUND  VERSUS  UATER- 
VAPOR  PRESSURE  FOR  FIVE  DIFFERENT  TEMPERATURES  AT  10  DEGREE 
INTERVALS  AROUND  THE  INPUT  TEMPERATURE. 


INPUT    ATOMOSPHERIC  PRESSURE  IN  MILLIBARS 

FREQUENCY  OF  ACOUSTIC  ENERGY  IN  HERTZ 
TEMPERATURE  IN  CELSIUS 

OUTPUT   PLOTS  OF  ATTENUATION  VERSUS  WATER-VAPOR  PRESSURE 

VARIABLES 

Atom_pres   INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

Atten      ATTENUATION  OF  ACOUSTIC  WAVE.  CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 

Att_temp    VARIABLE  USED  TO  DETERMINE  IF  OPERATOR  UANTS 
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Es 

Freq 

I 

J 

L 

Temp 

Temp_con 

Wat_pres 

X 
X$ 

Ynax 


TO  PLOT  ATTENUATION  VERSUS  UATER_VAPOR  PRESSURE 

FOR  VARIOUS  TEMPERATURE  AROUND  THE  INPUT  TEMPERATURE. 

THE  PLOT  IS  DONE  IN  THE  SUBPROGRAM  Att_tenp. 
SATURATION  VAPOR  PRESSURE  AT  GIVEN  TEMPERATURE. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 

MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POUER  UNTIL 
FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 
SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 
INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
USED  TO  INSURE  VALUE  OF  TEMPERATURE  PASSED  TO  SUBPROGRAM 

Att_temp  UAS  NOT  CHANGED. 
ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 
STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 
LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 


Yroax=0 

Temp_con=Tenp 

Temp=Tenp+Z® 

FOR   J=0   TO   4    STEP    . 1 

Wai;_pres=J 

CALL  Attenuat ion(  Atom_pres , Atten.Freq, Temp,Uat_pres ) 

IF  Atten>Ynax  THEN 
Ymax=Atten 

END  IF 
NEXT  J 

Ymax=PROUND< Ynax  + . 005 , -  2  ) 
GRAPHICS  ON 
VIEWPORT  15, 1Z0, 10,70 
UINDOU  0,15,0,Ymax 
AXES  .5,. 005, 0,0, 2,2 
CLIP  OFF 
! LABEL  PLOT 
CSIZE  4..B 
LDIR  0 

MOVE  7.5,Yi»iax*1.  1 
LORG  5 

LABEL  "ATTENUATION  VERSUS  UATER-VAPOR  PRESSURE" 
CSIZE  4..B 

! LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  B 
FOR  J=0  TO  15 

MOVE  J,0 

LABEL  J 
NEXT  J 
MOVE  7.5,-Ymax*. 1 
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CSIZE  4, .6 

LABEL  " WATER- VAPOR  PRESSURE   nib" 

! LABEL  VERTICAL  AXES 

LORG  8 


STEP  .01 


CSIZE  4..G 

FOR  J=0  TO 

Ynax 

MOVE  . 

,Z5,J 

LABEL 

J 

NEXT  J 

LDIR  PI/2 

LORG  B 

MOVE  -Z,Ymax«.5 

CSIZE  4..G 

LABEL  "ATTENUATION    1/m" 

! PRINT  INPUT  CONDITIONS 

LDIR  0 

LORG  1 

MOVE  G.Ynax 

CSIZE  4, .6 

LABEL  USING  "K"  ;  "ATMOSPHERIC  PRESSURE  =  "  ,  Atom_pres  ,  "  nb" 

MOVE  E,Ynax*.95 

CSIZE  4, .6 

LABEL  USING  "K" ; "FREQUENCY=" ,Freq, "  HZ" 

CLIP  ON 

J=0 

FOR  L=1  TO  5 

FOR  X=0.  TO  G+J  STEP  .05 
Uat_pres=X 
IF  INT(L/2)-L/Z  THEN 

LINE  TYPE  3 
ELSE 

LINE  TYPE  1 
END  IF 

CALL  Attenuation*  Atom_pres , At ten ,Freq , Tenp , Uat_pres ) 
PLOT  X.Atten 

NEXT  X 

LDIR  0 

LORG  Z 

CSIZE  3, .6 

MOVE  X.Atten 

LINE  TYPE  1 

LABEL  USING  "K" ; Tenp , "C" 

Temp=Temp- 10 

J=J+1 .6 
NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 
PAUSE 
GCLEAR 

OUTPUT  KBD;"K" ; 
FOR  J=0  TO  15 
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PRINT  "  " 
NEXT  J 
! 

PRINT  "UOULD  YOU  ARE  LIKE  TO  SEE  THE  SAME  PLOT  ONLY" 
PRINT  "WITH  ATTENUATION  VERSUS  RELATIVE  HUMIDITY" 
LINPUT  "IF  SO  ENTER  Y,  IF  NOT  ENTER  N",X$ 
Att_-tenp=FNYes(X$) 
IF  Att_temp=1  THEN 

OUTPUT  KBD;"K"; 

Temp=Temp_con 

UINDOU  0,  100,0,Y««iax 

AXES  10,.005,0,0,Z,Z 

CLIP  OFF 

(LABEL  PLOT 

CSIZE  4, .6 

LDIR  0 

MOVE  50, Ymax* 1 . 1 

LORG  5 

LABEL  "ATTENUATION  VERSUS  RELATIVE  HUMIDITY" 

! LABEL  HORIZONTAL  AXIS 

LDIR  0 

LORG  G 

FOR  J=0  TO  100  STEP  10 
MOVE  J,0 
LABEL  J 

NEXT  J 

MOVE  50, -Ymax*. 1 

CSIZE  4, .6 

LABEL  "RELATIVE  HUMIDITY" 

! LABEL  VERTICAL  AXES 

LORG  8 

CSIZE  4, .6 

FOR  J=0  TO  Ymax  STEP  .01 
MOVE  Z,J 
LABEL  J 

NEXT  J 

LDIR  PI/Z 

LORG  G 

MOVE  -10, Ymax*. 5 

CSIZE  4, .6 

LABEL  "ATTENUATION   1/n" 

!PRINT  INPUT  CONDITIONS 

LDIR  0 

LORG  G 

MOVE  50, Ymax 

CSIZE  4, .6 

LABEL  USING  "K" ; "ATMOSPHERIC  PRESSURE=" , Atom_pres , "  mb' 

MOVE  50, Ymax*. 95 

CSIZE  4, .6 

LABEL  USING  "K" ; "FREQUENCY-"  ,  Freq,  "  Hz" 
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CLIP  ON 
J=0 

Tenp=Tei»ip-Z0 
FOR  L=1  TO  5 
T=Temp+Z73 
£5-10^(9.4-2353/1) 
FOR  X»0  TO  90- J  STEP  .5 

Uat_pres  =  Es*A-ton_pres*X/<  (  fliom_pres-Es»(  1 +X/ 1 00) >* 100 ) 
IF  INT(L/Z)=L/Z  THEN 

LINE  TYPE  3 
ELSE 

LINE  TYPE  1 
END  IF 

CALL  Attenuation* Atom_pres ,  Atten.Freq ,  Temp , Uat_pres ) 
PLOT  X.Atten 
NEXT  X 
LDIR  0 
LORG  Z 
CSIZE  3,. 6 
MOVE  X.Atten 
LINE  TYPE  1 

LABEL  USING  "K" ;Temp,"  C" 
Temp=Tenp+ 10 
J=J+15 
NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 
PAUSE 
GCLEAR 

OUTPUT  KBD;"K"; 
ELSE 

OUTPUT  KBD;"K"; 
END  IF 

Tenp=Tenp_con 
SUBEND 


sub  Log_log<  Xmin.Xnax , Ynin, Ymax ,Xvar(  * ) , J, Range (  * ) ,L , Tit le$ ,Xlabel$ , Y label 


THIS  SUBROUTINE  MAKES  A  LOG-LOG  PLOT  OF  DATA  PASSED  FROM  THE 
MAIN  PROGRAM. 

INPUT    MINIMUM  VALUES  OF  X  AND  Y  FOR  PLOT 

MAXIMUM  VALUES  OF  X  AND  Y  FOR  PLOT 

X  AND  Y  VALUES  TO  BE  PLOTTED 

TITLE  OF  PLOT 

LABELS  FOR  X  AND  Y  AXIS 
!      OUTPUT   LOG-LOG  PLOT 
I 
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VARIABLES 

Dx  STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

Dy  STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Titles  STRING  PASSED  TO  SUBPROGRAM  PT  FOR  TITLE  OF  PLOT. 

XlabelS  LABEL  ON  X  AXIS  PASSED  TO  SUBPROGRAM  Pi  FOR  PLOTTING. 

Xmax  VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xmin  VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 

Xrange  VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 

Xvar<*>  VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Pt 

YlabelS  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Pt  FOR  PLOTTING. 

Ymax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Ynin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

Dx»1 

Dy=1 

Xrange=ABS( Xmax-Xwin ) 

Yrange=ABS( Ymax-Ymin) 

UINDOU  Xmn.Xmax,  Ynin.Ynax 

AXES  Dx,Dy,Xmin,Ynin,  1  ,  1 

CLIP  OFF 

! 

! LABEL  PLOT 

CSIZE  4, .6 

LDIR  0 

LORG  5 

MOVE  Xmin+.5*Xrange, 1 . 1* Yrange 

Title$=TRIM$(Title$) 

LABEL  Title« 

! LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  5 

FOR  J=Xmin  TO  Xnax  STEP  Dx 

CSIZE  4, .6 

MOVE  J-.013*Xrange,-.067*Yrange 

LABEL  "10" 

MOVE  J+.013*Xrange,-.033»Yrange 

CSIZE  Z 

LABEL  J 
NEXT  J 

MOVE  Xnin+. 5* Xrange,-. 1Z*Yrange 
LORG  5 
CSIZE  4, .6 

Xlabel$=TRIM$< XlabelS) 
LABEL  XlabelS 
'LABEL  VERTICAL  AXES 
LORG  8 
FOR  J=0  TO  Ynax  STEP  Dy 
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CSIZE  4, .6 

MOVE  Xmin-.0Z5*Xrange,J 

LABEL  "10" 

CSIZE  Z 

MOVE   Xniin-  .00Z5*Xrange ,  J+.03*Yrange 

LABEL    J 
NEXT   J 
LDIR  PI/Z 
LORG   B 

MOVE   Xroin-. 1*Xrange, .5*Yrange 
CSIZE   4, .6 

Ylabel$=TRIM$(Ylabel«> 
LABEL  Ylabel$ 
CLIP  ON 
FOR  J-1  TO  L 

PLOT  LGT(Xvar(  J)  )  ,LGT(Range(  J)) 
NEXT  J 
SUBEND 
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APPENDIX    B 
COMPUTER   PROGRAM    2 

FULLER,    ROBERT  PR0G_Z  FREQUENCY 

9    SEP    85 

*******PURPOSE******* 

THIS  PROGRAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  FREQUENCY.   THE  FOLLOWING  INPUTS 
ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1 )ATOMOSPHERIC  PRESSURE  in  millibars 

Z)SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 
a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 
3)TEMPERATURE  IN  DEGREES  CELSIUS 
4)UATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNDER  DATA 

5)ANTENNA  DIAMETER  IN  METERS 

G)POUER  TRANSMITTED  BY  ECHOSOUNDER  IN  WATTS 

7)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOUING  GRAPHS  TO  AN  EXTERNAL  PLOTTER: 
1 )RANGE(n. )  VERSUS  FREQUENCY  OF  ECHOSOUNDER 
2)RANGE(m.)  VERSUS  EXCESS  ATTENUATION  FOR  VARIOUS  FREQUENCUIES 

*****  VARIABLES***** 


Again      TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 
Ano_change  TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 

ANOUTHER  RUN. 
Ant_area   ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 

ANTENNA  DIAMETER. 
Ant_diam    INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 
Aton_pres   INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 
Atten      ATTENUATION  OF  ACOUSTIC  WAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
Bn         BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 

VERTICAL  VELOCITIES. 
C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

CI         VARIABLE  USED  TO  SCALE  INVERSION  HEIGHT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 
C3         VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 
CneZ<*)     ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 
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!  PARAMETER.   CALCULATED  BASED  ON  SELECTION  OF 

!  PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 

!  AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 

!  PARAMETER. 

!   CteZ       VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

!  THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

!   CveZ       VALUE  OF  THE  VELOCITY  STRUCTURE   PARAMETER. 

!  VALUES  BASED  ON  CALCULATION  USING  ASSUMED  DISSAPATION 

!  RATE. 

!   Dx         STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

!   Dy         STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

Epsilon    DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

Et         TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 

Er         EFFICIENCY  OF  ECHOSOUNDER  UHEN  RECEIVING  ACOUSTIC 
BACKSCATTER. 

Exc_att<*,*>  EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 

F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

Fmax       FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 
ATTENUATION. 

Freq       INPUT  FREQUENCY  OF  ECHOSOUNDER. 

G  ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 

POWER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 

IND        INDEX  USED  FOR  LOOP  FOR  DIFFERENT  FREQUENCIES 

Interval    MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 

Inver      HEIGHT  OF  INVERSION  LAYER. 

J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  UAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Mess_up    TEST  VALUE  TO  CHECK  RESPONSES  UHEN  OPERATOR  IS  INPUTTING 
RESPONSES. 

N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 

Neu_va     TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 
A  VARIABLE  BEFORE  A  NEW  RUN. 

Noise      ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 

Pou_back<*)  POWER  BACKSCATTERED  FROM  GIVEN  RANGE. 

Pou_ret<«)  POWER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

Pou_trans   INPUT  OF  POWER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 

Profile    OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 

Pstar      VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

Pulse      TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.   INPUT  BY 
OPERATOR. 

R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1         VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE 

Range**)    ARRAY  OF  RANGE  VALUES. 

Remainder   REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 
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NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 

ATTENUATION. 
Rge        VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 
Rho        CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 
Sigma      FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 
Speed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 
Sumpou_back  SUM  OF  BACKSCATTERED  ENERGY 
T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
TempS      VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Title$     STRING  PASSED  TO  SUBPROGRAM  Semi_log  FOR  TITLE  OF  PLOT. 
Tstar      INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
Var        USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  WISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 
Ua-fc_pres    ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 
XlabelS    LABEL  ON  X  AXIS  PASSED  TO  SUBPROGRAM  Semi_log  FOR  PLOTTING. 
Xmax       VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xmin       VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xrange     VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 
Xvar<»)     VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Semi. 

YlabelS  LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Semi_log  FOR  PLOTTING. 

Ymax  LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Ymin  SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 

Yrange  RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 

Ze  INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


DECLARE  VARIABLES 
INTEGER  I 

DIM   Pou_back< 1 500) ,Pou_ret<  1500),Range<  1500) 
DIM   CneZ(  1500) ,Exc_att(  1500,30) ,Xvar< 1500) ,Ran<  30) 
DIM   Title$[50] ,Xlabel$[50] ,Ylabel$[ 1G1 
! 

PLOTTER  IS  705,"HPGL" 

LINE  TYPE  1 

! 

(INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" ,Temp 

INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" , At om_pres 

INPUT  "ENTER  UATER  VAPOR  PRESSURE  IN  MILLIBARS" ,Uat_pres 

I 
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! INPUT  ECHOSOUNDER  DATA 

INPUT  "ENTER  ANTENNA  DIAM  IN  METERS" , Ant_di am 

INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS-" .Pulse 

INPUT  "ENTER  POUER  TO  TRANSMITTER  IN  UATTS" ,Pou_trans 

Et«.ZS        ! TRANSMIT  EFFICIENCY 

Er=.Z5        (RECEIVER  EFFICIENCY 

G=.40        ! ANTENNA  EFFECTIVE  APERATURE  FACTOR 


SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 
PRINT  "YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 
PRINT  "PROFILE  TO  BE  USED." 
PRINT  "YOUR  SELECTIONS  ARE" 

PRINT  "      1       A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS 
PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  UHICH  ' 

PRINT  "  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE- 

PRINT  "  " 

PRINT  "      Z      THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 
PRINT  "  UITH  A  FACTOR  OF  TUO  TO  APPROXIMATE  LOOKING" 

PRINT  "  UP  A  THERMAL  PLOOM" 

PRINT  "  " 

PRINT  "      3      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS' 
PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 

PRINT  "  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 

PRINT  "  IS  EXP(-.001»HEIGHT  ABOVE  65  METERS  AND  HEIGHT" 

PRINT  "  TO  THE  -1.4G  BELOU  B5  METERS- 

PRINT  "  " 

PRINT  "     4      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
PRINT  "  HEIGHT  TO  THE  -4/3  " 

! 

PRINT  "  " 
PRINT  "  " 
Mess_up= 1 
UHILE  Mess_up=1 

INPUT  "ENTER  THE  DESIRED  PROFILE  (  1  OR  Z  OR  3  OR  4)", Prof  lie 
IF  Profile=4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" , Inver 
END  IF 
IF  Profiled  THEN 

Mess_up=0 
ELSE 

IF  Profile=Z  THEN 

Mess_up=0 
ELSE 

IF  Profile=3  THEN 

Mess_up=0 
ELSE 

IF  Profile=4  THEN 

Mess_up=0 
ELSE 
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PRINT  Profile,"  WAS  NOT  ONE  OF  THE  OPTIONS!!!!" 
Mess_up=1 
END  IF 
END  IF 
END  IF 
END  IF 
END  WHILE 
OUTPUT  KBD;"K" ; 
I 
! 

Again=1 
WHILE  Again=1 
FOR  J-1  TO  I 

Pou_back< J>=0 
Pou_ret(  J)=0 
Range < J>=0 
CneZ( J>=0 
FOR  L=1  TO  Ind 

Exc_att<J,L)=0 
NEXT  L 
NEXT  J 
! 

(CONVERT  TEMPERATURE  TO  KELVIN 
T=Tenp+Z73 
I 

! CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE* CELCIUS) 

Speed_sound=Z0.  05*  (  T> *  .5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 

C=Z0.05*Z73".S 

I 

1=1     ! INDEX  FOR  LOOP 

Ind-1 

Pou_back<0>=0     (INITIALIZE  VARIABLE  FOR  POWER  BACKSCATERRED 

An-t_area=PI*(An-t_diam/ZrZ    (ANTENNA  AREA 

Interval  =  ( Speed_sound*Pulse* 1  .E-3)/Z 

! 

(INITIALIZE  PLOTTING  PARAMETERS  SO  CAN  COMPARE  DURING  RUN  TO  SCALE 

! PLOTS 

Ymin=1000 

Ynax=0 

(CALCULATE  RETURNED  SIGNAL  POWER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 

(FOR  EACH  FREQUENCY. 

Range(0)=0 

FOR  Freq=100  TO  3000  STEP  100 
Bn=Z*Freq*< 1  -  1 /(  6/Speed_sound+ 1 )) (BANDWIDTH  FOR  3M/S  VERTICAL  VELOCITY 
Noise=1 .38E-Z3*Bn*(T)+Z.E-14! MINIMUM  DETECTABLE  SIGNAL  AS 
!      JOHNSON  NOISE(NEGLIGABLE)+  ESTIMATED  BACKGROUND 
(CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
CALL  Attenuai ion< Atom_pres , Atten.Freq.Temp, Wat_pres ) 
K=Z*PI*Freq/Speed_sound     ! WAVENUMBER 
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1-1 

Supipou_back  =0 
REPEAT 

Range( I )=Range< 1-1 )+Z 
SELECT  Profile 
CASE  1 

•THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
ITULAROSA  BASIN,  NEU  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 
! AVERAGING  TIME. 

CteZ=Z. 1Z*Range( I)"<-1.16) 
CASE  Z 

ITHIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
ITULAROSA  BASIN,  NEU  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1.8  WAS  INCLUDED  AS  PER  DIFFERENT 
•AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  UAS 
•INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 

CteZ=Z*Z. 1Z*Range( X)*(-1.16) 
CASE  3 

ITHIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
I  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
I  NIGHT  CONDITIONS.   THE  DEPENDENCE  WITH  HEIGHT  " 
IIS  EXP(-.001*HEIGHT  ABOVE  E5  METERS  AND  HEIGHT" 
I  TO  THE  -1.4G  BELOW  65  METERS" 
IF  Range(I)<65  THEN 

CteZ=7S.5*Range(I)A(-1 . 4G ) 
ELSE 

C-teZ  =  3.66E-Z*EXP(-.001»Range(I)) 
END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
I     THAT  IT'S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POWER 
I     THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  Cm. /sec 
I  THE  DEPENDENCE  WITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
I  MY  LEFT  EAR 

•EQUATION  FROM  NEFF ,  1975 
SELECT  Range<I)/Inver 
CASE  <.9 

C3=(  (  .0Z4)*(Tr<  ,6G7)> 
CteZ=C3*(  Ranged  ))A(-1  .33) 
CASE  .9  TO  1 
CteZ=CteZ 
C1=CteZ 
R1=Range( I ) 
CASE  1  TO  1.3 

CteZ=10A((LGT(C3*Range( I ) )-LGT(C3*R1  )  )+LGT(C1  ) ) 
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RZ=Range(  I ) 

CZ=CteZ 
CASE  ELSE 

CteZ=C3*(Range(ir<-1.33)-RZA<-1  .33))  +  CZ 
END  SELECT 
END  SELECT 

! CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epsilon=< .Z8EE/Range< I ) )*( 1+.07*( Range< I )  )A( .6) )*( 1 .5) 
CveZ=Z*Epsilon"( .6E7) 

! CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 
! FORMULA  FROM  TETARSKI 
CneZ(I)  =  (CteZ/(Z.38E-i5))  +  (CveZ/(C*C)  ) 

! CALCULATE  THE  FRACTION  OF  POUER  BACKSCATTERED  FROM  INTERVAL. 
! EQUATION  FROM  NEFF,  1975 
Signa=( .0039*<K*<  1/3) )*CteZ  )/( T)"Z 
I 

"CALCULATE  THE  EXCESS  ATTENUATION.   EQUATION  FROM  CLIFFORD,  1980 
!THE  EXCESS  ATTENUATION  IS  Exc_att( I ,L ,X ) 
IF  K50  THEN 

Renainder=0 
ELSE 

IF  K150  THEN 

Remainder=Range( I )  MODULO  10 
ELSE 

Remainder=Range< I )  MODULO  Z0 
END  IF 
END  IF 

IF  Renainder=0  THEN 
Rho=0 
L=0 
H=0 

Rge=Range(I)     ! CONSTANT  IN  INTEGRAL 
R=0 
FOR  J=0  TO  Z*I 

F=CneZ( INT< J/Z+1 ) ) 

F=F*( 1-R/(Rge)"( 1 . G7)+( R/( Rge ) ) A ( 1 .67)) 

IF  J>0  THEN 

IF  J<Z*I  THEN 

IF  INT<J/Z)=J/Z  THEN 
L  =  L  +  F 
F=0 
ELSE 

H=H+F 
F=0 
END  IF 
END  IF 
END  IF 
Rho=Rho+F 
R=R+1 
NEXT  J 
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Rho=Rho+4*L+Z*H 

Rho=< <<Rho*.33)»K»K«1 .46)"  (-.6) ) 

N=(Ant_diani/Rho>AZ 

IF    N<=1    THEN 

Ze=1/( 1+N) 
ELSE 

Ze-1.5/(1+N) 

ISTEP  OF  1.5«=>SEE  CLIFFORD  1980 
END  IF 
Exc  att<I,Ind)=Ze*Ze 


ELSE 


Exc_att(  I  ,Ind)=Exc_att<  1-1  ,Ind) 
END  IF 
! 

(CALCULATE    THE   POUER   BACKSCATTERED 

Pou_back( I )=(Pou_trans«Et-Sumpou_back )*EXP( -At  ten* Range ( I ) ) 
Pou_back(  I )=Pou_back( I )*Interval*Exc_att(  I , Ind)»Sigma 
Supipou_back=Sunpou_back+Pou_back  (  I  ) 
(CALCULATE    THE    POUER   BACKSCATTERED   TO   THE    ANTENNA 

Pou_ret( I )-Pou_back( I )*EXP( -Atten*Range(  I ) )*Ant_area»G*Er/Range< I ) 
PRINT    "RAN6E  =  "  .Ranged) 
PRINT    "POUER   RETURNED=",Pou_ret(  I ) 
1=1  +  1 
UNTIL   Pou_ret(I-1  >>0   AND  Pou_ret(  I- 1 )<=Noise 
Ran< Ind)=Range< I-Z ) 
IF    Range<I-Z)>Ynax    THEN 

Ymax=Range( I-Z ) 
END    IF 
IF    Range(  I-ZKYmin    THEN 

Ymin=Range( I-Z ) 
END   IF 
Ind=Ind+1 
Again= 1 
FOR   J-1    TO    I 

Pou_back< J)=0 
Pou_ret< J)=0 
Range<J)=0 
CneZ< J)=0 
NEXT   J 

PRINT    "FREOUENCY=M ,Freq 
NEXT   Freq 
! 
! 

OUTPUT    KBD;"K"; 
PRINT    "INPUT    CONDITIONS" 
! 

PRINT   "    " 

PRINT  USING  "K";"TEMPERATURE=  " , Temp , "CELCIUS" 

PRINT  USING  "K"; "ATMOSPHERIC  PRESSURE=    " , Atom_pres , "nb" 
PRINT  USING  "K";"UATER  VAPOR  PRESSURE=    " , Uat_pres , "mb" 
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PRINT  USING  "K'VPULSE  LENGTH= 
PRINT  USING  "K";" ANTENNA  DIAMETER= 
PRINT  USING  "K'VPOUER  TRANSMITTED= 
PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED 
IF  Profile=4  THEN 

PRINT  "INVERSION  HEIGHT  ".Inver 
END  IF 

SELECT  Profile 
CASE  1 

PRINT 

PRINT 

PRINT 
CASE   Z 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 
CASE    3 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 
CASE    4 

PRINT 

PRINT 
END  SELECT 
! 

PRINT    USING 


.Pulse , "ns" 
,  Atrt_diam, " 
,Pou_trans , ' 
".Profile 


m . 


UATTS" 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 
WHICH  GIVES  A  HEIGHT  TO  THE  -1.1E  PROFILE" 

A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 
BUT  WITH  A  FACTOR  OF  TWO  TO  APPROXIMATE" 
LOOKING  UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398" 
FOR  NIGHT  CONDITIONS.   THE  DEPENDENCE  WITH" 
HEIGHT  IS  EXP<-.001»HEIGHT  ABOVE  B5  METERS" 
AND  HEIGHT  TO  THE  -  1 . 4G  BELOW  65  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
HEIGHT  TO  THE  -4/3  " 


"K'VTHE  MINIMUM  DETECTABLE  SIGNAL  WAS  SET  AT"  ,  Noise ,  "WATTS' 


PRINT  "  " 

I 

! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBD;"K"; 

I 

GRAPHICS  ON 

VIEWPORT  15,  1Z0, 10,70 

I 

! 

IPLOT  RANGE  VERSUS  FREQUENCY 

Xmin=0 

Xmax=3000 

Ymin=PROUND<  Ymin-50,Z ) 

Ymax=PROUND<  Ymax+50,Z ) 

Dx=500 

Dy= 1 00 

Xrange= Xmax - Xm l n 
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Yrange=Ynax-Ymin 

UINDOU  Xmin.Xmax.Ymin.Ymax 

AXES  Dx,DyfXnin,Ymin,Z,1 

CLIP  OFF 

! 

! LABEL  PLOT 

CSIZE  4..B 

LDIR  0 

LORG  5 

MOVE  .5*Xrange,Ymax*1 .05 

LABEL  "RANGE  VERSUS  FREQUENCY" 

(LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  5 

FOR  J-0  TO  Xnax  STEP  Dx 

CSIZE  4, .6 

MOVE  J,Ynin-.05*Yrange 

LABEL  J 
NEXT  J 

MOVE  . 5*Xrange,Ymin- . 1*Yrange 
CSIZE  4, .6 
LABEL  "FREQUENCY" 
•LABEL  VERTICAL  AXES 
LORG  8 
FOR  J=Ymin  TO  Ymax  STEP  Dy 

CSIZE  4, .6 

MOVE  Xmn-.00Z5»Xrange,J 

LABEL  USING  "K";J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  Xmin-. 15*Xrange, . 5* Yrange+Ymin 
CSIZE  4, .6 

LABEL  "RANGE  (METERS)" 
CLIP  ON 
Ind=1 
FOR  Freq=100  TO  3000  STEP  100 

PLOT  Freq,Ran< Ind) 

Ind=Ind+1 
NEXT  Freq 
! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBD;"K"; 


PLOT    RANGE    VERSUS    EXCESS    ATTENUATION 
Xcun=0 
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Ymax=4 

Xi"iax=  1 

Ymin=0 

Range(0)=0 

Title$=" RANGE  vs  EXCESS  ATTENUATION" 

Xlabel$="EXCESS  ATTENUATION" 

Ylabel$=" RANGE  (METERS)" 

Dx  =  .1 

Dy=1 

Xrange=ABS( Xmax-Xnin ) 

Yrange=ABS( Ymax-Ynin) 

UINDOU  Xmin.Xmax.Ymin.Ymax 

AXES  Dx,Dy,Xi»iin,Yeiin, 1 , 1 

CLIP  OFF 

! 

! LABEL  PLOT 

CSIZE  4..G 

LDIR  0 

LORG  5 

MOVE  Xnin+.5*Xrange, 1 . 1*Yrange 

Title$=TRIM$<Title$> 

LABEL  TitleS 

! LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  5 

FOR  J=Xmin  TO  Xmax  STEP  Dx 

MOVE  J, -.033* Yrange 

CSIZE  4..G 

LABEL  J 
NEXT  J 

MOVE  Xnin+.5*Xrange,-. 1Z*Yrange 
LORG  5 
CSIZE  4, .6 

Xlabel$=TRIM$<  XlabelS) 
LABEL  XlabelJ 
! LABEL  VERTICAL  AXES 
LORG  8 
FOR  J=0  TO  Ymax  STEP  Dy 

CSIZE  4, .6 

MOVE  Xmin-.0Z5*Xrange,J 

LABEL  "10" 

CSIZE  Z 

MOVE  Xrcin- . 00Z5*Xrange , J+ . 03* Yrange 

LABEL  J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  Xnin- . 1 *Xrange , .5* Yrange 
CSIZE  4, .6 
Ylabel$=TRIM$(  YlabelS) 
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LABEL  YlabelS 
FOR  L=5  TO  Ind-1  STEP  5 
X-1 
UHILE  Exc_att<X,L)>0 

Xvar(X)=Exc_att(X,L> 
Range(  X ) «Range<  X-1  )  +Z 
X=X+1 
END  UHILE 
IF  INT(L/2)-L/Z  THEN 

LINE  TYPE  1 
ELSE 

LINE  TYPE  3 
END  IF 
Freq=L*100 

CALL  Semi_log( Freq,Xvar(  * ) ,X-  1  ,Range(  * ) ,X- 
NEXT  L 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 
PAUSE 


1  ) 


KBD;"K"; 

"WOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?' 


GCLEAR 
OUTPUT 
L INPUT 
Again=FNYes<X$) 
Mess_up=1 
UHILE  Mess_up=l 
SELECT  Again 
CASE  1 

LINPUT  "WOULD  YOU  LIKE 
Neu_va=FNYes(X$> 
Ano_change= 1 
UHILE  Ano_change=1 
SELECT  Neu_va 
CASE  1 
PRINT  " 
PRINT  USING 
PRINT  USING 
PRINT  USING 
PRINT  "  " 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


X$ 


TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)7",X$ 


VARIABLE 


USING 
USING 
USING 
USING 


"K" 

"K" 
"K" 

"K" 
"K" 
"K" 
"K" 


"1 
"2 

"3 

"4 
"5 
"G 

"7 


TEMPERATURE 
ATOMOSPHERIC 


CURRENT  VALUE" 

", Temp, "CELSIUS" 


PRESSURE 


WATER  VAPOR  PRESSURE 


" , Aton_pres , " mb' 
, Wat_pres , "mb" 


ANTENNA  DIAMETER     " , Ant_diam , "  n." 
PULSE  LENGTH  ", Pulse,"  ns" 

POWER  TRANSMITTED  " ,Pou_trans , "  WATTS' 
ATMOSPHERIC  PROFILE   ".Profile 


"ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 
"CHANGE" 
Var 
SELECT  Var 

CASE  1 

INPUT 
CASE  Z 


'TEMPERATURE=" ,Temp 
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"ATMOSPHERIC  PRESSURE  IN  nb=" , Atom_pres 
"UATER  VAPOR  PRESSURE  IN  nb=" ,Uat_pres 

IN  n. =" , An-t_dian 
"PULSE  LENGTH  IN  pis=", Pulse 
"POUER  TRANSMITTED  IN  UATTS=" . Pou_trans 


INPUT 
CASE  3 

INPUT 
CASE  4 

INPUT  "ANTENNA  DIAMETER 
CASE  5 

INPUT 
CASE  G 

INPUT 
CASE  7 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


3   ==>  Ct"Z 


Ct"Z  PROFILE  OF  Z*(-l. IB)" 

FROM  UALTERS/KUNDEL  1981" 
SAME  AS  ONE  BUT  WITH  FACTOR" 
OF  TUO  TO  APPROXIMATE  LOOKING" 
UP  A  THERMAL  PLOOM" 

PROFILE  OF  EXP(-.001*Z>" 
FROM  UALTERS/KUNDEL  1981." 
CT'Z  PROFILE  OF  ZA<-4/3)" 
ENTER  NUMBER  OF  DESIRED  PROFILE" , Prof ile 
IF  Profile=4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS=" , Inver 
END  IF 
CASE  ELSE 

PRINT  Var,"IS  NOT  ONE  OF  THE  OPTIONS" 
END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE< ENTER  Y  OR  N)?",X$ 
Neu_va=FNYes(X$) 
Mess_up=Z 
CASE  Z 

Mess_up=Z 
Ano_change= Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  !  M  !  ! " 
Ano_change= 1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",XS 
Neu_va=FNYes(X$) 
END  SELECT 
END  UHILE 
CASE  Z 

Mess_up=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N   !!!!!!" 
Mes5_up= 1 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again=FNYes(X$) 
SELECT 


END 
END  UHILE 
END  UHILE 
PRINT  "THAT'S  ALL 
end 


FOLKS" 
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CALCULTE  THE  ATTENUATION 
sub  Attenuation( Atom_pres ,  At ten, Freq, Tenp ,Uat_pres ) 
THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 

INPUT   ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
FREQUENCY  OF  SOUND  UAVE  IN  HERTZ 
TEMPERATURE  IN  DEGREES  CELCIUS 
UATER-VAPOR  PRESSURE  IN  MILLIBARS 

OUTPUT   ATTENUATION  IN  1 /METERS 


VARIABLES 
Aton_pres 
Atten 

Att  max 


INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

ATTENUATION  OF  ACOUSTIC  UAVE.   CALCULATED  IN 
SUBPROGRAM  ATTENUATION. 

VARIABLE  IN  SUBPROGRAM  ATTENUATION.   IT  IS  THE 
ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 
ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 


Fmax 

Freq 

H 
Pstar 


RATIO  OF  FREQUENCY  TO  FREQUENCY 
FREQUENCY  OF  MAXIMUM  ATTENUATION. 

ATTENUATION. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 
VARIABLE  USED  IN  INTEGRATION  FOR 


AT  MAXIMUM  ATTENUATION. 
USED  IN  SUBPROGRAM 


EXCESS  ATTENUATION. 


Temp 
Tstar 

Uat_pres 


VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 
IN  SUBPROGRAM  ATTENUATION. 

ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 
OPERATOR. 
H= 100*Uat_pres/Atom_pres 
Tstar=( 1 .8*Temp+49Z>/519 
Pstar=Aton_pres/1014 

Fi»iax  =  <  10+6G00*H+44400*H»H)»Pstar/Tstar".8 
Att_i*iax=.0078»Fi»iax*Tstar"<-Z.5>*EXP<7.77*<  1-1/Tstar  )  ) 
F=Freq/Fmax 

Atten=(Att_i»iax/304.8>*<  (  .  1  8*F  )  *Z  +  <  Z*F*F/(  1+F«F>  )AZ)A  .5 
Atten=< Atten+1 . 74E- 1 0* Freq* Freq )/4. 35 
SUBEND 
l 

def  FNYes(X$) 

!THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 

! 

!     INPUT     X$ 

! 

!     OUTPUT    FNYes 
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VARIABLES 

TenpS      VARIABLE  STRING  USED  IN  FUNCTION  YES. 
X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 
!  RESPONSE  TO  YES  OR  NO  QUESTION. 

DIM  TempSI 1] 
TempSI 1 ,1]=TRIM$(X$) 
SELECT  Tewp$ 
CASE  "Y","y" 

RETURN  1 
CASE  "N",Mn" 

RETURN  Z 
CASE  "  " 

RETURN  1 
CASE  ELSE 

RETURN  -Z 
END  SELECT 
FNEND 
! 

! 
sub  Seni_log( Freq,Xvar(  * ) ,L , Range ( * ) , J ) 

THIS  SUBROUTINE  MAKES  A  SEMI-LOG  PLOT  OF  DATA  PASSED  FROM  THE 
MAIN  PROGRAM. 

INPUT 

X  AND  Y  VALUES  TO  BE  PLOTTED 
OUTPUT   SEMI-LOG  PLOT 

VARIABLES 

J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Xvar(»)  VALUE  OF  ARRAY  OF  X  VALUES  TO  BE  PLOTTED  IN  SUBPROGRAM  Semi 


CLIP  ON 

FOR  J=1  TO  L  STEP  1 

PLOT  Xvar< J) ,LGT(  Range(  J)) 
NEXT  J 
LDIR  PI/4 
LORG  Z 
CSIZE  3, .6 
LINE  TYPE  1 

MOVE  Xvar( J-1 ),L6T(Range(  J- 1  )) 
LABEL  USING  "K" ; Freq, "Hz" 
SUBEND 
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APPENDIX  C 


r-l 


COMPUTER  PROGRAM  3 


! FULLER,  ROBERT       PR0G_3       EFFICIENCY 

! 

! 10  SEP  85 


*•*••« #PURP0SE* ****** 

!      THIS  PROGRAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
! ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  EFFICIENCY  OF  THE  THE  TRANSDUCER. 
THE  FOLLOWING  INPUTS  ARE  REQUIRED: 

ATMOSPHERIC  DATA 

1 )ATOMOSPHERIC  PRESSURE  in  millibars 

Z)SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 
a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 
!     3)TEMPERATURE  IN  DEGREES  CELSIUS 
!     4)WATER  VAPOR  PRESSURE  IN  millibars 
! 

ECHOSOUNDER  DATA 

5)ANTENNA  DIAMETER  IN  METERS 

6  FREQUENCY 

I     7)P0UER  TRANSMITTED  BY  ECHOSOUNDER  IN  UATTS 
!     8)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 
! 

!THE  PROGRAM  OUTPUTS  THE  FOLLOWING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 
!     1)RANGE(m.)  VERSUS  EFFICIENCY  OF  ECHOSOUNDER 


***** VARIABLES***** 


Again      TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 
Ano_change  TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 

ANOUTHER  RUN. 
Ant_area   ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 

ANTENNA  DIAMETER. 
Ant_diam    INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 
Atom_pres   INPUT  OF  ATMOSPHERIC  PRESSURE  IN  nb. 
Atten      ATTENUATION  OF  ACOUSTIC  WAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
Bn         BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 

VERTICAL  VELOCITIES. 
C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

CI         VARIABLE  USED  TO  SCALE  INVERSION  HEIGHT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 
C3         VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 
CneZ<»)     ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 
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!  PARAMETER.   CALCULATED  BASED  ON  SELECTION  OF 

!  PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 

!  AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 

!  PARAMETER. 

!   CteZ       VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 
!  THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 

!   CveZ       VALUE  OF  THE  VELOCITY  STRUCTURE  PARAMETER. 
!  VALUES  BASED  ON  CALCULATION  USING  ASSUMED  DISSAPATION 

RATE. 

Dx         STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

Dy         STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 

Epsilon    DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 

Eff        TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 

Ef-f        EFFICIENCY  OF  ECHOSOUNDER  WHEN  RECEIVING  ACOUSTIC 
BACKSCATTER. 

Exc_att     EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 

F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 

Fnax       FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 
ATTENUATION. 

Freq       INPUT  FREQUENCY  OF  ECHOSOUNDER. 

G  ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 

POWER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 

IND        INDEX  USED  FOR  LOOP  FOR  DIFFERENT  EFFICIENCIES 

Interval    MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 

Inver      HEIGHT  OF  INVERSION  LAYER. 

J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  WAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Mess_up    TEST  VALUE  TO  CHECK  RESPONSES  WHEN  OPERATOR  IS  INPUTTING 
RESPONSES. 

N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 

Neu_va     TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 
A  VARIABLE  BEFORE  A  NEW  RUN. 

Noise      ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 

Pow_back<*>  POWER  BACKSCATTERED  FROM  GIVEN  RANGE. 

Pou_ret<»)  POWER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 

Pou_trans   INPUT  OF  POWER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 

Profile    OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 

Pstar      VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 
CALCULATION. 

Pulse      TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.   INPUT  BY 
OPERATOR. 

R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1         VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE 

Range**)    ARRAY  OF  RANGE  VALUES. 

Remainder   REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 
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NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 

ATTENUATION. 
Rge        VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 
Rho        CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 
Sigma      FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 
Speed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 
Sumpou_back 

T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
Temp«      VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Tstar      INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
Var         USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  WISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 
Uat_pres    ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 
Xmax       VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xmin       VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xrange     VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 
YlabelS    LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Seni_log  FOR  PLOTTING. 
Ymax       LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Ymin       SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Yrange     RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 
Ze         INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back<  1500) ,Pou_ret( 1500) , Range*  1500) 
DIM  Cne2( 1500) ,Ran( 30) 
! 

PLOTTER  IS  705,,,HPGL•, 
LINE  TYPE  1 


"INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" , Temp 

INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" ,Atom_pres 

INPUT  "ENTER  UATER  VAPOR  PRESSURE  IN  MILLIBARS" , Uat_pres 


INPUT  ECHOSOUNDER  DATA 
INPUT  "ENTER  ANTENNA  DIAM  IN  METERS" , Ant_diam 
INPUT  "ENTER  ECHOSOUNDER  FREQUENCY  IN  HERTZ", Freq 
INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS=" ,Pul se 
INPUT  "ENTER  POWER  TO  TRANSMITTER  IN  UATTS" ,Pou_trans 
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G=.40 


! ANTENNA  EFFECTIVE  APERATURE  FACTOR 


YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 
PROFILE  TO  BE  USED." 
YOUR  SELECTIONS  ARE" 


1 


l 

! SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  ' 
PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  WHICH  " 
GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 
UITH  A  FACTOR  OF  TUO  TO  APPROXIMATE  LOOKING" 
UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
PRESENTED  IN  WALTERS/KUNDEL  1381  PAGE  398  FOR" 
NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
IS  EXP<-.001*HEIGHT  ABOVE  B5  METERS  AND  HEIGHT" 
TO  THE  -1.4G  BELOU  G5  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
HEIGHT  TO  THE  -4/3  " 


PRINT  "  " 
PRINT  "  " 
Mess_up=1 
UHILE  Mess_up=1 

INPUT  "ENTER  THE  DESIRED  PROFILE  (1  OR  Z  OR  3  OR  4)", Prof  lie 
IF  Profile=4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" , Inver 
END  IF 
IF  Profiled  THEN 

Mess_up=0 
ELSE 

IF  Profile=Z  THEN 

Mess_up=0 
ELSE 

IF  Profile=3  THEN 

Mess_up=0 
ELSE 

IF  Profile=4  THEN 

Mess_up=® 
ELSE 

PRINT  Profile,"  UAS  NOT  ONE  OF  THE  OPTIONS!!!! 
Mess_up= 1 
END  IF 
END  IF 
END  IF 
END  IF 
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END  WHILE 

OUTPUT  KBD;"K"; 

! 

! 

Again=1 

WHILE  Again=1 

FOR  J=»1  TO  I 

Pou_back(J)«0 
Pou_ret(  J)=0 
Range< J)=0 
CneZ(J>=0 

NEXT  J 

! 

(CONVERT  TEMPERATURE  TO  KELVIN 

T=Temp+Z73 

! 

(CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE( CELCIUS ) 

Speed_sound=Z0. 05»( T ) " . 5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 

C  =  Z0.05*Z73'\5 

I 

1=1     (INDEX  FOR  INNER  LOOP  (POWER  RETURNED  TO  NOISE) 
Ind-1   (INDEX  FOR  OUTER  LOOP  (EFFICIENCY) 

Pou_back(0>»0     (INITIALIZE  VARIABLE  FOR  POWER  BACKSCATERRED 
Ant_area=PI»(Ant_diam/Z)AZ    (ANTENNA  AREA 
Interval=( Speed_sound* Pulse* 1 .E-3)/Z 
I 

(INITIALIZE  PLOTTING  PARAMETERS  SO  CAN  COMPARE  DURING  RUN  TO  SCALE 
(PLOTS 
Ynin=1000 
Ynax=0 

(CALCULATE  RETURNED  SIGNAL  POWER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 
(FOR  EACH  FREQUENCY. 
Range(0)=0 

FOR  Eff=.05  TO  .5  STEP  .05 
Bn=Z*Freq»< 1-1 /(  G/Speed_sound+ 1 )) (BANDWIDTH  FOR  3M/S  VERTICAL  VELOCITY 
Noise=1 .38E-Z3*Bn*(T)+Z.E-14!MINIMUM  DETECTABLE  SIGNAL  AS 
!      JOHNSON  NOISE(NEGLIGABLE)+  ESTIMATED  BACKGROUND 
(CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
CALL  Attenuat  ion(  Atom_pres , At ten , Freq , Temp , Wat_pres ) 
K=Z*PI*Freq/Speed_sound     ! UAVENUMBER 
1  =  1 

Sumpou_back  =0 
REPEAT 

Range( I )  =  Range(  I- 1 )  +  Z 
SELECT  Profile 
CASE  1 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
(AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
(WAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
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(TULAROSA  BASIN,  NEW  MEXICO. 

(AN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 

(AVERAGING  TIME. 

CteZ=Z.  1Z*Range<ir(-1.16> 
CASE  Z 

!THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  UEATHER  ABOVE  THE 
(TULAROSA  BASIN,  NEW  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 
(AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TUO  UAS 
(INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 

CteZ=Z*Z. 1Z*Range<I)~(-1. IB) 
CASE  3 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
(PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
(NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
(IS  EXP<-.001« HEIGHT  ABOVE  G5  METERS  AND  HEIGHT" 
(TO  THE  -1.4G  BELOU  G5  METERS" 
IF  Range(IXG5  THEN 

CteZ=75.5*Range< I)*(-1 .46) 
ELSE 

CteZ=3.66E-Z*EXP(-.001»Range(I)) 
END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!     THAT  IT'S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POWER 
!     THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  Cm. /sec 
i THE  DEPENDENCE  UITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
(MY  LEFT  EAR 

(EQUATION  FROM  NEFF ,  1975 
SELECT  Range< I )/Inver 
CASE  <.9 

C3=(<  .0Z4)*(Tr(  .GB7)> 
CteZ=C3*(  Ranged)  )A(-1.33) 
CASE  .9  TO  1 
CteZ=CteZ 
C1=CteZ 
R1=Range< I ) 
CASE  1  TO  1.3 

CteZ=10/>((LGT(C3*Range<  I  )  )-LGT<  C3»R1  ))  +  LGT(C1  )  ) 
RZ=Range( I ) 
CZ=CteZ 
CASE  ELSE 

CteZ=C3*(Range(ir<-1  .  33)-RZ"<  -  1  .  33  )  )  +  CZ 
END  SELECT 
END  SELECT 

(CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epsilon=( .Z866/Range< I ) )*( 1+ .07*( Range< I )>A(.E>)A<1.5> 
CveZ=Z*EpsilonA( .667) 
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(CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

! FORMULA  FROM  TETARSKI 

CneZ(I)«(C+.eZ/<Z.98E+5))  +  <CveZ/<C»C>) 

(CALCULATE  THE  FRACTION  OF  POUER  BACKSCATTERED  FROM  INTERVAL. 

(EQUATION  FROM  NEFF ,  1S75 

Sigma=< .0039*(K"(  1/3) >*CteZ >/< T) AZ 

I 

(CALCULATE  THE  EXCESS  ATTENUATION.   EQUATION  FROM  CLIFFORD.  1980 

I  THE  EXCESS  ATTENUATION  IS  Exc_att 

IF  K50  THEN 

Remainder =0 
ELSE 

IF    K150   THEN 

Reroai nder= Range ( I )  MODULO  10 
ELSE 

Remainder=Range< I )  MODULO  Z0 
END  IF 
END  IF 

IF  Renainder-0  THEN 
Rho=0 
L=0 
H=0 

Rge=Range( I )     (CONSTANT  IN  INTEGRAL 
R=0 
FOR  J=0  TO  Z*I 

F=CneZ(INT( J/Z+1 )) 

F=F»(  1-R/(Rger<  1  . G7>+< R/< Rge)  )*(1.67>) 

IF  J>0  THEN 

IF  J<Z«I  THEN 

IF  INT(J/Z)=J/Z  THEN 
L  =  L  +  F 
F=0 
ELSE 

H=H+F 
F=0 
END  IF 
END  IF 
END  IF 
Rho=Rho+F 
R=R+1 
NEXT  J 

Rho=Rho-U*L  +  2*H 

Rho=(  (  (  Rho* . 33  >*K*K* 1 . 46 ) A  < - . 6  > ) 
N=(Ant_diam/Rho)"Z 
IF  N<=1  THEN 

Ze=1/( 1+N) 
ELSE 

Ze=1 .5/( 1+N) 

(STEP  OF  l.5-->SEE  CLIFFORD  1980 
END  IF 
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Exc_att=Ze*Ze 
ELSE 

E x c_at t = E x c_at t 
END  IF 
! 

(CALCULATE  THE  POUER  BACK SCATTERED 

Pou_back(  I  )=<Pou_trans*Ef f-Sumpou_back  )*EXP(  -Atten*Range<  I  )  ) 
Pou_back(  I )=Pou_back( I ) *I nterval * Ex c_att* Sigma 
Sumpou_back  =  Sumpouj_back+Pou_back(  I ) 
! CALCULATE  THE  POUER  BACKSCATTERED  TO  THE  ANTENNA 
Pou_ret(  I  )=Pou_back<  I  )*EXP(  -At -ten*  Range  (  I  )  >»An+,_area*G*Ef  f /Range(  I  ) 
PRINT  "RANGE=",  Ranged) 
PRINT  "POUER  RETURNED-" ,Pou_ret< I) 
1  =  1  +  1 
UNTIL  Pou_ret(I-1 >>0  AND  Pou_ret(  I- 1 )<=Noise 
Ran( Ind)=Range( I-Z) 
IF  Ranged- Z)>Ymax  THEN 

Ymax=Range( I-Z  ) 
END  IF 
IF  Range<  I-ZKYmin  THEN 

Ymin=Range( I-Z ) 
END  IF 
Ind=Ind+1 
Again= I 
FOR  J-1  TO  I 

Pou_back( J)=0 
Pou_rei(  J)=0 
Range( J)=0 
CneZ( J)=0 
NEXT  J 
NEXT  Eff 
! 
! 

OUTPUT  KBD;"K"; 
PRINT  "INPUT  CONDITIONS" 
I 

PRINT  "  " 

PRINT  USING  "K";"TEMPERATURE= 

PRINT  USING  "K"; "ATMOSPHERIC  PRESSURE= 

PRINT  USING  "K" ;"UATER  VAPOR  PRESSURE= 

PRINT  USING  "K"; "PULSE  LENGTH= 

PRINT  USING  "K"; "ANTENNA  DIAMETER= 

PRINT  USING  "K" ; "ECHOSOUNDER  FREOUENCY= 

PRINT  USING  "K"; "POUER  TRANSMITTED= 

PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ".Profile 

IF  Profile=4  THEN 

PRINT  "INVERSION  HEIGHT  ".Inver 
END  IF 
SELECT  Profile 

CASE  1 


,Temp,"CELCIUS" 
, Atom_pres , "mb" 
,  Ua+._pres ,  "mb" 
.Pulse , "ms" 
, Ant_diam,"  m." 
,Freq, "Hz" 
,Pou_trans,"  UATTS" 
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PRINT 

PRINT 

PRINT 
CASE  Z 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 
CASE  3 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 
CASE  4 

PRINT 

PRINT 
END  SELECT 
! 

PRINT  USING  "K"i" 
! 

PRINT  "  " 


fl  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.1G  PROFILE- 
BUT  UITH  A  FACTOR  OF  TUO  TO  APPROXIMATE" 
LOOKING  UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398" 
FOR  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH" 
HEIGHT  IS  EXP<-.001«HEIGHT  ABOVE  65  METERS" 
AND  HEIGHT  TO  THE  -  1 . 46  BELOU  65  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
HEIGHT  TO  THE  -4/3  " 


THE  MINIMUM  DETECTABLE  SIGNAL  UAS  SET  AT" , Noise , "UATTS' 


! 

PRINT  "HIT  CONTINUE  TO  CONTINUE' 

PAUSE 

OUTPUT  KBD;"K"; 

! 

GRAPHICS  ON 

VIEUPORT  15, 1Z0, 10,7® 

! 

! 

IPLOT  RANGE  VERSUS  FREQUENCY 

Xmin=0 

Xnax= . 5 

Ymin=PROUND( Ynin-S0,Z) 

Ymax=PROUND( Ymax+50 , Z ) 

Dx=.05 

Dy= 1 00 

Xrange=Xmax-Xnin 

Yrange=Ymax- Ymin 

UINDOU  Xmin.Xmax.Ynin.Ymax 

AXES  Dx,Dy,Xmin,Ynin,Z,1 

CLIP  OFF 

! 

(LABEL  PLOT 
CSIZE  4, .6 
LDIR  0 
LORG  5 
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MOVE  .5*Xrange,Ymax*1 .05 

LABEL  "RANGE  VERSUS  TRANSDUCER  EFFICIENCY" 

(LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  5 

FOR  J=0  TO  Xnax  STEP  Dx 

CSIZE  4, .6 

MOVE  J,Ymin-.05*Yrange 

LABEL  J 
NEXT  J 

MOVE  .5*Xrange,Ynin-. 1*Yrange 
CSIZE  4,.B 
LABEL  "EFFICIENCY" 
(LABEL  VERTICAL  AXES 
LORG  8 
FOR  J=Ymin  TO  Ymax  STEP  Dy 

CSIZE  4, .6 

MOVE  Xmin-.00Z5*Xrange, J 

LABEL  USING  "K";J 
NEXT  J 
LDIR  PI/2 
LORG  B 

MOVE  Xnin- . 1 5*Xrange , . 5*Yrange+Ymin 
CSIZE  4,.G 

LABEL  "RANGE  (METERS)" 
CLIP  ON 
Ind=1 
FOR  Eff=.05  TO  .5  STEP  .05 

PLOT  Eff ,Ran(Ind) 

Ind=Ind-H 
NEXT  Eff 
I 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 
PAUSE 
GCLEAR 

OUTPUT  KBD;"K"; 
I 
I 

LINPUT  "UOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?" ,X$ 
Again=FNYes(X$> 
Mess_up=s1 
UHILE  Mess_up=1 
SELECT  Again 
CASE  1 

LINPUT  "UOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",XS 
Neu_va=FNYes<  X$ > 
Ano_change=1 
UHILE  Ano_change=1 
SELECT  Neu_va 
CASE  1 


o 


PRINT 


VARIABLE 


"1 
"Z 


TEMPERATURE 


CURRENT  VALUE" 

"  ,Temp,"CELSIUS" 


ATOMOSPHERIC  PRESSURE 
UATER  VAPOR  PRESSURE  ' 


, Atom_pres , "mb" 
Uat_pres , "mb" 


ANTENNA  DIAMETER 
PULSE  LENGTH 
POUER  TRANSMITTED 


"  ,  Ant_diam , "  m.  " 
"  ,Pul  se  ,  "  pis" 
,Pou_trans,"  UATTS" 


ATMOSPHERIC  PROFILE   ".Profile 


PRINT  USING  "K" 

PRINT  USING  "K" 

PRINT  USING  "K"|"3 

PRINT  "  " 

PRINT  USING  "K" ;"4 

PRINT  USING  "K";"S 

PRINT  USING  "K";MG 

PRINT  USING  "K";"7 

PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 

PRINT  "CHANGE" 

INPUT  Var 

SELECT  Var 

CASE  1 

INPUT  "TEMPERATURE=" , Temp 
CASE  Z 

INPUT  "ATMOSPHERIC  PRESSURE  IN  mb 
CASE  3 

INPUT  "UATER  VAPOR  PRESSURE  IN  mb=" , Uat_pres 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  m. =" , Ant_diam 
CASE  5 

INPUT  "PULSE  LENGTH  IN  ms=", Pulse 
CASE  G 

INPUT  "POUER  TRANSMITTED  IN  UATTS=" ,Pou_trans 


Atom_pres 


CASE  7 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
INPUT 


C-TZ  PROFILE  OF  ZA<-1.  IB)" 

FROM  UALTERS/KUNDEL  1981" 
SAME  AS  ONE  BUT  UITH  FACTOR" 
OF  TUO  TO  APPROXIMATE  LOOKING" 
UP  A  THERMAL  PLOOM" 

==>  C-TZ  PROFILE  OF  EXP(  -  .001  »Z  )  " 
FROM  UALTERS/KUNDEL  1981." 

==>  CT*Z  PROFILE  OF  ZA(-4/3)" 

Prof  i 1 e 


'ENTER  NUMBER  OF  DESIRED  PROFILE" 
IF  Profile=4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS-" 
END  IF 
CASE  ELSE 

PRINT  Var, "IS  NOT  ONE  OF  THE  OPTIONS- 
END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE( ENTER  Y  OR  N)?",X$ 
Neu_va=FNYes(X$) 
Mess_up=Z 
CASE  Z 

Mess_up=Z 
Ano_change=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N  Mill" 


Inver 


Ano_change=1 

LINPUT  "CHANGE  fi  VARIABLE  (ENTER  Y  OR  N)?",X$ 
New_va=FNYes(X$) 
END  SELECT 
END  UHILE 
CASE  2 

Mess_up=2 
CASE  ELSE 

PRINT  "YOUR  CHOICES  WHERE  Y  OR  N   MM!!" 
Mess_up=1 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again=FNYes(X$) 
END  SELECT 
END  UHILE 
END  UHILE 

PRINT  "THAT'S  ALL,  FOLKS" 
end 


CALCULTE  THE  ATTENUATION 
sub  Attenuat  ion( Atom_pres , At ten ,Freq , Tenp, Uat_pres  ) 
THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 

INPUT   ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
FREQUENCY  OF  SOUND  UAVE  IN  HERTZ 
TEMPERATURE  IN  DEGREES  CELCIUS 
UATER-VAPOR  PRESSURE  IN  MILLIBARS 

OUTPUT   ATTENUATION  IN  1 /METERS 


VARIABLES 
Aton_pres 
Atten 

Att  max 


IN  mb. 
CALCULATED  IN 


Fnax 

Freq 

H 

Pstar 

Temp 
Tstar 

Uat_pres 


INPUT  OF  ATMOSPHERIC  PRESSURE 
ATTENUATION  OF  ACOUSTIC  UAVE. 

SUBPROGRAM  ATTENUATION. 
VARIABLE  IN  SUBPROGRAM  ATTENUATION.   IT  IS  THE 

ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

ATTENUATION  FOR  THE  INPUT  CONDITIONS. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 
INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 
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!  OPERATOR. 

H=100*Uat_pres/fttom_pres 

Tstar«< 1 .8*Temp+49Z )/B19 

Pstar«Atom_pres/ 1014 

Fmax-(  10+G600*H+44400*H*H)*Ps-tar/TstarA.8 

A-tt_Piax=.0078*Fnax*TstarA(-Z.5)»EXP(7.77*<  1-1/Tstar)  ) 

F»Freq/Fi»iax 

A-tten=(Att_max/304.8)*<( . 1 8»F r Z+< Z*F»F/< 1 +F*F  )  ) "2  >  ' . 5 

Atten«< Atten+1 .74E- 10*Freq*Freq)/4. 3S 

SUBEND 

! 

def  FNYes(X$> 

ITHIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 

! 

!     INPUT     X$ 

OUTPUT    FNYes 
VARIABLES 

TeptpS      VARIABLE  STRING  USED  IN  FUNCTION  YES. 
X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 
RESPONSE  TO  YES  OR  NO  QUESTION. 
DIM  Temp$[ 1] 
Tenp$t  1  ,  U«TRIMS<X$> 
SELECT  Tenp$ 
CASE  "Y","y" 

RETURN  1 
CASE  "N" ,"n" 

RETURN  Z 
CASE  "  " 

RETURN  1 
CASE  ELSE 

RETURN  -Z 
END  SELECT 
FNEND 
l 
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APPENDIX    D 

COMPUTER   PROGRAM    4 
FULLER,    ROBERT  PR0G_4  ANTENNA   SIZE 

10  SEP    85 


****#*#PURPOSE******* 

THIS  PROGRAM  WILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  ANTENNA  SIZE  OF  THE  THE  ECHOSOUNDER. 
THE  FOLLOUING  INPUTS  ARE  REQUIRED: 
ATMOSPHERIC  DATA 

1 )ATOMOSPHERIC  PRESSURE  in  millibars 

Z)SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 
a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 
3)TEMPERATURE  IN  DEGREES  CELSIUS 
4)UATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNDER  DATA 

G)FREOUENCY 

G)POUER  TRANSMITTED  BY  ECHOSOUNDER  IN  UATTS 

7)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOUING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 
1 )RANGE<n. )  VERSUS  ANTENNA  DIAMETER  FOR  VARIOUS  FREQUENCIES 


*****  VARIABLES***** 


Again      TEST  VALUE  IN  WHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 
Ano_change  TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 

ANOUTHER  RUN. 
Ant_area   ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 

ANTENNA  DIAMETER. 
Ant_diam    INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 
Atom_pres   INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 
At-fcen      ATTENUATION  OF  ACOUSTIC  WAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
Bn         BANDWIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 

VERTICAL  VELOCITIES. 
C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

CI         VARIABLE  USED  TO  SCALE  INVERSION  HEIGHT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 
C3         VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 
CneZ<*>     ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 

PARAMETER.   CALCULATED  BASED  ON  SELECTION  OF 


PROFILE  FOR  TEMPERRTURE  STRUCTURE  PARAMETER 

AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 

PARAMETER. 
CteZ       VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 
CveZ       VALUE  OF  THE  VELOCITY  STRUCTURE   PARAMETER. 

VALUES  BASED  ON  CALCULATION  USING  ASSUMED  DISSAPATION 

RATE. 
Dx  STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 

Dy         STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 
Epsilon    DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 
Eff         TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 
Eff         EFFICIENCY  OF  ECHOSOUNDER  UHEN  RECEIVING  ACOUSTIC 

BACKSCATTER. 
Exc_att     EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 
F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
Fnax       FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 
Freq       INPUT  FREQUENCY  OF  ECHOSOUNDER. 
G  ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POWER  UNTIL 

POUER  UAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 
IND        INDEX  USED  FOR  LOOP  FOR  DIFFERENT  EFFICIENCIES 
Interval    MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 
Inver      HEIGHT  OF  INVERSION  LAYER. 
J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  UAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Mess_up    TEST  VALUE  TO  CHECK  RESPONSES  UHEN  OPERATOR  IS  INPUTTING 

RESPONSES. 
N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 
Neu_va     TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 

A  VARIABLE  BEFORE  A  NEU  RUN. 
Noise      ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 
Pou_back<*)  POUER  BACKSCATTERED  FROM  GIVEN  RANGE. 
Pou_ret<«>  POUER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 
Pou_trans   INPUT  OF  POUER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 
Profile    OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 
Pstar      VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 
Pulse      TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.   INPUT  BY 

OPERATOR. 
R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1         VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE 
Range* »>    ARRAY  OF  RANGE  VALUES. 
Remainder   REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 

NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 
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ATTENUATION. 
Rge        VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 
Rho        CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 
Sigma      FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 
Speed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 
Sumpou_back 

T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
Temp$      VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Tstar'      INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
Var        USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  UISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 
Uat_pres    ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 
Xmax       VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xmin       VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xrange     VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 
Ylabel$    LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Semi_log  FOR  PLOTTING. 
Ymax       LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Ymm       SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Yrange     RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 
Ze         INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back( 1500) ,Pou_ret< 1500) ,Range<  1500) 
DIM  CneZ( 1500), Ran< 50,6) 
! 

PLOTTER  IS  705,"HPGL" 

LINE  TYPE  1 

I 

! INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" , Temp 

INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" ,Atom_pres 

INPUT  "ENTER  UATER  VAPOR  PRESSURE  IN  MILLIBARS" ,Uat_pres 


INPUT  ECHOSOUNDER  DATA 
INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS=" , Pulse 
INPUT  "ENTER  POWER  TO  TRANSMITTER  IN  UATTS" ,Pou_trans 
G=.40        ! ANTENNA  EFFECTIVE  APERATURE  FACTOR 
Eff=.ZS 


! 

(SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 
PRINT  "YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES" 
PRINT  "PROFILE  TO  BE  USED." 
PRINT  "YOUR  SELECTIONS  ARE" 

PRINT  "      1       A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  " 
PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1381  PAGE  398  WHICH  " 

PRINT  "  GIVES  A  HEIGHT  TO  THE  -LIB  PROFILE" 

PRINT  "  " 

PRINT  "      Z      THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT- 
PRINT  "  UITH  A  FACTOR  OF  TUO  TO  APPROXIMATE  LOOKING" 
PRINT  "            UP  A  THERMAL  PLOOM" 
PRINT  "  " 

PRINT  "      3      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 

PRINT  "  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 

PRINT  "  IS  EXP<-.001*HEIGHT  ABOVE  G5  METERS  AND  HEIGHT" 

PRINT  "  TO  THE  -1.4E  BELOU  65  METERS" 

PRINT  "  " 

PRINT  "     4      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
PRINT  "  HEIGHT  TO  THE  -4/3  " 

I 

PRINT  "  " 
PRINT  "  " 
Mess_up= 1 
UHILE  Mess_up=1 

INPUT  "ENTER  THE  DESIRED  PROFILE  (1  OR  2  OR  3  OR  4)", Profile 
IF  Profile=4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" , Inver 
END  IF 
IF  Profiled  THEN 

Mess_up=0 
ELSE 

IF  Profile=Z  THEN 

Mess_up=0 
ELSE 

IF  Profile=3  THEN 

Mess_up=@ 
ELSE 

IF  Profile=4  THEN 

Mess_up=0 
ELSE 

PRINT  Profile,"  UAS  NOT  ONE  OF  THE  OPTIONS!!'" 
Mess_up=1 
END  IF 
END  IF 
END  IF 
END  IF 
END  UHILE 
OUTPUT  KBD;"K"; 
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I 
! 

Again= 1 
UHILE  Again=1 
FOR  J=1  TO  I 

Pou_back( J)=0 
Pou_ret< J)=0 
Range< J>=0 
CneZ( J>=0 
NEXT  J 
! 

(CONVERT  TEMPERATURE  TO  KELVIN 
T=Tercp+Z73 
! 

(CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE ( CELCIUS ) 
Speed_sound=Z0. 05* ( T )  * . 5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 
C=Z0.05*Z73'\5 
! 

1=1     ! INDEX  FOR  INNER  LOOP  (POUER  RETURNED  TO  NOISE) 
Ind=1   ! INDEX  FOR  SECOND  LOOP  (EFFICIENCY) 
Ifreq-1  ! INDEX  FOR  THIRD  LOOP  (FREQUENCY) 

Pou_back(0)=0     UNITIALIZE  VARIABLE  FOR  POUER  BACKSCATERRED 
Interval=( Speed_sound*Pul se* 1 . E-3)/Z 
! 

UNITIALIZE  PLOTTING  PARAMETERS  SO  CAN  COMPARE  DURING  RUN  TO  SCALE 
! PLOTS 
Yroin=1000 
Ymax=0 

(CALCULATE  RETURNED  SIGNAL  POUER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 
IFOR  EACH  FREQUENCY. 
Range<0>=0 

FOR  Freq=500  TO  1500  STEP  500 
Sunpou_back  =0 

FOR  Arvt_diam=.5  TO  B.Z5-Hreq  STEP  .Z5 
Ant_area=PI*(  Ant_diam/Z )"2 

Bn=Z*Freq*( 1  -  1 /< B/Speed_sound+ 1 ))!BANDUIDTH  FOR  3M/S  VERTICAL  VELOCITY 
Noise=1 .38E-Z3*Bn*(T)+Z.E-14!MINIMUM  DETECTABLE  SIGNAL  AS 
!      JOHNSON  NOISE(NEGLIGABLE)+  ESTIMATED  BACKGROUND 
! CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
CALL  Attenuat ion( Atom_pres , At  ten ,Freq, Temp, Uat_pres ) 
K=Z*PI*Freq/Speed_sound     ! UAVENUMBER 
1  =  1 

REPEAT 
Range* I >=Range< 1-1 )+Z 
SELECT  Profile 
CASE  1 

ITHIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
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ITULAROSA  BASIN,  NEU  MEXICO. 

!flN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 

! AVERAGING  TIME. 

CteZ-Z. 1Z*Range( I  r<-1 . 16) 
CASE  Z 

ITHIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  UEATHER  ABOVE  THE 
!TULAROSA  BASIN,  NEU  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1.8  UAS  INCLUDED  AS  PER  DIFFERENT 
! AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TUO  UAS 
! INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 

CteZ=Z*Z. 1Z»Range(I)A<-1. 1G) 
CASE  3 

!THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
(PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
! NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
IIS  EXP<-.001*HEIGHT  ABOVE  65  METERS  AND  HEIGHT" 
!TO  THE  -1.4G  BELOU  65  METERS" 
IF  Range(I)<65  THEN 

CteZ=75.5*Range(ir(-1  .46) 
ELSE 

CteZ=3.G6E-Z*EXP<-.001*Range< I ) ) 
END  IF 
CASE  4 

! CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!     THAT  IT'S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POUER 
!     THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  Cm.  /sec 
!THE  DEPENDENCE  UITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
!MY  LEFT  EAR 

! EQUATION  FROM  NEFF ,  1975 
SELECT  Ranged) /I nver 
CASE  <.9 

C3=((  .0Z4)»(TK(  .667)) 
CteZ=C3*(Range(I))"(-1 .33) 
CASE  .9  TO  1 
CteZ=CteZ 
C1=CteZ 
R1=Range(  I ) 
CASE  1  TO  1.3 

CteZ=10*((LGT(C3*Range< I > >-LGT< C3*R1 ))+LGT(C1 ) ) 
RZ=Range< I ) 
CZ  =  C-teZ 
CASE  ELSE 

CteZ  =  C3*(Range( I ) * ( -  1 . 33>-RZ"(  -  1 . 33 ) )+CZ 
END  SELECT 
END  SELECT 

! CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epsilon=(  .Z866/Range< I ) )*(  1  +  .07«( Range< I ) )  '  (  .  6 ) ) *  (  1  .  5  ) 
CveZ  =  Z»Epsilon'v(  .667) 
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(CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

(FORMULA  FROM  TETARSKI 

CneZ( I )  =  <  CteZ/( Z .98E+5) )  +  ( CveZ/( C*C ) ) 

(CALCULATE  THE  FRACTION  OF  POWER  BACKSCATTERED  FROM  INTERVAL. 

(EQUATION  FROM  NEFF ,  1975 

Sigma=< .0039*(KA<  1/3) )*CteZ )/< T) AZ 

I 

(CALCULATE  THE  EXCESS  ATTENUATION.   EQUATION  FROM  CLIFFORD,  1980 

(THE  EXCESS  ATTENUATION  IS  Exc_att 

IF  K50  THEN 

Remainder3® 
ELSE 

IF    K150  THEN 

Remainder=Range( I )  MODULO  10 
ELSE 

Remainder=Range< I )  MODULO  Z0 
END  IF 
END  IF 

IF  Remainder=0  THEN 
Rho=0 
L=0 
H=0 

Rge=Range<I)     (CONSTANT  IN  INTEGRAL 
R=0 
FOR  J=0  TO  Z*I 

F=CneZ(INT< J/2+1 )> 

F  =  F*(  1-R/(Rge)A( 1 .B7)  +  < R/(  Rge) ) A ( 1.G7) ) 

IF  J>0  THEN 

IF  J<Z«I  THEN 

IF  INT(J/Z)=J/Z  THEN 
L=L+F 
F=0 
ELSE 

H=H+F 
F=0 
END  IF 
END  IF 
END  IF 
Rho=Rho+F 
R=R+1 
NEXT  J 

Rho=Rho+4*L+Z*H 

Rho=(((Rho*.33)*K*K*1 .4B)"(-.B>  > 
N=(  Ant_diai*i/Rho)AZ 
IF  N<=1  THEN 

Ze=1/( 1+N) 
ELSE 

Ze=1 .5/( 1+N) 

(STEP  OF  1.5==>SEE  CLIFFORD  1980 
END  IF 
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Exc_att:=Ze*Ze 
ELSE 

Exc_att=Exc_att 
END   IF 
! 

•CALCULATE    THE   POUER   BACKSCATTERED 

Pou_back< I )=<Pou_trans*Ef f-Sumpou_back )*EXP( -Atten»Range( I ) ) 
Pou_back(  I  )=Pou_back(  I  )*Interval*Exc_att*Sigrna 
Surnpou_back=Sui»ipou_back+PoiJ_back<  I  ) 
•CALCULATE   THE   POUER   BACKSCATTERED   TO   THE    ANTENNA 
Pou_ret< I)-Pou_back(I >«EXP< -Atten*Range< I ) >»Ant_area»G*Ef f /Range( I )*Z 
PRINT    "RANGE  =  "  .Ranged  ) 
PRINT    "POUER   RETURNED-", Pou_ret( I ) 
PRINT    "FREQUENCY-", Freq 
1  =  1  +  1 
UNTIL   Pou_re+.U-1  )>0   AND  Pou_ret< I- 1 )<=Noise 
Ran( Ind,Ifreq)=Range(  I-Z) 
IF    Ran(Ind,Ifreq)>Ynax    THEN 

Ymax=Ran(  Ind, If req) 
END   IF 
IF    Randnd.IfreqKYnin    THEN 

Ymin=Ran(  Ind , If req) 
END    IF 
Ind=Ind+1 
Again= 1 
FOR   J-0   TO    I 

Pou_back<  J)=0 
Pou_ret(  J)=0 
Range( J)=0 
CneZ( J>=0 
NEXT   J 
NEXT    Ant_diam 
If req=If req+ 1 
Ind=1 
NEXT   Freq 
I 


OUTPUT   KBD;"K"; 

PRINT    "INPUT   CONDITIONS" 

I 
RINT    "    " 

"TEMPERATURE- 
"ATMOSPHERIC  PRESSURE- 
"WATER  VAPOR  PRESSURE- 
"PULSE  LENGTH- 
"TRANSDUCER  EFFICIENCY- 
"POUER  TRANSMITTED- 

PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED 

IF  Profile=4  THEN 

PRINT  "INVERSION  HEIGHT  ",Inver 


PRINT 

"  " 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

PRINT 

USING 

"K" 

,Temp,"CELCIUS" 

, Atom_pres , "mb" 

, Uat_pres , "mb" 

,Pul se , "ns" 

,Eff 

,Pou_trans,"  UATT5" 

" .Profile 
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END  IF 

SELECT  Profile 

CASE  1 

PRINT  " 
PRINT  " 
PRINT  " 

CASE  Z 

PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 

CASE  3 

PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 

CASE  4 

PRINT  " 
PRINT  " 
END  SELECT 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA"' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  " 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.1E  PROFILE" 

A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  " 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 
BUT  WITH  A  FACTOR  OF  TUO  TO  APPROXIMATE" 
LOOKING  UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA" 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398" 
FOR  NIGHT  CONDITIONS.   THE  DEPENDENCE  WITH" 
HEIGHT  IS  EXP(-.001*HEIGHT  ABOVE  B5  METERS" 
AND  HEIGHT  TO  THE  -  1 . 46  BELOU  65  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
HEIGHT  TO  THE  -4/3  " 


PRINT  USING  "K";"THE  MINIMUM  DETECTABLE  SIGNAL  WAS  SET  AT" , Noise , "UATTS" 


PRINT  "  " 

! 

I 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

OUTPUT  KBD;"K"; 

! 

GRAPHICS  ON 

VIEWPORT  15,120,10,70 


PLOT  RANGE  VERSUS  ANTENNA  DIAMETER  FOR  VARIOUS  FREQUENCIES 
Xroin=0 
Xmax=6 

Ynin=PROUND( Ynin-50,Z) 
Ynax=PROUND( Ymax+50 , Z ) 
Dx=1 
Dy= 1 00 

Xrange=Xmax-Xnin 
Yrange=Ymax-Ymin 
WINDOW  Xi*iin,Xmax,Ynin,Ynax 
AXES  Dx,Dy,Xmin,Ymin,1 ,1 
CLIP  OFF 
! 

"LABEL  PLOT 
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CSIZE  4..G 
LDIR  0 
LORG  5 

MOVE  .5*Xrange,Ymax*1.05 
LABEL  "RANGE  VERSUS  ANTENNA  DIAMETER" 
(LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  5 

FOR  J=0  TO  Xmax  STEP  Dx 
CSIZE  4, .6 

MOVE  J,Ynin-.05*Yrange 
LABEL  J 
NEXT  J 

MOVE  . 5*Xrange,Ymin-. 1*Yrange 
CSIZE  4, .6 

LABEL  "ANTENNA  DIAMETER  <m.>" 
! LABEL  VERTICAL  AXES 
LORG  8 

FOR  J=Ynin  TO  Ymax  STEP  Dy 
CSIZE  4, .6 

MOVE  Xnin-.00Z5*Xrange,J 
LABEL  USING  "K";J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  Xi»iin-  .  1  5*Xrange ,  .  5*  Yrange+Yrnin 
CSIZE  4, .6 

LABEL  "RANGE  (METERS)" 
CLIP  ON 
Ind=1 
If req=1 
FOR  Freq=500  TO  1500  STEP  500 

FOR  Ant_diam=.5  TO  G.25-Ifreq  STEP  .25 
IF  INT(Ifreq/Z)=Ifreq/2  THEN 

LINE  TYPE  3 
ELSE 

LINE  TYPE  1 
END  IF 

PLOT  An-t_diam,Ran<  Ind.Ifreq) 
Ind=Ind+1 
NEXT  Ant_dian 
LDIR  0 
LORG  Z 
CSIZE  3..B 

MOVE  Ant_diam-. 1 ,Ran( Ind-1 .Ifreq) 
LINE  TYPE  1 

LABEL  USING  "K" ; Freq, "Hz" 
If req=If req+ 1 
Ind=1 
NEXT  Freq 
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! 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBD;"K"; 


I 

LINPUT  "WOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again=FNYes(X$) 
Mess_up=1 
WHILE  Mess_up=1 
SELECT  Again 
CASE  1 

LINPUT  "WOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?" 
Neu_va=FNYes(X$) 
Ano_change= 1 
WHILE  Ano_change=l 
SELECT  Neu_va 
CASE  1 
PRINT  "       VARIABLE  CURRENT  VALUE" 


X$ 


TEMPERATURE        ", Temp , "CELSIUS" 
ATOMOSPHERIC  PRESSURE  " , Aton_pres , "mb" 
WATER  VAPOR  PRESSURE 


ANTENNA  DIAMETER 

PULSE  LENGTH 

POWER  TRANSMITTED  " ,Pou_trans , "  WATTS' 


, Wat_pres , "mb" 

,  Ant_diam  ,  "  it. ' 
.Pulse,"  ns" 


ATMOSPHERIC  PROFILE 


Profile 


PRINT  USING  "K'V'1 

PRINT  USING  "K'V'Z 

PRINT  USING  "K";"3 

PRINT  "  " 

PRINT  USING  "K";"4 

PRINT  USING  "K";"5 

PRINT  USING  "K";"G 

PRINT  USING  "K";"7 

PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 

PRINT  "CHANGE" 

INPUT  Var 

SELECT  Var 

CASE  1 

INPUT  "TEMPERATURE  =  ",TePip 
CASE  Z 

INPUT  "ATMOSPHERIC  PRESSURE  IN  nb=" , Atom.pres 
CASE  3 

INPUT  "WATER  VAPOR  PRESSURE  IN  mb=" ,Uat_pres 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  m. =" , Ant_dian 
CASE  S 

INPUT  "PULSE  LENGTH  IN  ns=", Pulse 
CASE  G 

INPUT  "POWER  TRANSMITTED  IN  WATTS="  Pou  trans 


CASE  7 
PRINT 
PRINT 
PRINT 
PRINT 


1   ==>  C-TZ  PROFILE  OF  ZA(-1.1G)" 
FROM  WALTERS/KUNDEL  1981" 

Z   ==>  SAME  AS  ONE  BUT  WITH  FACTOR" 

OF  TWO  TO  APPROXIMATE  LOOKING" 
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PRINT  "         UP  A  THERMAL  PLOOM" 

PRINT  "  3   =•=>  CfZ  PROFILE  OF  EXP(  -  .  001  «Z  )  " 

PRINT  "        FROM  UALTERS/KUNDEL  1981." 

PRINT  "  4   «>  CT"Z  PROFILE  OF  Z*<-4/3>" 

INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" , Prof  lie 

IF  Profile=4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS-" , I nver 
END  IF 
CASE  ELSE 

PRINT  Var,"IS  NOT  ONE  OF  THE  OPTIONS" 
END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE( ENTER  Y  OR  N)?",X$ 
Neu_va=FNYes(X$) 
Mess_up=Z 
CASE  Z 

Mess_up,=  Z 
Ano_change=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  Ml!  I" 
Ano_change*=  1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X$ 
Neu_va=FNYes(X$) 
END  SELECT 
END  WHILE 
CASE  Z 

Mess_up=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N   MM!!" 
Mess_up= 1 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again=FNYes(X$) 
END  SELECT 
END  WHILE 
END  WHILE 

PRINT  "THAT'S  ALL,  FOLKS" 
end 


CALCULTE  THE  ATTENUATION 
sub  Attenuat  ion(  Atoi*i_pres  ,  At  ten  ,  Freq  ,  Tenp  ,  Uat_pres  ) 
THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 

INPUT    ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
FREQUENCY  OF  SOUND  WAVE  IN  HERTZ 
TEMPERATURE  IN  DEGREES  CELCIUS 
UATER-VAPOR  PRESSURE  IN  MILLIBARS 

OUTPUT   ATTENUATION  IN  1 /METERS 
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! VARIABLES 


I 


Atoi*i_pres 
Atten 

Att  max 


Fnax 

Freq 

H 

Pstar 

Temp 
Tstar 


INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 
ATTENUATION  OF  ACOUSTIC  WAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
VARIABLE  IN  SUBPROGRAM  ATTENUATION.   IT  IS  THE 

ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 

ATTENUATION  FOR  THE  INPUT  CONDITIONS. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 

VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 
VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 
INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 


IN  SUBPROGRAM  ATTENUATION. 
Uat_pres    ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
H=  l00*Ua-t_pres/Aton_pres 
Tstar=( 1 .8»Tenp+43Z)/519 
Pstar=AtoM_pres/ 1014 

Fnax=(  10+6G00*H-f44400»H»H)*Pstar/Tstar"  .8 
A-t-t_max  =  .0078*Fmax»TstarA(-Z.5)*EXP(7.77»<  1-1 /Tstar)  ) 
F=Freq/Fmax 

Atten=( Att_Rax/304.8)*(  (  . 1 8»F ) *Z  +  ( Z*F»F/<  l+F^F ) >*Z )* .5 
A-tten=(  Atten+I  .74E-  1  0* Freq* Freq  )/4.  35 
SUBEND 


def  FNYes(X$) 
THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 


INPUT 


OUTPUT 
VARIABLES 
Temp$ 
X$ 


X$ 


FNYes 


VARIABLE  STRING  USED  IN  FUNCTION  YES. 
STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 
RESPONSE  TO  YES  OR  NO  QUESTION. 
DIM  Temp$[ 1 ] 
Temp$[ 1 , 1 ]=TRIMI(X$) 
SELECT  Temp$ 
CASE  "Y'V'y" 
RETURN  I 


CASE 


CASE 


CASE 


"N"  ,"n' 
RETURN 

RETURN 
ELSE 
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RETURN  -2 
END  SELECT 
FNEND 
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APPENDIX    E 

COMPUTER   PROGRAM    5 

(FULLER,    ROBERT  PR0G_5  POUER 

! 

! 10  SEP   85 

! 

I 

! #*#»#*#PURPOSE******* 

! 

!      THIS  PROGRAM  UILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 

(ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  POUER  TRANSMITTED  BY  THE  ECHOSOUNDER 

IFOR  A  RANGE  OF  FREQUENCIES 

!THE  FOLLOUING  INPUTS  ARE  REQUIRED: 

!     ATMOSPHERIC  DATA 

!     1  >ATOMOSPHERIC  PRESSURE  in  millibars 

Z)SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 
a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 

3)TEMPERATURE  IN  DEGREES  CELSIUS 

4)WATER  VAPOR  PRESSURE  IN  millibars 

ECHOSOUNDER  DATA 

5)FREQUENCY 

G) ANTENNA  DIAMETER  IN  METERS 

7)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

THE  PROGRAM  OUTPUTS  THE  FOLLOUING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 

1)RANGE<m.)  VERSUS  POUER  TRANSMITTED  IN  UATTS  FOR  VARIOUS  FREQUENCIES 


****♦ VARIABLES***** 


Again      TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

UISHES  TO  MAKE  ANOUTHER  RUN. 
Ano_change  TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

UISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 

ANOUTHER  RUN. 
Ant_area   ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 

ANTENNA  DIAMETER. 
Ant_diam    INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 
Atom_pres   INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 
Atten      ATTENUATION  OF  ACOUSTIC  WAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
Bn         BANDUIDTH  OF  RECEIVER  NEEDED  TO  RECEIVE  SIGNALS 

DOPPLER  SHIFFED  BY  THREE  METER  PER  SECOND 

VERTICAL  VELOCITIES. 
C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

CI         VARIABLE  USED  TO  SCALE  INVERSION  HEIGHT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 
C3         VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 
CneZ<*)     ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 
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!  PARAMETER.   CALCULATED  BASED  ON  SELECTION  OF 

!  PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 

!  AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 

!  PARAMETER. 

!   C-teZ        VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 
CveZ       VALUE  OF  THE  VELOCITY  STRUCTURE   PARAMETER. 

VALUES  BASED  ON  CALCULATION  USING  ASSUMED  DISSAPATION 

RATE. 
Dx         STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 
Dy         STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 
Epsilon    DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 
Eff        TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 
Eff        EFFICIENCY  OF  ECHOSOUNDER  WHEN  RECEIVING  ACOUSTIC 

BACKSCATTER. 
Exc_att     EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 
F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
Fmax       FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 
Freq       INPUT  FREQUENCY  OF  ECHOSOUNDER. 
G  ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POUIER  UNTIL 

POWER  UAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 
Ifreq      INDEX  FOR  INCREMENTS  OF  FREQUENCY 

IND        INDEX  USED  FOR  LOOP  FOR  DIFFERENT  TRANSMITTED  POUER 
Interval    MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 
Inver      HEIGHT  OF  INVERSION  LAYER. 
J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  UAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Mess_up    TEST  VALUE  TO  CHECK  RESPONSES  UHEN  OPERATOR  IS  INPUTTING 

RESPONSES. 
N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 
Neu_va     TEST  VALUE  IN  CHECKING  IF  OPERATOR  UISHES  TO  CHANGE 

A  VARIABLE  BEFORE  A  NEU  RUN. 
Noise      ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 
Pou_back<*>  POUER  BACKSCATTERED  FROM  GIVEN  RANGE. 
Pou_ret(»)  POUER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 
Pou_trans   INPUT  OF  POUER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 
Profile    OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 
Pstar      VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 
Pulse      TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.   INPUT  BY 

OPERATOR. 
R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1         VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE 
RangeU)    ARRAY  OF  RANGE  VALUES. 
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Remainder   REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 

NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 
ATTENUATION. 

Rge        VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 

Rho        CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 
ATTENUATION. 

Sigma      FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 

Speed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 

Sumpou_back 

T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 

Temp$      VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Ts-tar      INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
Var        USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  WISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 
Uat_pres    ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 
Xmax       VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xmin       VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xrange     VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 
YlabelS    LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Semi_log  FOR  PLOTTING. 
Ymax       LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Ymin       SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Yrange     RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 
Ze         INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back( 1500) ,Pou_ret<  1500) ,Range(  1500) 
DIM  CneZ(  1500), Ran( 50,6) 
! 

PLOTTER  IS  705,MHPGL" 
LINE  TYPE  1 

! 

•INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" , Temp 

INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS" ,Atom_pres 

INPUT  "ENTER  UATER  VAPOR  PRESSURE  IN  MILLIBARS" ,Uat_pres 


INPUT  ECHOSOUNDER  DATA 
INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS=" , Pulse 
INPUT  "ENTER  THE  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS" , Ant_diam 
G=.40        ! ANTENNA  EFFECTIVE  APERATURE  FACTOR 


3o 


Eff=.Z5 


(SELECT  THE  PROFILE  FOR  THE 


YOU  HAVE  A  CHOICE  OF 
PROFILE  TO  BE  USED." 
YOUR  SELECTIONS  ARE" 


1 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

I 

PRINT 

PRINT 

Mess_up=1 

UHILE    Mess_up=1 

INPUT  "ENTER 
IF  Profile=4 


TEMPERATURE 
TEMPERATURE 


STRUCTURE 
STRUCTURE 


PROFILES" 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  ' 
PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  UHICH  " 
GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 
UITH  A  FACTOR  OF  TUO  TO  APPROXIMATE  LOOKING" 
UP  A  THERMAL  PLOOM" 

A  TEMPERATURC  STRUCTURE  PROFILE  BASED  ON  DATA  AS- 
PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
IS  EXP(-.001*HEIGHT  ABOVE  G5  METERS  AND  HEIGHT" 
TO  THE  -1.4G  BELOU  65  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
HEIGHT  TO  THE  -4/3  " 


THE  DESIRED  PROFILE 
THEN 


( 1  OR  2  OR  3  OR  4)" .Profile 


INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS" , Inver 
END  IF 
IF  Profiled  THEN 

Mess_up=0 
ELSE 

IF  Profile=Z  THEN 

Mess_up=0 
ELSE 

IF  Profile=3  THEN 

Mess_up=0 
ELSE 

IF 


END 
END  IF 


END 
IF 


Profile=4  THEN 
Mess_up=0 
ELSE 

PRINT  Profile 
Mess_up= 1 
END  IF 
IF 


UAS  NOT  ONE  OF  THE  OPTIONS!!!! 


o 


END  UHILE 
OUTPUT  KBD;"K" ; 


Again=1 
UHILE  Again- 1 
FOR  J-1  TO  I 

Poui_back(  J>=0 

Pou_ret<J)=0 

Range< J)=0 

CneZ< J)=0 
NEXT  J 
! 

(CONVERT  TEMPERATURE  TO  KELVIN 
T=Temp+273 
! 

! CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE* CELCIUS ) 
Speed_sound=Z0. 05* ( T ) " . 5 

(CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 
C=Z0.05*Z73A.5 
! 

1  =  1     I  INDEX  FOR  INNER  LOOP  ( POUER  RETURNED  TO  NOISE) 
Ind-1   'INDEX  FOR  SECOND  LOOP  (EFFICIENCY) 
Ifreq=l  'INDEX  FOR  THIRD  LOOP  (FREQUENCY) 

Pou_back<0>=0     UNITIALIZE  VARIABLE  FOR  POWER  BACKSCATERREO 
Interval=( Speed_sound* Pulse* 1 . E-3)/Z 
Ant_area=PI*<  Ant_dian/Z )"Z 
! 

UNITIALIZE  PLOTTING  PARAMETERS  SO  CAN  COMPARE  DURING  RUN  TO  SCALE 
(PLOTS 
Ymin=1000 
Ynax=0 

'CALCULATE  RETURNED  SIGNAL  POUER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 
IFOR  EACH  FREQUENCY. 
Range* 0)=0 

FOR  Freq=500  TO  1500  STEP  500 
FOR  Pou_trans=50  TO  500  STEP  50 
Bn=Z*Freq*< 1  -  1 /<  67Speed_sound+ 1 )) (BANDWIDTH  FOR  3M/S  VERTICAL  VELOCITY 
Noise=1.38E-Z3*Bn*(T)+Z.E-14!MINIMUM  DETECTABLE  SIGNAL  AS 
!      JOHNSON  NOISE(NEGLIGABLE)*  ESTIMATED  BACKGROUND 
(CALCULATE  THE  ATTENUATION  COEFFICIENT  BASED  ON  ATMOSPHERIC  DATA 
CALL  Attenuat  ion(  Atoro_pres  ,  At  ten  ,Freq  ,  Tenp  ,Uat_pres  ) 
K=Z*PI*Freq/Speed_sound     ! UAVENUMBER 
1  =  1 

Sumpou_bacl<  =0 
REPEAT 
Range( I )  =  Range(  I- 1 )+Z 
SELECT  Profile 
CASE  1 
'THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
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IAS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 

I  UAS  FOR  MIDDAY  AND  CLEAR  UEATHER  ABOVE  THE 

ITULAROSA  BASIN,  NEU  MEXICO. 

IAN  ADDITONAL  FACTOR  OF  1 . 8  WAS  INCLUDED  AS  PER  DIFFERENT 

! AVERAGING  TIME. 

CteZ=Z. 1Z«Range< I  )A<-1 . 16) 
CASE  Z 

ITHIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
I  UAS  FOR  MIDDAY  AND  CLEAR  UEATHER  ABOVE  THE 
ITULAROSA  BASIN,  NEU  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 
'AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TUO  UAS 
(INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 

CteZ=Z»Z. 1Z*Range( I )"<-!. IB) 
CASE  3 

ITHIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
(PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
I  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
IIS  EXP(-.001*HEIGHT  ABOVE  65  METERS  AND  HEIGHT" 
I  TO  THE  -1.4G  BELOU  G5  METERS" 
IF  Range(  I  XB5  THEN 

C-teZ=75.5*Range(  I  )*<  -1  .46) 
ELSE 

CteZ=3.66E-Z»EXP(-.001*Range( I ) ) 
END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
I     THAT  IT'S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POUER 
I     THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  Cm. /sec 
ITHE  DEPENDENCE  UITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
I  MY  LEFT  EAR 

I  EQUATION  FROM  NEFF ,  1975 
SELECT  Range( I)/Inver 
CASE  <.9 

C3=(( .©Z4)»(T>*(.BB7)) 
CteZ=C3*( Range( I ) >*< -  1 . 33) 
CASE  .9  TO  1 
CteZ=CteZ 
C1=CteZ 
R1=Range( I ) 
CASE  1  TO  1.3 

CteZ=10"((LGTCC3*Range(I ))-LGT(C3»R1 )>+LGT(C1  )  ) 
RZ=Range( I ) 
CZ=CteZ 
CASE  ELSE 

CteZ=C3»(Range( I ) " ( -  1 . 33 )-RZ"< -  1 . 33 ) )  +  CZ 
END  SELECT 
END  SELECT 
"CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 


33 


Epsilon=(  .Z866/Range< I ) )*(  1  +.07*< Range(  I ) )A( .6) r< 1.5) 

CveZ=Z*Epsilon*<  .GB7) 

! CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 

! FORMULA  FROM  TETARSKI 

CneZ(  I )  =  (  CteZ/( Z . 98E+5) )  +  ( CveZ/( C*C  ) ) 

•CALCULATE  THE  FRACTION  OF  POWER  BACKSCATTERED  FROM  INTERVAL. 

! EQUATION  FROM  NEFF .  1975 

Sigma=< .0039«<KA< 1/3) >«CteZ)/< T)*Z 

! 

(CALCULATE  THE  EXCESS  ATTENUATION.   EQUATION  FROM  CLIFFORD,  1980 
!THE  EXCESS  ATTENUATION  IS  Exc_att 
IF  K50  THEN 

Remainder =0 
ELSE 

IF  K150  THEN 

Remainder=Range(  I )  MODULO  10 
ELSE 

Remainder=Range( I )  MODULO  Z0 
END  IF 
END  IF 

IF  Remainder=0  THEN 
Rho=0 
L-0 
H=0 

Rge=Range(I)     (CONSTANT  IN  INTEGRAL 
R=0 
FOR  J=0  TO  Z*I 

F=CneZ(INT( J/Z  +  1  )) 

F=F*(  1-R/(Rge)"(  1  .  57)  +  (  R/(  Rge  )  >*<  1  .  67  )  ) 

IF  J>0  THEN 

IF  J<Z»I  THEN 

IF  INT(J/Z)=J/Z  THEN 
L=L+F 
F=0 
ELSE 

H=H+F 
F=0 
END  IF 
END  IF 
END  IF 
Rho=Rho+F 
R=R+1 
NEXT  J 

Rho=Rho+4*L4Z*H 

Rho=< ( (Rho*.33)»K*K*1 .4B)*(-.B) ) 
N=(An-t_diam/Rho)AZ 
IF  N<=1  THEN 

Ze=1/<  1+N) 
ELSE 

Ze=1 .5/( 1+N) 
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!STEP  OF  1.5-=>SEE  CLIFFORD  198® 
END  IF 

Exc_att=Ze*Ze 
ELSE 

Exc_att=Exc_att 
END  IF 
! 

! CALCULATE  THE  POWER  BACKSCATTERED 

Pou_back( I )=<Pou_trans*Ef f -Sunpou_back )*EXP(  -Atten»Range( I > ) 
Pou_back(  I  )=»Pou_back(  I  )*Interval*Exc_att«Sigma 
Sunpou_back  =  Sunpou_back4Pou_back(  I ) 
(CALCULATE  THE  POUER  BACKSCATTERED  TO  THE  ANTENNA 
Pou_ret<  I  )-Pou_back<  I  )*EXP(  -At ten* Range (  I  )  )*An-t_area»G*Ef  f  /Range<  I  )  *  2 
PRINT  "RANGE  =  "  .Ranged) 
PRINT  "POUER  RETURNED-" ,Pou_ret(  I ) 
PRINT  "FREQUENCY=",Freq 
1  =  1  +  1 
UNTIL  Pou_ret<  1-1  )>0  AND  Pou»_ret(  I-  1  )<=Noise 
Ran( Ind, If req)=Range<  1-2 ) 
IF  Ran<Ind,Ifreq)>Ymax  THEN 

Ynax=Ran(  Ind ,  If req) 
END  IF 
IF  Randnd.IfreqKYmin  THEN 

Ymin=Ran(  Ind, If req) 
END  IF 
Ind=Ind4l 
Again=1 
FOR  J-0  TO  I 

Pou_back(  J)=0 
Pou_ret(  J>=0 
Range( J)=0 
Cne2( J)=0 
NEXT  J 
NEXT  Pou_trans 
If req=If req+1 
Ind=1 
NEXT  Freq 
I 
I 

OUTPUT  KBD;"K"; 
PRINT  "INPUT  CONDITIONS" 


PRINT  "  " 

PRINT  USING  "K" ; "TEMPERATURE = 
PRINT  USING  "K"; "ATMOSPHERIC  PRESSURE- 
PRINT  USING  "K"; "WATER  VAPOR  PRESSURE- 
PRINT  USING  "K"; "PULSE  LENGTH- 
PRINT  USING  "K"; "TRANSDUCER  EFFICIENCY- 
PRINT  U5ING  "K"; "POWER  TRANSMITTED- 
PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED 


,Temp,"CELCIUS" 

, Atom_pres , "mb" 

,Uat_pres , "nb" 

,Pul se , "ms" 

,Eff 

,Pou_trans,"  UATTS" 

"  .Profile 
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IF  Profile=4  THEN 

PRINT  "INVERSION  HEI 
END  IF 

SELECT  Profile 
CASE  1 

PRINT  "     1 
PRINT 
PRINT 
CASE  Z 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
CASE  3 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
CASE  4 

PRINT 
PRINT 
END  SELECT 


GHT  ".Inver 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 

A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  ' 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 
BUT  UITH  A  FACTOR  OF  TWO  TO  APPROXIMATE" 
LOOKING  UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA' 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398" 
FOR  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH" 
HEIGHT  IS  EXP(-.001*HEIGHT  ABOVE  G5  METERS" 
AND  HEIGHT  TO  THE  -1.46  BELOU  G5  METERS" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
HEIGHT  TO  THE  -4/3  " 


PRINT  USING  "K'VTHE  MINIMUM  DETECTABLE  SIGNAL  UAS  SET  AT"  .Noise ,  "UATTS' 
! 

PRINT  "  " 


I 

PRINT    "HIT    CONTINUE    TO    CONTINUE" 

PAUSE 

OUTPUT    KBD;"K"; 

I 

GRAPHICS   ON 

VIEUPORT    15,120,10,70 


PLOT    RANGE    VERSUS   POUER   TRANSMITTED   FOR   VARIOUS    FREQUENCIES 
Xnin=0 
Xrnax=600 

Ymin=PROUND<  Ymin-50,Z ) 
Ynax=PROUND( Ymax+50, Z ) 
Dx=100 
Dy= 1 00 

Xrange=Xi»iax-Xmin 
Yrange=Ymax-Ymin 
UINDOU   Xmin,Xnax,Yi«iin,Ynax 
AXES    Dx.Dy.Xmin.Ymin.l  ,1 
CLIP    OFF 
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! LABEL  PLOT 
CSIZE  4..G 
LDIR  0 
LORG  5 

MOVE  .5*Xrange,Ynax*1 .05 
LABEL  "RANGE  VERSUS  POWER  TRANSMITTED" 
I  LABEL  HORIZONTAL  AXES 
LDIR  0 
LORG  5 

FOR  J-0  TO  Xnax  STEP  Dx 
CSIZE  4, .6 

MOVE  J,Ynin-.05*Yrange 
LABEL  J 
NEXT  J 

MOVE  . 5*Xrange,Ynin-. 1*Yrange 
CSIZE  4, .6 

LABEL  "POWER  TRANSMITTED  (WATTS)" 
! LABEL  VERTICAL  AXES 
LORG  8 

FOR  J=Ymin  TO  Ytnax  STEP  Dy 
CSIZE  4, .6 

MOVE  Xmin-.00Z5*Xrange,J 
LABEL  USING  "K";J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  Xmn-.  15*Xrange,  .  5*Yrange+Ymin 
CSIZE  4, .6 

LABEL  "RANGE  (METERS)" 
CLIP  ON 
Ind=1 
If req=1 

FOR  Freq=500  TO  1500  STEP  500 
FOR  Pou_trans=50  TO  500  STEP  50 
IF  INT(Ifreq/Z)=Ifreq/Z  THEN 

LINE  TYPE  3 
ELSE 

LINE  TYPE  1 
END  IF 

PLOT  Pou_trans,Ran< Ind.Ifreq) 
Ind=Ind+1 
NEXT  Pou_-trans 
LDIR  0 
LORG  Z 
CSIZE  3, .6 

MOVE  Poui_trans-50,Ran(  Ind-  1  ,  Ifreq) 
LINE  TYPE  1 

LABEL  USING  "K" ; Freq , "Hz" 
I f req=If req+ 1 
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Ind=1 

NEXT  Freq 
i 

PRINT  "HIT  CONTINUE 

TO  CONTINUE" 

PAUSE 

GCLEAR 

OUTPUT  KBD^'K"; 
! 

! 

LINPUT  "UOULD  YOU  LIKE  TO 

MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 

Again=FNYes(X$> 

Mess_up=1 

■ 

WHILE  Mess_up=1 

SELECT  Again 

CASE  1 

LINPUT  "UOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X$ 

Neu_va=FNYe£ 

,(X$) 

Ano_change=1 

WHILE  Ano_change=1 

SELECT  Neu 

j_va 

CASE  1 

PRINT 

" 

VARIABLE           CURRENT  VALUE" 

PRINT 

USING 

"K";"1   TEMPERATURE        " ,Tenp, "CELSIUS" 

PRINT 

USING 

"K";"Z   ATOMOSPHERIC  PRESSURE  " , Atom_pres , "nb" 

PRINT 

USING 

"K";"3   WATER  VAPOR  PRESSURE  " , Wat_pres , "nb" 

PRINT 

■■  ■■ 

PRINT 

USING 

"K";"4   ANTENNA  DIAMETER     ",Ant_diam,"  m." 

PRINT 

USING 

"K";"5  PULSE  LENGTH         ".Pulse,"  ms" 

PRINT 

USING 

"K";"B   POWER  TRANSMITTED  " ,Pou_trans , "  WATTS" 

PRINT 

USING 

"K";"7   ATMOSPHERIC  PROFILE   ".Profile 

PRINT 

•i  •• 

PRINT 

"ENTER 

THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 

PRINT 

"CHANGE" 

INPUT 

Var 

SELECT  Var 

CASE    1 

INPUT    "TEMPERATURE=",Tenp 
CASE    2 


INPUT 

"ATMOSPHERIC  PRESSURE  IN  nb=" , Atom_pres 

CASE  3 

INPUT 

"WATER  VAPOR  PRESSURE  IN  mb=" ,Wat_pres 

CASE  4 

INPUT 

"ANTENNA  DIAMETER  IN  n. =" , Ant_diam 

CASE  S 

INPUT 

"PULSE  LENGTH  IN  ms=" , Pulse 

CASE  G 

INPUT 

"POWER  TRANSMITTED  IN  WATTS=" ,Pou_trans 

CASE  7 

PRINT  " 

1   ==>  CtAZ  PROFILE  OF  ZA(-1.16)" 

PRINT  " 

FROM  WALTERS/KUNDEL  1981" 
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PRINT  "  Z   -->  SAME  AS  ONE  BUT  UITH  FACTOR- 
PRINT  "        OF  TUO  TO  APPROXIMATE  LOOKING" 
PRINT  "        UP  A  THERMAL  PLOOM" 
PRINT  "  3   -=>  Ct'Z  PROFILE  OF  EXP( - . 001 «Z ) " 
PRINT  "        FROM  UALTERS/KUNDEL  1981." 
PRINT  "  4   — >  CT"Z  PROFILE  OF  Z*(-4/3)" 
INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" , Prof  lie 
IF  Profile-4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS-" , Inver 
END  IF 
CASE  ELSE 

PRINT  Var,"IS  NOT  ONE  OF  THE  OPTIONS" 
END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE( ENTER  Y  OR  N)?",XS 
Neu_va=FNYes(X$) 
Mess_up=Z 
CASE  Z 

Mess_up=Z 
Ano_change=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  !  M  ! ! " 
Ano_change= 1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N)?",X3! 
Neu_va=FNYes(X$) 
END  SELECT 
END  WHILE 
CASE  Z 

Mess_up=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N   I ! I ! ! ! " 
Mess_up=1 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?",X$ 
Again=FNYes(X$) 
END  SELECT 
END  UHILE 
END  UHILE 

PRINT  "THAT'S  ALL,  FOLKS" 
end 
! 
! 

!CALCULTE  THE  ATTENUATION 

sub  Attenuat ion( Atom_pres , At ten ,Freq, Temp,Uat_pres  ) 
ITHIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
!IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 
I 

!  INPUT   ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
!  FREQUENCY  OF  SOUND  UAVE  IN  HERTZ 

!  TEMPERATURE  IN  DEGREES  CELCIUS 

!  UATER-VAPOR  PRESSURE  IN  MILLIBARS 

I 
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!    OUTPUT   ATTENUATION  IN  1 /METERS 


VARIABLES 
Atom_pres 
Atten 

Att  max 


INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

ATTENUATION  OF  ACOUSTIC  UAVE.   CALCULATED  IN 
SUBPROGRAM  ATTENUATION. 

VARIABLE  IN  SUBPROGRAM  ATTENUATION.   IT  IS  THE 
ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 
ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 


Fmax 

Freq 

H 

Pstar 


RATIO  OF  FREQUENCY  TO  FREQUENCY 
FREQUENCY  OF  MAXIMUM  ATTENUATION 

ATTENUATION. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 
VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS 


AT  MAXIMUM  ATTENUATION. 
USED  IN  SUBPROGRAM 


ATTENUATION. 
FOR  INTERMEDIATE 


VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION 
CALCULATION. 
Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
Tstar      INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
Uat_pres    ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 
OPERATOR. 
H=  100»Uat_pres/Atom_pres 
Tstar=< 1 .8*Temp+49Z)/5l9 
Pstar=Atom_pres/10!4 

Fmax=( 10+6G00«H+44400*H»H)*Pstar/TstarA . 8 
Att_max=.0078*Fmax»Tstar"<-Z.5)»EXP(7.77»<  1- 1 /Tstar  )  ) 
F=Freq/Fmax 

Atten=(  At  t_max/  304.8  )•(  (  .  1  8»F  )  "2  +  (  Z«F*F/(  1+F»F )  )*Z ■>*  .5 
Atten=( Atten+1 .74E- 1 0»Freq»Freq  )  /4 . 35 
SUBEND 
i 

def  FNYes(XS) 
THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 


INPUT 


OUTPUT 
VARIABLES 
Tempi 
X$ 


X$ 


FNYes 


VARIABLE  STRING  USED  IN  FUNCTION  YES. 

STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 
RESPONSE  TO  YES  OR  NO  QUESTION. 
DIM  Temp$[ 1 ] 
Temp$[ 1 , 1 ]=TRIM$( X$> 
SELECT  TempS 


CASE 


CASE 


CASE 


••Y",'V 
RETURN 
"N"  ,"n' 
RETURN 
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RETURN  1 
CASE  ELSE 

RETURN  -Z 
END  SELECT 
FNEND 
! 
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APPENDIX    F 
COMPUTER   PROGRAM    6 


FULLER,    ROBERT  PR0G_G  NOISE 

10  SEP   85 

*#»#*»#PURPOSE* ****** 

THIS  PROGRAM  UILL  TAKE  PARAMETERS  OF  AN  ECHOSOUNDER  AND 
ESTIMATE  THE  RANGE  AS  A  FUNCTION  OF  BACKGROUND  NOISE 
THE  FOLLOWING  INPUTS  ARE  REQUIRED: 
ATMOSPHERIC  DATA 

1 )ATOMOSPHERIC  PRESSURE  in  Millibars 

Z)SELECT  ONE  OF  FOUR  ATMOSPHERIC  PROFILES  FOR  THE  TEMPERATURE 
STRUCTURE  PROFILE. 
a)FOR  ONE  PROFILE  THE  HEIGHT  OF  THE  INVERSION  LAYER  IS  ALSO  INPUT 
3)TEMPERATURE  IN  DEGREES  CELSIUS 
4)UATER  VAPOR  PRESSURE  IN  Millibars 

ECHOSOUNDER  DATA 

5)FREOUENCY 

G) ANTENNA  DIAMETER  IN  METERS 

7)PULSE  LENGTH  OF  THE  TRANSMITTED  ACOUSTIC  ENERGY 

8)P0UER  TRANSMITTED 

THE  PROGRAM  OUTPUTS  THE  FOLLOUING  GRAPH  TO  AN  EXTERNAL  PLOTTER: 
1 )RANGE(m. )  VERSUS  BACKGROUND  NOISE  IN  UATTS 

#**♦♦ VARIABLES***** 


Again      TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

WISHES  TO  MAKE  ANOUTHER  RUN. 
Ano_change  TEST  VALUE  IN  UHILE  LOOP  TO  DETERMINE  IF  OPERATOR 

UISHES  TO  CHANGE  THE  VALUE  OF  INPUT  DATA  BEFORE  MAKING 

ANOUTHER  RUN. 
Ant_area   ECHOSOUNDER  ANTENNA  AREA  CALCULATED  FROM  INPUT  OF 

ANTENNA  DIAMETER. 
Ant_dian    INPUT  OF  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS. 
Atoi«i_pres   INPUT  OF  ATMOSPHERIC  PRESSURE  IN  nb. 
Atten      ATTENUATION  OF  ACOUSTIC  WAVE.   CALCULATED  IN 

SUBPROGRAM  ATTENUATION. 
C  VELOCITY  OF  SOUND  CALCULATION  AT  ZERO  CELSIUS. 

CI         VARIABLE  USED  TO  SCALE  INVERSION  HEIGHT  IN 

TEMPERATURE  STRUCTURE  PROFILE  THREE. 
C3         VARIABLE  USED  IN  CALCULATING  THE  FOURTH  CteZ  PROFILE. 
CneZ<»)     ARRAY  OF  VALUES  OF  THE  ACOUSTIC  REFRACTIVE  INDEX 

PARAMETER.   CALCULATED  BASED  ON  SELECTION  OF 

PROFILE  FOR  TEMPERATURE  STRUCTURE  PARAMETER 

AND  ASSUMED  PROFILE  FOR  THE  VELOCITY  STRUCTURE 
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PARAMETER. 
CteZ       VALUE  OF  THE  TEMPERATURE  STRUCTURE  PARAMETER. 

THE  OPERATOR  SELECTS  A  PROFILE  FROM  SEVERAL  PROVIDED. 
CveZ       VALUE  OF  THE  VELOCITY  STRUCTURE   PARAMETER. 

VALUES  BASED  ON  CALCULATION  USING  ASSUMED  DISSAPATION 

RATE. 
Dx         STEP  SIZE  FOR  X  AXIS  FOR  VARIOUS  PLOTS. 
Dy         STEP  SIZE  FOR  Y  AXIS  FOR  VARIOUS  PLOTS. 
Epsilon    DISSAPATION  RATE  USED  IN  CALCULATION  OF  CveZ. 
Eff        TRANSMISSION  EFFICIENCY  OF  ECHOSOUNDER. 
Eff         EFFICIENCY  OF  ECHOSOUNDER  UHEN  RECEIVING  ACOUSTIC 

BACKSCATTER. 
Exc_att     EXCESS  "ATTENUATION"  AT  GIVEN  RANGE. 
F  VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 

RATIO  OF  FREQUENCY  TO  FREQUENCY  AT  MAXIMUM  ATTENUATION. 
Fmax       FREQUENCY  OF  MAXIMUM  ATTENUATION.  USED  IN  SUBPROGRAM 

ATTENUATION. 
Freq       INPUT  FREQUENCY  OF  ECHOSOUNDER. 
G  ANTENNA  EFFECTIVE  APERATURE  FACTOR. 

H  VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS  ATTENUATION. 

I  MAIN  INDEX  FOR  LOOP  TO  CHECK  RETURNED  POUER  UNTIL 

POWER  WAS  LESS  THAN  MINIMUM  DETECTABLE  SIGNAL. 
IND        INDEX  USED  FOR  LOOP  FOR  DIFFERENT  NOISE  CUT  OFF  LEVELS 
Interval    MAXIMUM  EFFECTIVE-SCATTERING  VOLUME  THICKNESS. 
Inver      HEIGHT  OF  INVERSION  LAYER. 
J  FIRST  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

K  UAVENUMBER 

L  SECOND  ORDER  INDEX  FOR  ASSORTED  LOOPS. 

Mess_up    TEST  VALUE  TO  CHECK  RESPONSES  UHEN  OPERATOR  IS  INPUTTING 

RESPONSES. 
N  TEST  VALUE  IN  CALCULATION  OF  COHERENCE  LENTH  IN  CALCULATION 

OF  EXCESS  ATTENUATION. 
Neu_va     TEST  VALUE  IN  CHECKING  IF  OPERATOR  WISHES  TO  CHANGE 

A  VARIABLE  BEFORE  A  NEW  RUN. 
Noise      ASSUMED  MINIMUM  DETECTABLE  SIGNAL. 
Pou_back(*>  POWER  BACKSCATTERED  FROM  GIVEN  RANGE. 
Pou_ret<*)  POWER  BACKSCATTERED  TO  ECHOSOUNDER  FROM  GIVEN  RANGE. 
Pou_trans   INPUT  OF  POUER  SUPPLIED  TO  TRANSDUCER  OF  ECHOSOUNDER. 
Profile    OPERATOR  INPUT  OF  CteZ  PROFILE  FROM  AVAILABLE  PROFILES. 
Pstar      VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION  FOR  INTERMEDIATE 

CALCULATION. 
Pulse      TRANSMITTED  PULSE  LENGTH  IN  MILLISECONDS.   INPUT  BY 

OPERATOR. 
R  RANGE  FOR  INPUT  TO  INTEGRAL  OF  EXCESS  ATTENUATION. 

R1         VARIABLE  USED  IN  SCALING  INVERSION  HEIGHT  FOR  PROFILE  THREE 
Range<*)    ARRAY  OF  RANGE  VALUES. 
Remainder   REMAINDER  OF  MODULO  FUNCTION  USED  TO  DECREES  THE 

NUMBER  OF  PASSES  THROUGH  THE  INTEGRATION  FOR  EXCESS 

ATTENUATION. 
Rge        VALUE  OF  RANGE  IN  INTEGRAL  FOR  EXCESS  ATTENUATION. 
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Rho        CORRELATION  LENGTH  USED  IN  CALCULATION  OF  EXCESS 

ATTENUATION. 
Sigma      FRACTION  OF  ENERGY  BACKSCATTERED  FROM  GIVEN  RANGE. 
Speed_sound  SPEED  OF  SOUND  IN  AIR  AT  INPUT  TEMPERATURE. 
T  INPUT  TEMPERATURE  IN  DEGREES  KELVIN. 

Temp       INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
5umpou_bacl< 
Temp$  "     VARIABLE  STRING  USED  IN  FUNCTION  YES. 

Tstar      INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
Var        USED  TO  SELECT  THE  VARIABLE  THE  OPERATOR  UISHED  TO 

CHANGE  BEFORE  MAKING  ANOTHER  RUN. 
Wat_pres    ATMOSPHERIC  WATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 

OPERATOR. 
X  THIRD  ORDER  INDEX  USED  IN  VARIOUS  LOOPS. 

X$         STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 

RESPONSE  TO  YES  OR  NO  QUESTION. 
Xnax       VALUE  OF  LARGEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xnin       VALUE  OF  SMALLEST  VALUE  OF  X  FOR  VARIOUS  PLOTS. 
Xrange     VALUE  OF  RANGE  OF  X  VALUES  FOR  VARIOUS  PLOTS. 
YlabelS    LABEL  ON  Y  AXIS  PASSED  TO  SUBPROGRAM  Semi_log  FOR  PLOTTING. 
Ymax       LARGEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Ynin       SMALLEST  VALUE  OF  Y  VARIABLE  FOR  VARIOUS  PLOTS. 
Yrange     RANGE  OF  Y  VARIABLES  FOR  VARIOUS  PLOTS. 
Ze         INTERMEDIATE  VALUE  IN  CALCULATION  OF  EXCESS  ATTENUATION. 


DECLARE  VARIABLES 
INTEGER  I 

DIM  Pou_back<5000) ,Pou_ret(  5000) ,Range(5000) 
DIM   CneZ(5000),Ran<50) 


PLOTTER  IS  70S,"HPGL" 

LINE  TYPE  1 

! 

(INPUT  ATMOSPHERIC  DATA 

INPUT  "ENTER  TEMPERATURE  IN  DEGREES  CELCIUS" , Temp 

INPUT  "ENTER  ATMOSPHERIC  PRESSURE  IN  MILLIBARS"  ,  At oivi_pres 

INPUT  "ENTER  WATER  VAPOR  PRESSURE  IN  MILLIBARS" , Wat_pres 


INPUT  ECHOSOUNDER  DATA 
INPUT  "ENTER  ECHOSOUNDER  FREQUENCY  IN  HERTZ", Freq 
INPUT  "ENTER  ECHOSOUNDER  PULSE  LENGTH  IN  MILLISECONDS=" .Pulse 
INPUT  "ENTER  THE  ECHOSOUNDER  ANTENNA  DIAMETER  IN  METERS" .Ant.diam 
INPUT  "ENTER  THE  POWER  TRANSMITTED  IN  WATTS" ,Pou_trans 
G=.40        (ANTENNA  EFFECTIVE  APERATURE  FACTOR 
Eff=.2S 
I 
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! 

! SELECT  THE  PROFILE  FOR  THE  TEMPERATURE  STRUCTURE 
PRINT  "YOU  HAVE  A  CHOICE  OF  TEMPERATURE  STRUCTURE  PROFILES- 
PRINT  "PROFILE  TO  BE  USED." 
PRINT  "YOUR  SELECTIONS  ARE" 

PRINT  "      1      A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS  " 
PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  UHICH  " 

PRINT  "  GIVES  A  HEIGHT  TO  THE  -1.1G  PROFILE" 

PRINT  "  " 

PRINT  "      Z      THE  TEMPERATURE  STRUCTURE  PROFILE  ABOVE  BUT" 
PRINT  "  UITH  A  FACTOR  OF  TUO  TO  APPROXIMATE  LOOKING" 

.PRINT  "  UP  A  THERMAL  PLOOM" 

PRINT  "  " 

PRINT  "      3      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS- 
PRINT  "  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
PRINT  "  NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
PRINT  "  IS  EXP(-.001*HEIGHT  ABOVE  G5  METERS  AND  HEIGHT" 
PRINT  "            TO  THE  -1.4B  BELOU  65  METERS" 
PRINT  "  " 

PRINT  "     4      A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A" 
PRINT  "  HEIGHT  TO  THE  -4/3  " 

! 

PRINT  "  " 
PRINT  "  " 
Mess_up=1 
WHILE  Mess_up=1 

INPUT  "ENTER  THE  DESIRED  PROFILE  (1  OR  Z  OR  3  OR  4)" .Profile 
IF  Profile=4  THEN 

INPUT  "ENTER  THE  HEIGHT  OF  THE  INVERSION  IN  METERS"  , Inver 
END  IF 
IF  Profiled  THEN 

Mess_up=0 
ELSE 

IF  Profile=Z  THEN 

Mess_up=@ 
ELSE 

IF  Profiled  THEN 

Mess_up=0 
ELSE 

IF  Profile=4  THEN 

Mess_up=0 
ELSE 

PRINT  Profile,"  UAS  NOT  ONE  OF  THE  OPTIONS!!!! 
Mess_up= 1 
END  IF 
END  IF 
END  IF 
END  IF 
END  UHILE 
OUTPUT  KBD;"K" ; 
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! 
I 

Again=1 
UHILE  Again=1 
FOR  J=1  TO  I 

Pou_back< J)=0 

Pou_re-t(  J)=0 

Range< J)=0 

Cne2(J)=0 
NEXT  J 
! 

(CONVERT  TEMPERATURE  TO  KELVIN 
T=Tenp+273 
I 

! CALCULATE  SPEED  OF  SOUND  BASED  ON  INPUT  TEMPERATURE* CELCIUS ) 
Speed_sound=20. 05* (  T ) A . 5 

! CALCULATE  THE  SPEED  OF  SOUND  AT  0  DEGREES  CELCIUS 
C  =  20. 05*273". 5 
! 

1=1     ! INDEX  FOR  INNER  LOOP  ( POUER  RETURNED  TO  NOISE) 
Ind-1   ! INDEX  FOR  SECOND  LOOP  (NOISE) 

Pou_back<0>=0     'INITIALIZE  VARIABLE  FOR  POUER  BACKSCATERRED 
Interval =( Speed_sound*Pulse* 1 . E-3)/Z 
Ant_area=PI*< Ant_diam/2 )"2 
! 

(CALCULATE  RETURNED  SIGNAL  POUER  UNTIL  IT  IS  LESS  THAN  BACKGROUND 

IFOR  EACH  FREQUENCY. 

Range(0)=0 

Noise=1 .E-18 

FOR  Ind=1  TO  5 

CALL  Attenuat  ion(  Atorn_pres  ,  Atten  ,  Freq  ,  Temp  ,  Uat_pres  ) 
K=2*PI*Freq/Speed_sound     (UAVENUMBER 
1  =  1 

REPEAT 
Range< I )=Range(  1-1  >+2 
Sunpou_back  =0 
SELECT  Profile 
CASE  1 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
!UAS  FOR  MIDDAY  AND  CLEAR  WEATHER  ABOVE  THE 
ITULAROSA  BASIN,  NEW  MEXICO. 

!AN  ADDITONAL  FACTOR  OF  1 . 8  UAS  INCLUDED  AS  PER  DIFFERENT 
(AVERAGING  TIME. 

Cte2=2. 12*Range(I )*(-1 . IB) 
CASE  2 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  IS  BASED  UPON  DATA 
!AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398. 
IUAS  FOR  MIDDAY  AND  CLEAR  UEATHER  ABOVE  THE 
ITULAROSA  BASIN,  NEU  MEXICO. 
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(AN  ADDITONAL  FACTOR  OF  1.8  UfiS  INCLUDED  AS  PER  DIFFERENT 
! AVERAGING  TIME  AND  AN  ADDITIONAL  FACTOR  OF  TWO  UAS 
'INCLUDED  TO  APPROXIMATE  LOOKING  UP  A  THERMAL  PLOOM 

CteZ=Z*Z. lZ*Range< I )A<-1 . IB) 
CASE  3 

(THIS  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA  AS" 
(PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  FOR" 
(NIGHT  CONDITIONS.   THE  DEPENDENCE  UITH  HEIGHT  " 
(IS  EXP<-.001»HEIGHT  ABOVE  65  METERS  AND  HEIGHT" 
(TO  THE  -1.4B  BELOW  G5  METERS" 
IF  Range(I><65  THEN 

CteZ=75.5»Range< I >A<-1 .46) 
ELSE 

C-teZ  =  3.66E-Z*EXP<-.001«Range<  I  )  ) 
END  IF 
CASE  4 

(CALCULATE  THE  TEMPERATURE  STRUCTURE  FACTOR  BASED  ON  ASSUMPTION 
!     THAT  IT'S  PROPORTIONAL  TO  HEIGHT  TO  THE  -4/3  POUER 
!     THAT  THE  SURFACE  VERTICAL  HEAT  FLUX  IS  .095  Cm. /sec 
(THE  DEPENDENCE  UITH  HEIGHT  ABOVE  THE  INVERSION  I  PULLED  FROM 
(MY  LEFT  EAR 

(EQUATION  FROM  NEFF ,  1975 
SELECT  Ranged  )/Inver 
CASE  <.9 

C3=(< .0Z4)*(T)A( .667)) 
CteZ=C3*<  Ranged  >)A<-1.33) 
CASE  .9  TO  1 
CteZ=CteZ 
C1=CteZ 
R1=Range( I ) 
CASE  1  TO  1.3 

CteZ=10A( (LGT(C3*Range(  I ) )-LGT(C3*R1 ) )+LGT(C1  ) ) 
RZ=Range( I ) 
CZ=CteZ 
CASE  ELSE 

CteZ=C3*(Range< I ) A( -  1 . 33>-RZA(  -  1 . 33 ) )  +  CZ 
END  SELECT 
END  SELECT 

(CALCULATE  VELOCITY  STRUCTURE  FACTOR.  FORMULA  FROM  GAYNOR  77 
Epsilon=( .Z866/Range( I ) )»(  1 + .07*< Range( I ) ) A(  . 6 ) ) " ( 1 . 5 ) 
CveZ  =  Z*EpsilonA(  .667) 

(CALCULATE  ACOUSTIC  REFRACTIVE  INDEX  STRUCTURE  FACTOR. 
(FORMULA  FROM  TETARSKI 
CneZ( I )  =  (CteZ/(Z.98E+5) )  +  (  CveZ/(  C»C ) ) 

(CALCULATE  THE  FRACTION  OF  POUER  BACKSCATTERED  FROM  INTERVAL. 
(EQUATION  FROM  NEFF,  1975 
Sigtia=(  .0039*<K"(  1 /3  )  )«CteZ  )/(  T)  AZ 
I 

(CALCULATE  THE  EXCESS  ATTENUATION.   EQUATION  FROM  CLIFFORD,  1980 
(THE  EXCESS  ATTENUATION  IS  Exc  ati 
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IF   K50  THEN 

Remainder =0 
ELSE 

IF    K150  THEN 

Remainder=Range(  I )  MODULO  10 
ELSE 

Remainder=Range( I )  MODULO  20 
END  IF 
END  IF 

IF  Remainder=0  THEN 
Rho=0 
L=0 
H=0 

Rge=Range(I)     ! CONSTANT  IN  INTEGRAL 
R=0 
FOR  J=0  TO  2*1 

F=Cne2(INT( J/2+1 )) 

F=F*< 1-R/(Rge)A< 1 . G7)  +  (  R/(  Rge) ) "<  1 . 67 ) ) 

IF  J>0  THEN 

IF  J<2*I  THEN 

IF  INT(J/2)=J/2  THEN 
L  =  L  +  F 
F=0 
ELSE 

H=H+F 
F=0 
END  IF 
END  IF 
END  IF 
Rho=Rho+F 
R=R+1 
NEXT  J 

Rho=Rho+4*L+2*H 

Rho=(((Rho*.33)*K*K*1 .4G)A(-.G) ) 
N=(Ant_diam/Rho)"2 
IF  N<=1  THEN 

2e=1/( 1+N) 
ELSE 

2e=1 ,5/( 1+N) 

!STEP  OF  1.5==>SEE  CLIFFORD  1980 
END  IF 

Exc_att=2e*2e 
ELSE 

Exc_att=Exc_att 
END  IF 
! 

(CALCULATE  THE  POUER  BACKSCATTERED 

Pou_back< I )=(Pou_trans*Ef f-Sumpou_back )*EXP( -Atten*Range( I ) ) 
Pou_back  (  I  )=Pou_back(  I  )*Interval*Exc_att*Sigi»ia 
Sunpou_back=Sunpou_back+Pou_back ( I ) 
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I  CALCULATE  THE  POUER  BACKSCATTERED  TO  THE  ANTENNA 
Pou_ret(  I  >-Pou_back<  I  )*EXP( -Atten*Range(  I  )  )»Ant_area*G*Ef f  /Range<  I  )*2 
PRINT  "RANGE-" , Range* I > 
PRINT  "POUER  RETURNED-" ,Pou_ret< I ) 
PRINT  "NOISE-" .Noise 
1  =  1  +  1 
UNTIL  Pou_ret(I-1 )  >0  AND  Pou_ret(  I- 1 )<  =  Noise 
Ran( Ind)-Range(  I-Z ) 
Again- 1 
FOR  J=0  TO  I 

Pou_back(J>=0 
Pou_ret(J)=0 
Range(  J>=0 
CneZ< J)=0 
NEXT  J 

Noise-Noise* 10 
NEXT  Ind 
) 
! 

OUTPUT  KBD;"K"; 
PRINT  "INPUT  CONDITIONS" 


PRINT  "  " 

PRINT  USING  "K"; "TEMPERATURE- 
PRINT  USING  "K" 


PRINT  USING  "K" 
PRINT  USING  "K" 
PRINT  USING  "K" 


"ATMOSPHERIC  PRESSURE- 
PRINT  USING  "K'V'UATER  VAPOR  PRESSURE- 
PRINT  USING  "K"; "FREQUENCY- 

"PULSE  LENGTH- 
"TRANSDUCER  EFFICIENCY 
"POUER  TRANSMITTED- 
PRINT  "TEMPERATURE  STRUCTURE  PROFILE  USED  ".Profile 
IF  Profile-4  THEN 

PRINT  "INVERSION  HEIGHT  ".Inver 
END  IF 

SELECT  Profile 
CASE  1 

PRINT  "      1 
PRINT  " 
PRINT  " 
CASE  Z 

PRINT  "  Z 
PRINT  " 
PRINT  " 
PRINT  " 
PRINT  " 
CASE  3 

PRINT  "      3 
PRINT  " 
PRINT  " 
PRINT  " 


,Tenp,"CELCIUS" 

, Aton_pres , "nb" 

, Uat_pres , "mb" 

,Freq,"  Hz" 

.Pulse , "ns" 

,Eff 

,Pou_-trans,"  UATTS" 


A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398  " 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.1E  PROFILE" 

A  TEMPERARURE  STRUCTURE  PROFILE  BASED  ON  DATA" 
AS  PRESENTED  IN  WALTERS/KUNDEL  1981  PAGE  398  " 
UHICH  GIVES  A  HEIGHT  TO  THE  -1.16  PROFILE" 
BUT  WITH  A  FACTOR  OF  TUO  TO  APPROXIMATE- 
LOOKING  UP  A  THERMAL  PLOOM" 

A  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  DATA" 
AS  PRESENTED  IN  UALTERS/KUNDEL  1981  PAGE  398" 
FOR  NIGHT  CONDITIONS.   THE  DEPENDENCE  WITH" 
HEIGHT  IS  EXP<-.001*HEIGHT  ABOVE  GS  METERS" 
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PRINT  "  AND  HEIGHT  TO  THE  -1.4E  BELOU  E5  METERS" 

CASE  4 

PRINT  "4ft  TEMPERATURE  STRUCTURE  PROFILE  BASED  ON  A' 

PRINT  "  HEIGHT  TO  THE  -4/3  " 

END  SELECT 
I 
! 

PRINT  "  " 
I 
! 

PRINT    "HIT   CONTINUE    TO   CONTINUE" 
PAUSE 

OUTPUT   KBD;"K"; 
! 

GRAPHICS    ON 

VIEWPORT    15, 1Z0, 10,70 

!PLOT    RANGE    VERSUS    NOISE 

! 

Xnin=-18 

Xmax=- 1 3 

Yniax=4 

Ymin=Z 

Dx  =  1 

Dy=1 

Yrange=Z 

Xrange=5 

UINDOU    Xmin,Xmax,Ymin,Yiiax 

AXES  Dx,Dy,Xmn,Ynin,  1  ,  1 

CLIP  OFF 

! 

(LABEL  PLOT 

CSIZE  4,.E 

LDIR  0 

LORG  4 

MOVE  Xmin+ . 5*Xrange , Ynin+ 1 *Yrange 

LABEL  "RANGE  vs  NOISE" 

! LABEL  HORIZONTAL  AXES 

LDIR  0 

LORG  5 

FOR  J=Xmin  TO  Xmax  STEP  Dx 

CSIZE  4, .6 

MOVE  J- . 01 3*Xrange , Ymi n- . 0G7* Yrange 

LABEL  "10" 

MOVE  J+.01 3*Xrange, Ynin- .033* Yrange 

CSIZE  Z 

LABEL  J 
NEXT  J 

MOVE  Xmin+.5*Xrange,Ynin-. 1Z*Yrange 
LORG  5 
CSIZE  4, .6 
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LABEL  "NOISE   (UATTS)" 
•LABEL  VERTICAL  AXES 
LORG  8 

FOR  J=Ymin  TO  Ymax  STEP  Dy 
CSIZE  4,.G 

MOVE  Xnin-.01*Xrange,J 
LABEL  10 
CSIZE  Z 

MOVE  Xmin- . 00Z5* Xrange , J+ . 03» Yrange 
LABEL  J 
NEXT  J 
LDIR  PI/Z 
LORG  G 

MOVE  Xnin- . 1  * Xrange ,Ymin+. 5* Yrange 
CSIZE  4, .6 

LABEL  "RANGE   (METERS)" 
CLIP  ON 
Ind=1 

Noise=1 .E-18 
FOR  Ind=1  TO  5 

PLOT  LGT(Noise),LGT(Ran(Ind)> 
Noise=Noise* 10 
NEXT  Ind 
I 

PRINT  "HIT  CONTINUE  TO  CONTINUE" 
PAUSE 
GCLEAR 

OUTPUT  KBD;"K"; 
! 
! 

LINPUT  "UOULD  YOU  LIKE  TO  MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)?" ,X$ 
Again=FNYes(X$) 
Mess_up=1 
UHILE  Mess_up=1 
SELECT  Again 
CASE  1 

LINPUT  "UOULD  YOU  LIKE  TO  CHANGE  A  VARIABLE  (ENTER  Y  OR  N>?" ,XS 
Neu_va=FNYes( X$) 
Ano_change= 1 
UHILE  Ano_change=1 
SELECT  Neu_va 
CASE  1 
PRINT  "       VARIABLE  CURRENT  VALUE" 

PRINT  USING  "K";  "1   TEMPERATURE        " , Temp , "CELSIUS" 
PRINT  USING  "K'V'Z   ATOMOSPHERIC  PRESSURE  "  ,  Aton_pres  ,  "mb' 
PRINT  USING  "K";"3   UATER  VAPOR  PRESSURE  " ,Uat_pres , "mb" 
PRINT  "  " 

PRINT  USING  "K";"4   ANTENNA  DIAMETER      " , Ant_diam , "  n." 
PRINT  USING  "K";"5   PULSE  LENGTH  ".Pulse,"  rns" 

PRINT  USING  "K";"G   POUER  TRANSMITTED  "  ,Poui_trans,"  UATTS' 
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PRINT  USING  "K";"7   ATMOSPHERIC  PROFILE   " . Profile 

PRINT  "  " 

PRINT  "ENTER  THE  NUMBER  OF  THE  VARIABLE  YOU  WISH  TO" 

PRINT  "CHANGE" 

INPUT  Var 

SELECT  Var 

CASE  1 

INPUT  "TEMPERATURE  =  "  ,Tei*ip 
CASE  Z 

INPUT  "ATMOSPHERIC  PRESSURE  IN  nb=" , Aton_pres 
CASE  3 

INPUT  "UATER  VAPOR  PRESSURE  IN  nb=" , Uat_pres 
CASE  4 

INPUT  "ANTENNA  DIAMETER  IN  m. =" , Arvt_dian 
CASE  5 

INPUT  "PULSE  LENGTH  IN  ms=" , Pulse 
CASE  G 

INPUT  "POUER  TRANSMITTED  IN  UATTS="  ,Pou_"trans 
CASE  7 

PRINT  "  1   ==>  C-TZ  PROFILE  OF  Z*<-1.16->" 
PRINT  "        FROM  UALTERS/KUNDEL  1981" 
PRINT  "  Z   ==>  SAME  AS  ONE  BUT  KITH  FACTOR" 
PRINT  "        OF  TUO  TO  APPROXIMATE  LOOKING" 
PRINT  "        UP  A  THERMAL  PLOOM" 
PRINT  "  3  ==>  CtAZ  PROFILE  OF  EXP< - .001 *Z )" 
PRINT  "        FROM  UALTERS/KUNDEL  1981." 
PRINT  "  4   ==>  CT'Z  PROFILE  OF  ZA(-4/3)" 
INPUT  "ENTER  NUMBER  OF  DESIRED  PROFILE" , Prof ile 
IF  Profile=4  THEN 

INPUT  "HEIGHT  OF  INVERSION  IN  METERS=" , Inver 
END  IF 
CASE  ELSE 

PRINT  Var. "IS  NOT  ONE  OF  THE  OPTIONS- 
END  SELECT 

LINPUT  "MADE  ANOTHER  CHANGE( ENTER  Y  OR  N)?",X$ 
Neu_va=FNYes(X$) 
Mess_up=Z 
CASE  Z 

Mess_up=Z 
Ano_change=Z 
CASE  ELSE 

PRINT  "YOUR  CHOICES  UHERE  Y  OR  N  MM!" 
Ano_change=1 

LINPUT  "CHANGE  A  VARIABLE  (ENTER  Y  OR  N  )?",><$ 
Neu_va=FNYes(X$) 
END  SELECT 
END  WHILE 
CASE  Z 

Mess_up=Z 
CASE  ELSE 
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PRINT  "YOUR  CHOICES  UHERE  Y  OR  N   I ■ !  M  " " 

Mess_up= 1 

LINPUT  "MAKE  ANOTHER  RUN  (ENTER  Y  OR  N)7",X$ 

Again=FNYes(X$> 
END  SELECT 
END  WHILE 
END  UHILE 

PRINT  "THAT'S  ALL,  FOLKS" 
end 
I 
I 

ICALCULTE  THE  ATTENUATION 

sub  Attenuat  ion( Atom_pres , At  ten ,Freq , Temp ,Uat_pres ) 
I  THIS  SUBPROGRAM  CALCULATES  THE  ATTENUATION  OF  SOUND 
!IN  AIR  BASED  UPON  EQUATIONS  IN  NEFF  1975 
I 

!    INPUT    ATMOSPHERIC  PRESSURE  IN  MILLIBARS 
I  FREQUENCY  OF  SOUND  UAVE  IN  HERTZ 

I  TEMPERATURE  IN  DEGREES  CELCIUS 

!  UATER-VAPOR  PRESSURE  IN  MILLIBARS 
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OUTPUT   ATTENUATION  IN  1 /METERS 


i VARIABLES 
!   Atom_pres 
I   Atten 
I 

I   Att  max 


Fnax 

Freq 

H 

Pstar 

Temp 
Tstar 

Uat_pres 


INPUT  OF  ATMOSPHERIC  PRESSURE  IN  mb. 

ATTENUATION  OF  ACOUSTIC  UAVE.   CALCULATED  IN 
SUBPROGRAM  ATTENUATION. 

VARIABLE  IN  SUBPROGRAM  ATTENUATION.   IT  IS  THE 
ATTENUATION  AT  THE  FREQUENCY  OF  THE  MAXIMUM 
ATTENUATION  FOR  THE  INPUT  CONDITIONS. 

VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION.   IS  THE 


RATIO  OF  FREQUENCY  TO  FREQUENCY 
FREQUENCY  OF  MAXIMUM  ATTENUATION 

ATTENUATION. 
INPUT  FREQUENCY  OF  ECHOSOUNDER. 
VARIABLE  USED  IN  INTEGRATION  FOR  EXCESS 


AT  MAXIMUM  ATTENUATION. 
USED  IN  SUBPROGRAM 


ATTENUATION. 
FOR  INTERMEDIATE 


VARIABLE  USED  IN  SUBPROGRAM  ATTENUATION 

CALCULATION. 
INPUT  TEMPERATURE  IN  DEGREES  CELSIUS. 
INTERMEDIATE  VALUE  USED  IN  CALCULATION  OF  ATTENUATION 

IN  SUBPROGRAM  ATTENUATION. 
ATMOSPHERIC  UATER  PRESSURE  IN  MILLIBARS.   INPUT  BY 


I  OPERATOR. 

H=100*Uat_pres/Atom_pres 

Tstar=( 1 . 8«TempM49Z)/519 

Pstar=Atom_pres/1014 

Fmax  =  (  10-tG600«H444400»H«H)«Pstar/Tstar*  .8 

Att_max=.0078*Fmax»Tstar*<-Z.5)»EXP(7.77»<  1-1 /Tstar  )  ) 

F=Freq/Fmax 

Atten=( Att    max/304 . 8 )»( <  . 1 8»F  )" Z  +  ( Z*F»F/(  1 +F»F  )  r Z r  . ! 


53 


Atten=(Atten-H  .74E- 10*Freq*Freq)/4. 35 
SUBEND 
! 

def  FNYes(X$) 
THIS  FUNCTION  INTERPRETS  THE  OPERATORS  RESPONSES  TO  YES  NO  QUESTIONS 


INPUT 


OUTPUT 
VARIABLES 
Tei»ip$ 
X$ 


X$ 


FNYes 


VARIABLE  STRING  USED  IN  FUNCTION  YES. 

STRING  PASSED  TO  FUNCTION  YES  TO  DETERMINE  OPERATOR 
RESPONSE  TO  YES  OR  NO  QUESTION. 
DIM  Temp$[ U 
Temp$[ 1 ,1]=TRIM$(X$) 
SELECT  Temp$ 
CASE  "Y","y" 


RETURN 
CASE  "N"  ,"n' 

RETURN 
CASE  "  " 

RETURN 
CASE  ELSE 

RETURN 
END  SELECT 


FNEND 
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